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Quadrupole  perturbed  nuclear  magnetic  resonance  (NMR)  has  been  used 

2 
to  determine  the  quadrupole  coupling  constant  e'qQ/h,  asymmetry  para- 
meter n,  and  electric  field  gradient  (EFG)  tensor  in  single  crystal 

2 
glycine.   The  results  at  room  temperature  are  e  qQ/h  =  1190  kHz,  0  =  0.505, 

with  the  principal  axis  roughly  parallel  to  the  C-N  bond  direction.   As 

2 
the  temperature  increased,  e  qQ/h  decreased.   These  results  are  consis- 
tent with  earlier  work,  except  there  is  no  evidence  of  a  significant 
change  in  the  orientation  of  the  EFG  tc-nsor  at  higher  temperatures  as 

reported  earlier. 

15  15 

The  NMR  spectrum  of  "'  N  in  o-phase  solid   N,  with  paramagnetic  0^ 

impurity  added  has  been  d^terniined  as  a  function  of  magnetic  field  and 
0-  concentration.   A  classical  Pake  doublet  was  obtained  at  650  G  as 
expected;  however,  <;he  spectrum  became  increasingly  asymmetric  as  the 
field  was  increased  to  a  maximi'm  of  9  kG.   Increasing  the  0.,   concentra- 
tion broadened  Ihe  line  li\A.   did  not  affect  the  asymmetry.   Various 


mechanisms  may  affect  the  line  shape  but  none  have  been  found  to  satis- 
factorily explain  the  asymmetry  in  all  quantitative  aspects. 

The  tJMR  spectrum  of   N  in  an  a-phase  solid  containing  isotopic 

15   15   14   14 
molecular  species  in  the  ratio   N  :   N-  N:   N  =  1:4:4  has  been  obtain- 
ed as  a  function  of  magnetic  field  and  compared  with  computer  simulations. 

Besides  the  Pake  doublet  due  to  the   N  molecules,  an  asymmetric  triplet 

15   14  14 

results  from  the   N-   N  mixed  molecule.   In  the  latter  case  the   N 

nucleus,  which  e::pericnces  a  quadrupolar  interaction  comparable  with  the 

15 
Zeeman  Harailtonian,  produces  a  dipolar  interaction  at  the   N  nucleus 

which  varies  with  the  magnetic  field  in  accordance  with  theoretical  ex- 
pectations.  Additional  effects  occur  due  to  level  crossing  and  to  the 
asymmetry  of  the  Pake  doublet  mentioned  above. 

The  nuclear  spin-lattice  relaxation  time  T  of   N  has  been  measured 
in  liquid   N_  for  the  temperature  interval  63  <  T  <  77  K  and  in  the  6- 
phase  solid  for  38  <  T  <  63  K.   The  relaxation  is  attributed  to  the 
spin-rotational  interaction  in  both  the  liquid  and  6  solid.   Any  possible 

discontinuity  in  T  across  the  triple  point  is  obscured  by  scatter  in  the 

14 
data.   When  combined  with  previous  T^  data  for   N  ,  it  is  concluded  that 

a  rotational  diffusion  nod?.?  developed  for  liquids  may  hold  in  the  B 

solid,  and  that  the  -.nolecuies  reorif:nt  r.ore  slowly  in  the  liquid  than  in 

the  solid  near  the  triple  point. 


CHAPTER  I 
INTRODUCTION 


1.1   High  Field  NQR  of   N  in  Single  Crystal  Glycine 


Nitrogen  is  an  important  atom  in  biological  material,  and  nuclear 

14 
quadrupole  resonance  of   N  nuclei  in  biological  molecules  can  provide 

information  about  the  structure  of  the  molecules.   Problems  arise  due 

to  the  usually  low  atomic  density  and  the  low  resonant  frequencies 

which  combine  to  give  a  poor  signal-to-noise  ratio.   Glycine,  for  ex— 

lA 
ample,  has  a   N  nuclear  spin  density  of  one  nucleus  per  five  C,  0, 

y 

and  N  atoms,  not  counting  the  five  hydrogen  atoms.   Furthermore,  there 

are  two  inequivalent  sites,  reducing  the  density  of  those  partici- 

14 
pating  in  a  resonance  by  a  factor  of  two.   The   N  quadrupole  coupling 

constant  is  about  1200  kHz  at  room  temperature  and  increases  only 
slightly  at  77  K  which,  combined  with  the  low  density,  results  in  a 
weak  signal  that  is  difficult  to  observe.   To  date,  the  pure  quadru- 
pole resonance  has  not  been  observed  to  the  knowledge  of  the  author. 

Three  research  groups,  Andersson,  Gourdji,  Guibe,  and  Proctor 
[l],  Blinc,  Mali,  Osredkar,  Prelesnik,  Zupancic,  and  Ehrenberg  [2], 

and  Edmonds  and  Speight  [3],  have  reported  the  quadrupole  coupling 

">  14 

constant   e*'qQ/h  and  asymmetry  parameter  q  of   N  in  glycine 

obtained  from  quadrupole  perturbed  ISTMR.  References  [l]  and  [2]  also 
reported  the  quadrupole  coupling  tensor,  or  eouivalently,  the  direc- 
tion cosines  of  the  principal  axes  of  the  ti;nsor.   Edrionds  and  Speight 


.1- 


-2- 


[3]  used  a  piowder  sample  and  therefore  were  unable  to  obtain  the  com- 
plete tensor. 

Andersson  £t  al.  [l]  used  a  wide-band  Varian  spectrometer  and 
were  able  to  obtain  enough  data  to  deduce  the  complete  quadrupole 
coupling  or  electric  field  gradient  (EFG)  tensor  at  room  temperature. 
They  reported  that  the  principal  axis  was  roughly  parallel  to  the 
axis  of  the  NK  group.   Blinc  e^  al-  [2]  used  pulsed  double  resonance 
methods  at  a  temperature  of  140°  C,  and  reported  the  principal  axis 
was  about  60°  away  from  the  direction  obtained  by  Andersson  et  al. ; 
furthermore  a  much  lower  quadrupole  coupling  constant  was  deduced. 

The  purpose  of  the  present  experiment  was  two-fold.   First,  the 
laboratory  had  recently  acquired  a  superconductive  magnet  capable  of 
producing  up  to  a  100-kG  field  and  it  was  decided  to  employ  it  in  an 
investigation  of  the  use  of  a  strong  magnetic  field  to  overcome  the 
weak  signal.   Next,  an  attempt  would  be  made  to  explain  the  differences 
between  the  data  of  Andersson  et^  £l .  [l]  and  those  of  Blinc  £t  al^.  [2]. 

1.2   NHR  in  Liquid  and  Solid  Nitrogen 

Nitrogen  is  an  important  and  convenient  substance  for  NMR  studies 
for  a  number  of  reasons.   Liquid  nitrogen  as  a  cryogenic  fluid  is 
readily  available  at  low  cost  and  requires  no  special  handling  other 
than  that  due  to  its  low  temperature.   Aside  from  the  attractive 
cost  factor,  it  is  advantageous  to  have  the  sample  and  the  cryogenic 
fluid  both  the  same,  as  it  facilitates  temperature  and  pressure  con- 
trol.  A  temperature  region  from  about  50  K  to  77  K  in  solid  and 
liquid  phases  is  easily  attainable  using  liquid  nitrogen. 


-3- 


More  important  than  the  practical  aspects  mentioned  above  are 

the  physical  properties.   Two  nonradioactive  isotopes  are  available, 

14  15 

the  abundant   N  with  spin  1=1  and   N  with  I  =  h   which  has  only 

about  1/3%  natural  abundance.   Three  species  of  diatomic  molecules 

15     15   14       14 
are  therefore  available,    N  ,    N-   N,  and   N^. 

Only  three  nonradioactive  elements  have  spin  1=1,  viz.,   N, 

7        2 
Li,  and   H.   Lithium  is  a  metal  and  deuterium  is  such  a  light  atom 

that  quantum  effects  are  considerably  more  important  than  in  the 

14  14 

heavier   N.   Therefore,   N  occupies  a  r.-'ther  unique  position.   As 

nuclei  with  spin  '5  have  no  quadrupole  moment,  nitrogen  can  be  chosen 

with  or  without  quadrupole  interactions  by  selecting   N„  or   N  , 

15   14 
respectively,  and   N-   N  mixed  molecules  permit  a  study  of  unlike 

spin  coupling.   Three  very  different  spin  systems  are  therefore 
available  with  essentially  tne  same  lattice  dynamics. 

As  implied  by  the  name,  nuclear  spin-lattice  relaxation  depends 
on  the  coupling  between  nuclear  spins  and  the  lattice.   The  possi- 
bility of  selecting  either  I  =  1  or  I  =  J^  nuclei  in  essentially  the 
same  lattice  turns  out  to  be  even  more  convenient  in  spin-lattice 

relaxation  studies  than  the  simple  fact  that  the  spins  are  different. 

14 
It  happens  that  the  quadrupolar  interaction  is  dominant  in   N  nu- 
clear relaxation  and  spin-rotation  interaction  is  dominant  in   N^ 
relaxation  as   N„  experiences  no  quadrupolar  interaction.   Hence, 
two  fundamental  relaxation  processes  can  be  isolated  in  the  same 
environment.   Furthermore,  we  will  find  the  particular  two  processes 
are  uniquely  related  in  that  a  change  in  molecular  dynamics  due  to  a 
temperature  change  has  opposite  effects  on  the  two  processes. 

It  is  also  interesting  to  note  that  3  nitrogen,  stable  from  35  K 


-4- 

to  the  triple  point,  63  K,  is  a  plastic  crystal  with  low  molecular 
orientational  order.   Thus  rotational  and  orientational  dependent 
properties  in  B-solid  nitrogen  have  some  similarities  to  the  same 
properties  in  liquid  nitrogen. 

Finally,  there  are  two  crystalline  forms  readily  obtained  at  low 
pressures,  3  nitrogen  mentioned  earlier,  and  a   nitrogen,  a  much  more 
rigid  form,  stable  below  35  K.   In  a  nitrogen,  molecules  undergo 
librational  motion  about  fixed  directions  in  the  lattice.   The  effect 
of  the  librational  motion  for  some  purposes  is  to  reduce  the  magnitude 
of  the  interaction  by  motional  averaging.   The  resulting  line  shape 
is  therefore  unchanged  except  for  an  overall  reduction  in  width  due 
to  the  motional  averaging.   We  therefore  can  conduct  line  shape 
studies  of  nuclei  in  molecules  with  like  spin  (I  =  -2)  or  with  unlike 

y 

spin  (I  =  l,is)  in  the  same  crystalline  structure. 

TYie   purpose  of  the  present  v;oi"k  is  to  study  the  line  shape  of   N 
in   N   and  in  a  mixture,  of    N^  and    N--   N  samples,  and  to  study  the 
spin-lattice  relaxation  times  in  liquid  and  6-solid   X,.  and  thereby 
fill  a  gap  in  existing  nitrogen  data.   In  so  doing,  contributions  will 
also  be  made  to  the  understanding  of  fundamental  processes,  such  as 
spin-rotation  interaction  in  solids,  which  will  have  applications  to 
other  sj'stems. 


CHAPTER  ri 
THEORY 


2.1  General 


A  nucleus  with  a  magnetic  moment  p  in  a  uniform  external  magnetic 
field  H  has  an  energy  given  by 


E  =  -u-H     .  (2.1.1) 

o 


The  quantum  mechanical  Hamiltonian  is  the  well-known  Zeeman  term 


V^   =  -y-H  =  -yfil-H   =  -YtiH  I  (2.1.2) 

Z       o         o       o  z 


where  I  is  the  nuclear  spin,  y   is  the  gyromagnetic  ratio,  and 

K   =  zH  . 
o     o 


The  Zeeman  energy  levels  are 


E      =    <m  r?/^  ku>   =   -  vhH  m        , 
m  Z  ■      o 


or 


E     =  -fti.j  m        ,  (2.1.3) 

m  o 


where  o)     =  -yH   . 
o  o 


-5- 
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Transitions  between  levels  may  be  induced  by  application  of  a 

->• 
small  rotating  magnetic  field,  perpendicular  to  H^, 


H  =  H,cosa)t,  H  =  H,  sinwt   ,  (2.1.  A) 

X    1     '   y    1 


as  a  perturbation  of  H  .   In  practice,  an  alternating  magnetic  field 
is  applied  which  can  be  thought  of  as  the  sura  of  two  counter-rotating 
components.   It  will  be  seen  later  that  only  one  component  is  effec- 
tive in  inducing  transitions. 

The  perturbing  Hamiltonian  is 


-iu)t      iwt,     (2.1.5) 
.1  =  -iimH,  u.e 

y  y 


y^    =  -rri(H  I   +  H  I  )  =  -IsYftH^d  e  """  +  I  e""") 
1       XX     yy        1  + 


where  uhe  raising  and  lowering  operators   1+  ^  I  i  il  have  been 

_    X     y 

introduced.   These  operators  have  nonvanishing  matrix  elements  only 

between  states  for  which  m  differs  by  one,  so  the  only  transitions 

allowed  will  differ  in  energy  by  AE  =  ficj  .   llius,  the  frequency  of  the 

perturbing  alternating  magnetic  field  must  be  co^  in  order  to  induce 

transitions.   The  counter-rotating  component  -o)   is  2(o  off  resonance 

o       o 

and  can  be  ignored  for  most  purposes.   The  resonant  condition,  to  =  o)^, 
allows  an  exchange  of  energy  between  the  spin  system  and  the  alter- 
nating field,  and  this  exchange  can  be  obaerved  electronically  if  the 
net  exchange  is  net  zero. 

The  nucleus  experiences  not  only  the  external  field  H  ,  but  also 
fields,  both  electric  and  magnetic,  due  to  its  surroundings.   Neigh- 
boring nuclei  with  magnetic  moments  influence  the  nucleus  of  interest 
through  dipole-dipole  coupling.   Net  electrostatic  forces  may  be 
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exerted  on  a  nucleus  with  a  nonspherical  charge  distribution.   While 
a  molecule  is  rotating,  its  charges  cause  a  magnetic  field  at  the 
sites  of  its  own  nuclei,  giving  rise  to  spin-rotation  interaction. 
There  are  other  interactions  which  are  less  important  in  the  present 
work,  and  will  be  disregarded. 

An  isolated  nucleus  would  resonate  at  cu  ,  providing  no  informa- 
tion other  than  the  Zeeman  frequency.   In  practice,  effects  such  as 
those  mentioned  above  modify  the  behavior  of  the  nucleus  by  adding  to 
the  applied  field  or  by  causing  other  forces  on  the  nucleus.   The 
nuclei  therefore  resonate  at  different  frequencies,  depending  on  their 
surroundings.   The  frequency  spectrum  therefore  provides  information 
about  the  surroundings  and  about  the  interactions  taking  place. 

In  order  for  there  to  be  a  net  exchange  of  energy  between  the 
alternating  field  and  the  spin  system,  the  spins  must  have  a  net  ex- 
change of  energy  with  a  third  system,  the  lattice.   This  takes  place 
through  the  interactions  mentioned  above  as  well  as  through  some  other 
interactions  less  important  here,  and  the  rate  of  the  interaction  is 
characterized  by  a  spin-lattice  relaxation  time  T  . 

Much  of  the  discussion  below  will  follow  Abragam  [4],  Andrew  [5], 
Cohen  and  Reif  [b],  or  possibly  other  standard  references,  and  will  not 
be  specifically  referenced  at  each  step. 

2.2   Intramolecular  Dipolar  Coupling 

Consider  one  nucleus   i  at  the  origin  P,  and  a  second  nucleus 

j   at  r,  with  magnetic  moments  p.  and  \i . ,    respectively.   The  magnetic 

■■»- 
potential  at  P  due  to  u,    is 
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<1).  .  =  -u.-V(l/r  ) 


(2.2.1) 


and  the  magnetic  field  is 


(2.2.2) 


Performing  the  indicated  operations  results  in 


(2.2.3) 


The  dipolar  Hamiltonian,  obtained  by  letting  H^  ^  H_  in  Eq.  (2.1.2), 


IS 


\  ■  -h-\i  -   f^-j  -  «"i-ij><":-ii'/<ll^'J  •  "•'•'' 


as  H   interacts  with  the  nucleus   i   the  same  way  as  H  .   If  an  ex- 
ij  ° 

ternal  field  H   =  zH   is  now  applied,  where  H   >>  H   ,  the  total 
o     o  o     J-j 

Hamiltonian  of  the  two  nuclei  is 


?/  =  - 


Y.-hH  I. 
1   o  iz 


^2 


3 
r.  , 
13 


h-^i 


3(I.-r..)(I.T  ) 

JL—iO. J ?J_ 

2 


(2.2.5) 


Letting 


W. 


1.1 


I.-I. 


3(I.-r..)(I.T..) 


2 
r. . 

2-3 


(2.2.6) 
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it  is  customary  to  write 


W..  =A+B+C+D+E+F  {2.1.1) 


where 


y 


A  =  I.  I.  (1  -  3cos^e) 

12  J  z 


B  =  -('<:)  (I. ^I._  +  I._I.^)(1  -  3cos^9) 


C  =  -(3/2)(I.  I.,  +1.  I. ,)sinecos9e  ^* 

IZ  J+      jZ  1+ 


D  =  -(3/2) (T.  I.   +  I.  I.  )sinecosee^*  =  C* 

12  J-     JZ  1- 


E  =  -(3/4)(I^^I.^)sin^ee  '^^ 


F  =  -(3/4)(I._I._)sia-ee^^*  =  E*   .  (2.2.8) 


The  field  (2.2.3")  H,.  is  on  the  order  of  1  G.  whereas  H  is  10 

or  10  G,  so  //  may  be  treated  as  a  perturbacicn  of  ?/^,   Restricting 

the  problem  to  like  spin  =5  nuclei,  the  zero-order  energy  level  is 

E  =  -yhH  (m.  +  m.),  and  according  to  first-order  perturbation  theory, 
Z       o  1    J 

only  those  parts  of  the  perturbing  Hamiitonian  which  induce  no  net 
chanee  in  M  =  m.  +  m.  contribute  to  the  energy.   Inspection  of  Eq. 
(2.2.8)  shows  only  terms  A  and  B  fit  this  condition,  and  terms  C 
through  F  may  be  disregarded  here.   These  tci-nis  actually  permit  weak 


-10- 


resonances  at  0  and  2cj   frequencies,  which  are  of  no  interest  in  the 
present  work.   The  truncated  Hamiltonian,  good  to  an  excellent  approxi- 
mation, is  therefore 

2  2 

}/=  -Y^H    (I.   +  I.  )  +^4^  (l-3cos^e)[l.  I.   -  (%)(!. ^I.   +  I.  I.J] 
o   iz    jz     3  iz  jz        x+  J-    1-  j+ 

(2.2.9) 

It  is  noted  that  term  A  can  be  interpreted  as  nucletis  i  experi- 
encing the  average  field  due  to  nucleus  j,  whereas  term  B  is  the 
simultaneous  'flipping'  of  both  spins,  which  can  occur  v/hen  i  and  j 
are  like  nuclei  and  therefore  in  resonance. 

There  are  four  states  for  the  case  of  two  spin  ^  nuclei,  repre- 
sented by  !++>,  |+->.  |-+>>  I — >•   Term  B  has  no  diagonal  matrix 
elements  in  this  representation;  however,  liivta.:  combinations  may'De 
used  to  define  a  triplet  state. 


l+>  =  !-H-> 


|0>  =  (l//?)(|-l->  +  i-+>) 

|->  =  I— >  (2.2.10) 


and  a  singlet  state,  1 0>  =  (]./v2)  (\+->   -    |— r>).   The  singlet  state 
has  no  Zeeman  energy  and  isn't  coupled  to  the  triplet  state  by  either 
JL    or  ^  ,  and  may  be  ignored. 

Energy  levels  are  readily  calculated  and  are 


E,  =  -YfiH  +  (Y^h^/4r^)(l  -  3cos^e) 
1       o 
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Eq  =       -  (Y^ti^/2r^)(l  -  3cos^9) 


E   =  YfiH   +  (Y^1i^/4r^)(l  -  3cos^9)    .       (2.2.11) 
-1      o 


In  the  remainder  of  this  section  only  identical  spin  %  nuclei  in 
diatomic  molecules  will  be  considered,  whose  intramolecular  separation 
is  significantly  less  than  the  intermolecular  nuclear  separation. 
Thus,  the  i,j  subscripts  may  be  dropped,  as  the  l/r"*  factor  insures 
intramolecular  effects  dominate. 

Two  resonant  frequencies  may  be  observed  at  constant  H  , 


fiw'    =  E   ,    -   E^  =   Yhln     +  a(l   -   3cos^9)] 
-10  o 


y 


9 

-Rw"   =   E„   -   E,    =   YRfH      -   a(l   -    Scos^e) ]  (2.2.12) 

0  1  o 


•3 

where   a  =   3Yfi/4r    .      Or  H      could  be  varied   and   w  =  cj'    =  cj"   held   fixed, 

o 

the  resonant  condition  being  found  at  two  values  of  H  , 


H   =  H*  ±  a(3u^  -  1)    ,  (2.2.13) 

o 


where  H*  =  co/v  and  u  ■■=   cosO.   H*  is  the  constant  central,  or  Zeeman, 
field,  and  u  is  the  cosine  of  the  angle  9  betv/een  H  and  r. 

In  a  polycrystalline  sample,  any  orientation  is  equally  likely, 
and  for  each  orientation,  there  will  be  two  resonant  lines,  above  or 
below  H*  by  an  amount 
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h  =  H  -  H*  =  ±(3u^  -  1)    .  (2.2.14) 

o 


The  spectral  distribution  is  just  a  plot  of  the  density  of  lines  (or 

density  of  orientations)  versus  h,  so  the  density  of  orientations  must 

be  known  as  a  function  of  u. 

A  sphere  can  be  defined  by  all  possible  positions  of  the  second 

nucleus  if  the  first  is  held  fixed,  and  the  density  of  positions  is 

constant  over  the  surface  of  the  sphere.   Remembering  rcos9  =  ru  is 

the  projection  of  r  along  H  =  zH  ,  it  is  noted  that  any  value  of  u 

defines  a  circle  on  the  sphere,  and  a  plane  containing  the  circle  is 

normal  to  z  at  a  distance  ru  above  the  center  of  the  sphere.   Another 

plane  and  circle  are  defined  by  u  +  Au,  the  planes  being  separated  by 

rAu.   The  surface  of  a  sphere  between  two  planes  separated  by  d  is 

y 

2Tfrd,  so  the  surface  area  between  the  planes  defined  by  u  and  u  +  Au 

2 
is  2rr  Au,  independent  of  u  itself.   Therefore  any  value  of  u  is 

equally  probable,  and  the  spectrum  can  be  considered  a  plot  of  the 

density  of  u's  versus  h,  or  du/dh. 

From  (2.2.14) 


u  "^  (^  +  l)^  (2.2.15) 


and  the  spectral  distribution,  f(h),  is 


'"-"^'-^-^^.^'y'"        '  (2.2.16) 


In  the  case  of   N,"  y  ='  ~|Y|.  so  (2.2.15)  becomes 
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-?Jal 


Z\a\     h-> 


Figure  1.   Pake  powder  pattern  doublet.   The  theoretical  unbroadened 
spectrum  (dashed  line)  is  the  sum  of  the  two  components 
(dotted  lines)  resulting  from  the  transitions  shown  which  are 
for  the  case  of  a  negative  gyromagnetic  ratio.   A  typical 
spectrum  including  intermolecular  broadening  is  shown  by  the 
solid  line. 
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f(h)  cc  [-_-^.   +  i)-^  (2.2.17) 

the  upper  (lower)  sign  corresponding  to  the  [->  ->  [0>(|0>  ->-  |+>) 
transition.   Fron  (2.2.14),  we  note 

-2|a|  1  h  £  |a|        |->  ->  |0> 

-hllh£2|a|  l0>  ->-    |+>        .  (2.2.18) 

The  two  curves  overlap  in  the  region  -|a|  £  h  _<  |a|,  and  result 
in  a  spectrum  indicated  by  the  dashed  line  in  Fig.  1.   Intermolecular 
dipolar  coupling  causes  each  line  to  be  broadened,  resulting  in  a 
curve  such  as  that  shown  in  Fig.  1  by  the  solid  line.   The  spectrum 
sho^^i  in  Fig.  1  is  known  as  the  Pake  doublet  after  G.E.  Pake  [7].   It 
should  be  noted  that  the  curve  will  be  reversed  if  frequency  units 
vice  magnetic  field  units  are  used,  in  addition  to  a  scale  factor  set 
by  the  gyromagnetic  ratio. 

2.3  Nuclear  Quadrupole  Interaction 

The  charge  distribution  in  a  nonspherical  nucleus  exhibits  an 
orientational- dependent  interaction  with  the  surrounding  electric 
field  due  to  the  external  charge  distribution.   Classically,  this  may 
be  written  as 

E  =  /  p(r)(t)(r)d\  (2.3.1) 

where  p(r)  is  the  nucleai  charge  density  and  :J)(r)  is  the  potential  due 
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to  external  charges. 

One  convenient  method  of  treating  the  probleir  is  a  Taylor  series 
expansion  of  <p(v)    about  the  origin: 


E  =  9(0)/pdr+  I   <i>    J  X   pd   r  +  ^-r     I      (f-^/x  x  od  r  + 

^  a  •'  a       2!   '-^  a3    a  B 
a  a, 6 


3   .  1   r   ,    f      ,3 

(2.3.2) 


where  x  (a  =  1,2,3)  stands  for  x,  y,  z  respectively,  and  where 


a   dx 
a 


3^(j) 


'  ^aS,        dx   3x„ 
r=0  "  ^ 


etc.       (2.3.3) 


The  first  term  in  the  expansion  is  the  electrostatic  energy  of 
the  nucleus  taken  as  a  point  charge,  and  is  of  no  interest  as  it  is 

y 

just  a  constant  with  no  orientational  dependence.   The  second,  or 
dipole,  term  vanishes  as  nuclear  theory  tells  us  the  wave  function  has 
definite  parity  so  p(r)  =  p(-r),  and  since  x  changes  sign  in  opposite 
quadrants,  the  dipole  integral  vanishes.   For  the  same  reason,  all  the 
even  terms  vanish. 

The  terms  in  the  expansion  are  referred  to  as  monopole,  dipole, 
quadrupole,  octopole.  hexadecapole,  etc.,  as  they  are  identified  with 
the  type  of  charge  distribution  which  would  cause  the  corresponding 
potential.   We  hav'e  shown  the  dipole,  octopole,  and  every  other  higher 
term  vanishes  and  the  monopole  term  exists  but  is  of  no  interest. 
This  leaves  the  quadrupole,  hexadecapole,  and  higher  terms: 


1   V       r  3 

E=T  /   <}'   ^Jx  x^cd  r  +  hexadecapole  and  higher   .     (2.3.4) 

a, 3 
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Now  *(0)  is  roughly  e/r  ,  where  r  is  of  the  order  of  the  radius 
■  ^  '       °  •'    e         e 

of  an  electron  orbit,  so  the  monopole  term  >f>(0)  J  ad  r  is  on  the 

9  ^  o 

order  of  Ze  /r  .   The  next  nonzero  terra,  3  (4/3x  9x„  J  x  x^pd  r,  is 
e  a   B  "^   a  p 

3   2        ?   2    2 
on  the  order  of  (e/r  )(r  Ze)  =  (r'/r  )(ZeVr  )  where  r   is  on  the 

en        n   e       e         n 

2   2     -8 
order  of  a  nuclear  radius.   As  r  /r   ^10   ,  the  quadrupole  term  is 

n   e 

—  8 
on  the  order  of  10   times  the  monopole  term.   Likewise,  the  hexa- 

_  Q 

decapole  term  is  on  the  order  of  10    times  the  quadrupole  term  and 

so  oa  for  succeeding  terms.   Thus,  the  quadrupole  term  is  the  only 

one  of  interest  within  the  resolution  of  m.ost  spectrometers. 

If  we  replace  p  with  its  quantum  mechanical  equivalent  p    and 

op 
o 

define  the  symmetric  tensor  Q'   -  f  x  x„p   d  r,  we  have  the  quadrupolar 

aQ        •"      a   3  op 

Ikimiltonian   to   an  excellent   approximation, 


01,3 


It  is  convenient  to  substitute  the  symmetric,  tracelesG  tensor 


Q  o  =  3Q'   -  <5    y  Q'  (2.3.6) 

a3     a3    a(3  ^     aa 
a 


in  Eq.  (2.3.5)  which  then  becomes 


'H-  =  Z   I      ^    .0.   a+M)   ^     )[)   Q'  ]     •      (2.3.7) 

a ,  3  a      a 


The  potential  is  due  to  charges  outside  the  nucleus,  and  therefore 

satisfies  Laplace's  equation,  to  i>        is  traceless  and  hence  the  second 

aa 

terra  in  Eq.  (2.3.7)  vanishes.   Here  we  are  neglecting  any  electronic 
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charge  distribution  which  is  overlapping  the  nucleus  as  it  is  very 
small,  especially  after  subtracting  the  spherically  symmetric  part 
which  is  of  no  interest  in  orientational  studies.   We  are  therefore 
left  with 


%  =  i  I  VaB       •  ^2.3.8) 


a,B 


From  Eq.  (2.3.6)  and  the  definition  of  Q\,   we  have 

aB' 


2x    ,3 


•^uB  =  -^  ^^Vb  -  'aB*^  ^  V^    •  ^'-'-'^ 


Putting  in  P^„(r)  =  V  q  5(r  -  r  ),  and  performing  the  integration 
op      k 

results  in 


QaS  =  ^^  ^^^k^Bk"  V>  ^2.3.10) 


where  e  is  the  charge  of  a  proton  and  the  sum  is  over  the  protons  in 
the  nucleus.   The  expression  in  tne  parentheses  could  be  written 


„    _  3(^k^3k-*-  V.'^ak^    ,   .2  ,,  .^  ,^, 


We  now  look  at  another  operator 
and  assert 
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<ImlQ  „|lm'>  =  C<Iml(3/2)(I  I^  +  I^I  )  -  6  „I^Ilm'>   .   (2.3.13) 


We  notice  the  two  operators  are  constructed  in  the  same  way  and 

that  Q  „  is  a  linear  combination  of  the  Q  ^  's.   The  two  operators 

aB  aSk 

obey  three  conditions:   1)  they  are  synnnetric,  2)  they  are  traceless, 
and  3)  they  transform  under  rotation  of  coordinate  axes  in  the  same 
way  as  second  order  spherical  harmonics.   This  last  condition  is  equiv- 
alent to  saying  they  are  second  rank  tensor  operators,  or  that  they 
satisfy  the  same  commutation  rules  with  respect  to  I.   These  three 
conditions  are  sufficient,  using  group  theory,  to  prove  Eq.  (2.3.13) 
[S].   A  more  sophisticated  proof  calls  on  the  Wigner-Eckart  theorem 
which  shifts  the  complexity  to  the  proof  of  the  theorem.   A  simpler, 
more  direct  proof,  but  also  more  tedious,  comes  from  direct  matrix, 
multiplication  [8], 

Remembering  Q   was  related  to  the  quadrupole  portion  of  the 
nuclear  charge  distribution,  we  define  the  quadrupole  moment  Q  by 


eQ  =  <II  Q   II>  (2.3.14) 

7.Z. 


which  is  used  in  Eq.  (2.3.13)  to  obtain 


eQ  =  C<Il|31^  -  I^IlI>  -  CI(2I  -  1)  (2.3.15) 


from  which  we  obtain  C  -  eQ/I(2T.  -  1). 
Equation  (2.3.8)  can  be  rewritten 
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u  = 


-  i    eQ 


6  1(21  -  1)   ^-  *^aB 
a,  ii 


■f  (I  I.  +  I-I  ) 

2   a  6    S  a 


(2.3.16) 


We  can  diagonalize  (j)  ^  and  hereafter  assume  we  have  done  so  and 


use  a,B  without  primes  to  indicate  the  coordinate  system  in  which  <^ 


aS 


is  diagonal.   Thus, 


\   ^  ^  Ti^^-1)    ^  "^aa^'^l  -   ^'^         '  ^'^•^•^^> 


Labeling  the  axes  is  arbitrary,  so  we  define  a  principal  axes 
frame  by 


'  zz  '  —  '  yy '  —  '  XX ' 


(2.3.18) 


and  in  this  frame 


'Wr.     = 


1    eQ 


'Q    6  1(21  -  1) 


?    >-9  2    -v2  2    ->-2 

t.    (31    -  r)  +  cp    (3r  -  I")  +  t    (31    -  I  ) 

XX   X         yy   y         'zz   z 


(2.3.19) 


With  the  help  of  Laplace's  equation,  ^        +  <j)   -y-   ^        =  0, 
*^      '^        ^         XX    yy    zz 


■H  =  -  -r-^ r  r*  (31^  -  1^^  -^  (*   -  *  )(i2  -  1^)1 

"q   4  1(21  -  1)  L^zz.^   z     ■    ^xx   ^yy''  x    y  J 


(2.3.20) 


Introducing  A  =  e  qQ/4I(2I  -  1),  n  -  (^   -  (j)  )/^  _,  eq  =  4)  ^,  and 

the  raising  and  lowering  operators  I^  =  (I  ±  il  ),  we  obtain 

—     X     y 


T  ? 


?/q  =  A  31^  -  r  ^  (n/2)(i;  +  !_) 


(2.3.27.) 
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If  we  were  to  calculate  matrix  elements  for  I  =  0  or  I  =  ^  using 
Eq.  (2.3.17),  we  would  find  all  terms  vanish,  which  means  the  electric 
quadrupole  interaction  exists  only  for  I  >_  1.   Semiclassicaily,  we 
would  expect  no  orieiitational  dependence  for  1=0  which  has  no  z  com- 
ponent, and  I  =  ^2  has  two  states,  +^2  and  -%,  which  correspond  to  a 
reversal  of  spin  direction  but  no  change  in  charge  distribution  and 
again  no  orientational  dependence. 

We  have  introduced  the  electric  quadrupole  moment  Q,  the  electric 

field  gradient  (in  the  direction  of  maximum  gradient)  q,  the  asymmetry 

2 
parameter  n,  and  a  quantity  A  =  e  qQ/4I(2I  -  1).   We  notice  n  =  0  when 

(})   =  4)   ,  which  occurs  in  spherical,  cylindrical,  or  cubic  symmetiry, 
XX    yy 

hence  the  name  asymmetry  parameter.   The  maximum  value  n  can  have  is 

unity  when  *   =0  and  I A   I  =  I i   |.   We  also  note  that  a  and  Q  can 

^       XX  yy    '  zz' 

not  be  determined  individually  in  a  nuclear  quadrupole  resonance  ex- 

2 
periment  as  they  appear  only  as  a  product.   The  quantity  e  qQ'h  is 

called  the  quadrupole  coupling  constant,  and  it  is  convenient  to 

include  it  in  the  constant  A. 

The  quantization  direction  of  the  nuclear  quadnjpcle  interaction 

is  set  by  the  z  principal  axis  which  is  the  direction  of  the  maximum 

electric  field  gradient.   If  other  interactions  are  present,  they  in 

general  will  have  different  quantization  directions.   In  the  presence 

of  an  external  magnetic  field,  for  instance,  the  Zeejcan  Hamiltonian 

Jj     (discussed  in  Sec.  2.2)  is  diagonal  in  a  coordinate  system  with 

the  z  axis  parallel  to  the  external  field  which  in  general  does  not 

coincide  with  the  z  principal  axis.   Low  magnetic  field  and  high 

magnetic  field  cases  can  be  treated  by  first-order  perturbation  theory. 

The  intermediate  case  with  the  tvo  interactions  of  comparable 
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magnitude  must  be  solved  exactly  and  results  in  a  much  more  complex 
orientational  dependence.   The  low  magnetic  field  case  is  of  little 
interest  in  the  present  experiments  and  will  not  be  discussed.   The 
high  field  case  will  be  treated  in  Sec.  2.4  and  the  intermediate  re- 
gion in  Sec.  2.5. 

2.4   High  Field  Nuclear  Ouadrupole  Interaction 

If  a  system  of  identical  spins  is  placed  in  a  magnetic  field, 
the  well-known  Zeeman  Hamiltonian  is 


\   =  -^"o^z- 


(2.4.1) 


from  Eq.  (2.1.2),  the  prime  being  used  to  distinguish  the  coordinate 
system  in  which  7/     is  diagonal. 

If  nuclear  quadrupole  interactions  are  also  present,  we  must  add 
the  quadrupolar  Hamiltonian,  Eq.  (2.3.21): 


'V     =  y^  +  X   =   -Y^H  I  ,  +  A 
Z    Q       o  z 


2   ->2 
31  -  I 
z 


f  (^+  +  ^b 


(2.4.2) 


The  unprimed  coordinate  system  refers  to  the  principal  axes  frame. 

Restricting  the  discussion  to  the  high  field  approximation  where 
first-order  perturbation  theory  is  valid  for  "}!  ,   we  are  only  inter- 
ested in  the  diagonal  matrix  elements,  <m|?/|m>.   It  is  convenient  at 

this  point  to  return  to  Eq.  (2.3.15)  which,  using  ;{)   +  t.  .  "^  '^   =  0, 

ji.j'L  yy  Lt£» 

becomes 
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V^   =  ^77t|? 7Y  (6  I^  +  *  I^  +  *  I^)    .     (2.4.3) 

Q   21(21  -  1)   XX  X   yy  y    zz  z 


Writing  the  direction  cosine  of  x  with  resoect  to  x'  as  a  „, 

a         *        6     a3' 

whare  x  ,  a  =  1,2,3  stands  for  x,  y,  z,  and  x'  B =  1,2,3  stands  for 
a  t$ 

X '  ,  y '  ,  z  '  ,  we  have 


^n    ^    ^      ^«R'^R'    •  (2.4.4) 

a     ,      an      a 


Thus , 


2    2   2     2   2     2   2 
I  =  a   T  ,  +  a^„I^,  +  a,„I  ,  +  cross  terms  (2.4.5) 

X    11  X     12  y     13  z 


2      2 

and  likewise  for  I  and  I  .   We  are  considering  only  diagonal  matrix 

elements  and  cross  product  terms  like  <mla,,a,^(I  , I  ,  -r  I  , I  ,)!m> 

11  12  X  y     y  X 

vanish  and  may  be  dropped  for  our  purposes. 

Putting  Eq.  (2.4.4)  in  Eq.  (2.4.3)  results  in 


=    eQ 
Xj  "  21(21  -  1) 


^xx^^l^x'  +  42^J-  -^  43^z') 


-   *yy^^21^x'  ^  ^22^;-  ^  ^23^;-^  "  ^zz^^n^x'  ^  ^32^y'  ^  ^33^z'>. 

(2.4.6) 

It   is   easj'   to   show  that 

<m|  I_  ,  lra>   =    <m|  I",  |m>   =  i2<m|  I      -    I    ,  [m>   =  Js[l(I  +  1)    -  m"]      . 

(2.4.7) 
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As  the  only  allowed  transitions  are  those  for  which  Am  =  il,  we 
calculate 


<mll^,  |m>  -  <m  -  l|l^,  [m  -  1>  =  -!i(2ra  -  1) 


<m|l'',  |m>  -  <Tn  -  ijl^,  Im  -  1>  =  (2m  -  1)   .        (2,4.8) 


The  transitions  AE^,        ,,  =  <n\y^\m>   -  <m  -  l|*V^!™  "  1>  a'^e  then 

Q(m  ^-  m  -  1)      '  Q '  ' '  Q' 

evaluated,  yielding 


=  ^Q(2n-  -  I). 
Q(m  -^  m  -  1)    21(21  -  1) 


2     2  2    2 

■   r  ^11   ^12  ^  2  .  _^,  (    ^21  ^22^  2  >. 


2     2  '' 

a. 


.  r      ^-31   ^32     2  ,] 
^zz  ^-  "1 2"  "  ^33^  J 


eQ(2m  -  1)  ,  3,    2  ^^  3,    2    3,    2  .       /i  /  on 
21(21  -  1)   2  XX  13   2  yy  23   2  zz  33 


From  Eq.  (2.4.8)  v;e  see  that  E^,    ,  t  ^   \  =  ~E  ,        ^^ 
^  Q(-m  +  1  -^  -m)     Q(m  ^  m  -  1) 


and  from 


Sec.  2.1,  AE^,        ,,  =  AE,.   _^  .      .  =  -yhH  .   Hence  the  ob- 
Z(m  -^  m  -  1)      Z(-m  +  1  ->-  -m)       o 

served  transitions  are  in  pairs. 


AE+  =  AE  ±  AE     ,  (2.4.10) 

—     ^     4 


where  the  upper  sign  signifies  the  (m  -*■  m  -  1)  transition  and  the  lower 
sign  the  (-m  +  1  ^  -ra)  transition. 

We  are  not  ccncerned  with  the  Zeemau  term  which  depends  only  on 
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Projection 
on  z  axis 
+  or  -. 


Projection 
on  z  axis 
always  +. 


a.  Vector 


Tensor 


c 

o 


Projections  of  vector  (  

component  (  )  on  z  axis. 


)  and  tensor 


Figure  2.   Comparison  of  zr-axis  projections  of  rotating  vectors  and 

tensor  ccrrponents.   The  period  of  che  tensor  component  is  h 
that  of  the  vector. 
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H  and  the  nuclear  species.   It  is  sufficient  to  measure  the  difference 
o 

in  energies  between  the  two  transitions,  which  is  2AE  .   In  frequency 
units. 


3eQ(2m  -  1)  , ,    2   ,  ^    2   ,  ,    2  .        /o  /  -nx 
2^^  =  yi(2f^^T)h  ^*xx^3  -^  *yy^23  "^  *zz"33^   '     ^^-''-^^^ 


This  is  a  convenient  point  at  which  to  note  the  physical  inter- 
pretation of  the  frequency  separation  2Av   in  terms  of  the  orientation 

of  the  external  field  H  with  respect  to  the  principal  axes  x,  y,  z. 

o 

We  recall  a^ _  is  the  direction  cosine  of  the  x  principal  axis  with 

respect  to  the  direction  of  H   (which  is  parallel  to  z')  and  likewise 
•^  o 

2 
for  a„_  and  a__.   Thus,  <J)  ,a    is  the  magnitude  of  the  x  component  of 

the  electric  field  gradient  tensor  times  the  square  of  the  cosine  of 

the  angle  between  the  x  component  and  H  ,  and  so  on  for  the  other 

terms.   We  note  d)   and  ^        have  the  same  sign,  which  is  opposite  to 
XX      yy 

the  sign  of  ^ 

zz 

The  observed  frequency  shift  Av  is  proportional  to  the  sum  of 
the  three  'projections',  remembering  the  z  component  is  opposite  in 
sign  to  the  others. 

We  also  note  hei a  the  significance  of  the  square  of  the  cosine. 
If  the  EFG  components  were  timple  vectors  which  have  a  direction  and 
sense  (or  sign),  wb  would  s-ipect  any  effects  to  reverse  signs  every 
180°  as  shown  in  Fig.  2.   Noting  in  Fig.  2(c)  that  the  tensor  'pro- 
jection' goes  through  two  complete  cycles  in  a  360"^  rotation,  we 
expect  a  26  to  appeal  which  will  come  from  the  cosine  squared  factors. 

In  general,  the  orientations  of  the  principal  axes  are  unknown 
as  are  the  magnitudes  of  the  4   's.   We  therefore  express  the  EFG 
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tensor  in  the  lab  frame.   Using  a^^,    as  before  and  the  fact  that  <J)^g 


is  diagonal  (<p^^   =   <5^Q*„fl)  we  have 


o 


and  therefore 


o 


Finally,  Eq.  (2.4.11)  becomes 


^ 


o»    3eO(2m  -  ])  ,     _  y.  r?  A  lA') 

2AV  =  21(21  -T)h  '^z'z'  -  ^*z'z'  (-4.14) 


with  the  definition  of  K  obvious. 

2  2 

We  are  interested  in  finding  e  qQ/h  (or  A  =   e   qQ/4I(2I  -  l))j  n, 

and  the  orientation  of  the  principal  axes.   By  measuring  2Av  at 

enough  different  orientations  and  usiMg  the  transformation  Eq.  (2.4.13), 

the  desired  quantities  csn  be  worked  out.   Equation  (2.4.13)  is  not 

very  convenient,  however,  -ie  the  a   ,  a   ,  a^   terms  are  different 

for  each  orientation.   We  ther^.fcie  pick  a  crystal  frame  X,  Y,  Z  which 

is  known  and  may  contain  some  of  the  crystal  axes.   The  transfonnatioa 

from  the  lab  x' ,  y' ,  z'  ^ramd  to  the  X,  Y,  Z  frame  is  known,  and  the 

transformation  from  the  X,  Y,  Z  frame  to  the  x,  y,  =  principal  axis 

frame  will  be  constsr.t  af  both  frames  are  fixed  in  the  crystal. 
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The  procedure  is  to  find  Che  EFG  tensor  in  the  crystal  frame  by 
experimental  methods,  then  diagonalize  the  tensor.   Diagonalizatioa 
will  yield  the  EFG  tensor  in  the  principal  axis  frame  as  well  as  its 
orientation  with  respect  to  the  crystal  axes.   Volkoff,  Fetch  and 
Smellie  [9]  outlined  the  procedure  in  1952  and  it  has  become  a  rather 
standard  method  since  then.   The  Volkoff  method  more  recently  was 
shown  to  have  possible  ambiguities  [lO].   Furthermore  it  requires 
more  data  than  is  necessary  in  many  cases  as  discovered  in  the  present 
work  (discussed  in  Chap.  4)  and  as  recently  reported  by  El  Saffar 
[ll].   Nevertheless,  the  Volkoff  irethod  provides  a  very  convenient 
starting  point  and  it  is  easy  to  modify  it  to  reduce  the  quantity  of 
data  required. 


The  tensor 


y 


A      4 


'xz 


'XY 


*YY   *YZ 


*XZ    * 


YZ   ^ZZ 


(2.4.15) 


is  symmetric  as  displaj^ed  and  has  only  5  independent  elements  using 
the  vanishing  trace  property.   Diagonalization  gives 


XX 

0 

0 

0 

yy 

0 

0 

0 

* 

zz 


-eq(l-n)/: 


^ 


-eq(l+n)/2   0 


eq 


(2.4.16) 


and  the  direction  cosines  oL   the  principal  axes.   We  recall  the 
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y',x 


X 
Y 


z'cosOn^  +  x'sinSj^ 


-z'  3iri6 


X 


x'cos9. 


Figure  3.   Transforrr.atiou  between  the  lab  axes  and  crystal  axes  during  ar 
X  rotation.   Cyclic  permutations  of  X,  Y,  Z  result  in  the 
transformations  for  Z  and  Y  rotations. 
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elements  all  contain  the  factor  K  =  3eQ(2m  -  1)/2I(2I  -  l)h  so  only 

2 
e  qQ  is  actually  determined,  rather  than  eq  and  eQ. 

We  start  experimentally  by  selecting  crystal  axes  X,  Y,  Z  and 
aligning  the  crystal  in  a  magnetic  field  H  =  z  H  with  X,  Y,  Z  paral- 
lel to  y',  z',  x'  initially.   The  crystal  is  then  rotated  about  the 
X  =  y'  axis  through  an  angle  6   and  the  frequency  separation  2Av 

A  X 

observed  and  plotted  versus  9  ,  the  plot  being  termed  the  X  rotation 

A 

pattern.   Figure  3  shows  the  relation  between  the  crystal  and  lab 
frames  for  an  X  rotation. 

Using  the  transformation  equations  in  Fig.  3,  we  obtain 


=  [cose^  —  -  sme^  —J  <}. 


d  ,d  ,  ^      '-   X  9Y     X  az-* 

z   z 


=  '2(*YY  -^  *ZZ>  -^  '^(*YY  -  ^ZZ^-^^^^e^  -  \,__sin2e^  (2.4.17) 


which  is  put  in  Eq.  (2.4.14),  yielding 


2Av^  =  A^  +   B^cos29j,  +  C^sin29^ 


=  A^  +  (bJ  +  cJ)\os2(e,.  -  6.^)  (2.4.18) 


where 


Ax  =  ^K(*YY+  V  ^   -''^^*XX 


Bx  =  ^K(*yY-  W 
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S  ==  -"^^YZ 


tan26^  =  -^     .  (2.4.19) 

X   B^ 


Thus  the  X  rotation  yields  the  diagonal  elements  and  <{>„.,•   Similar 
equations  may  be  obtained  by  cyclic  permutation,  so  (J>   and  ^        may 

A^         Al 

be  obtained  from  Y  and  Z  rotationfi,  respectively. 

Diagonalization  of  the  EFG  tensor  requires  solving  the  cubic 
equation 


Y^  -  ay  -  b  =  0   ,  (2.4.20) 


where  y 


11  1  2 

^^XY    YZ    ^XZ    ^XX^YY    ^YY^ZZ    ^ZZ^XX"^ 


^  -  ^  ^*XX*YY'*ZZ  ■"  ^*XY*YZ*XZ  "  *XY*YZ  "  *Yy4z  "  *ZZ*L^  '  ^"•'^•^^^ 


We  find 


Y  =  2(a/3)^cos(a  -  2Tin/3),  n=  1,2,3  (2.4.22) 


cos3a  =  (|b|/2)(3/a)^^^     0  <  a  <  tt/6 


Y,  =  K.t   ,      Y_  =  K(}.   ,      Y_  =  K(J) 
1     XX      '2    ^yy      '3     zz 
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These  choices  insure 


Y3  =  IY3I.    IY3I  >  IY2!  >  IYj^I   .  (2.4.23) 


and  finally  we  obtain 


n  =  (Yj^  -  Y2)/Y3  =  3\aaa  (2.4.24) 


and 


le^qQ/h|  =  |eQ(^^^/h|  =  [21(21  -  l)/3(2m  -  1)]-1K({>^^I   .     (2.4.25) 


The  diagonalization  process  yields  the  direction  cosines  of  the 
X,  y,  7.   principal  axes  with  respect  to  the  X,  Y,  Z  crystal  axes: 


A     p     V  . , 

"     "  -^  (2.4.26) 


^1     D-,     D-    ,^2  ,  _2     2  ^i. 
In    2n    3n   (D,  +  Do   +  0_  ) 

xn    2n    3n 


where  X  ,  y  ,  v  are  the  direction  cosines  of  x  with  respect  to  X,  Y,  Z 

and  likewise  for  n  =  2  and  3. 

The  relative  signs  of  A  ,  u  ,  v   for  any  n  are  fixed,  but  all 

n   n   n 

three  signs  may  be  reversed  without  reversing  the  signs  for  the  other 
values  of  n.   Thus,  a  right-  or  left-handed  coordinate  system  may  be 
chosen. 

Volkoff  notes  that  only  five  readings  are  required  in  the  most 
general  case,  such  as  S„  =  9^  =  ^   =  45°  and  any  two  of  9   =  8„  = 

A       I       i.  A       I 

9„  -  0°.   This  would  require  a  very  accurate  knowledge  of  the 

4C1 
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orientations  and  in  the  case  of  inequivalent  sites,  identification  of 
the  transitions  might  be  difficult. 

Additional  data  increase  accuracy  and  may  be  helpful  in  detect- 
ing orientation  errors  and  in  identifying  the  transitions.   Some  of 
these  effects  as  well  as  the  single  rotation  method  will  be  discussed 
in  Chap.  4. 


2.5   Intramolecular  Dipolar  Splitting  of  the  NMR  Spectrum  of   N  in 

15   14 

N-   N  Mixed  Molecules 

15   14 
We  next  consider  the  case  of   N-  N  mixed  molecules  with  nuclear 

spins  1  ~  h   and  1=1  for  the  two  isotopes,  respectively.   The  total 

Hamiltonian  is  taken  as 


V     =  ^(15)  +  \{U)   +   7/(14)  +  ^^(15-14)   .       (2.5.1) 


15  14 

The  region  where  the   N  Zeeman  and   N  quadrupolar  energies  are  of 

15   14 
the  same  order  will  be  considersd.   The   N-  N  intramolecular  di- 
polar energy  is  much  less  and  may  be  treated  by  first-order  pertur- 
bation theory.   Wi  note  at  this  poi.nt  that  the  asymmetry  parameter  n 
is  zero  for  N-. 

We  proceed  in  the  following  manner.   The  ^1(1^)    term  in  zero 
field  gives  rise  to  two  energy  levels  corresponding  to  m  =  0  and 
m  =  il  with  state  functions  |l>,  |0>,  |-1>.   Adding  an  external  mag- 
netic  field  H   re.iicves  the  ciegeneracy  and  the  state  functions  become 
linear  combinations  of  |l>,  |0>,  |-1>: 


IV  ="ln>^^-^^2n'°^-'^3nl-^^    '         (^.5.2) 
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Working  in  the  principal  axes  (PA)  frame,  we  calculate  the  a   's 

mn 

using 'H(l'i)  and  V  (1^) .      Next,  T/(15)  (with  state  functions  \+h>y 

\-\>    in  a  frame  where  H  =  zH  )  is  expressed  in  the  PA  frame,  and 

o     o 

the  state  functions  are  linear  combinations  of  I+*2>,  |-%>: 


U-j->  =  Sj^+K'2>  +  82+l-'-5>    -  (2.5.3) 


To  first  order,  the  total,  state  functions  are  taken  as  products 
of  the  exact  uncoupled  state  functions. 


5  +>  =  k  >k+>   .  (2.5.4) 


■^  (15-14)  is  treated  as  a  perturbation  causing  a  first-order  - 

14 
shift  which  is  dependent  on  the  state  of  the   N  nuclei. 

From  Eqs.  (2.3.21)  and  (2,1.2)  we  have 


"^Q.Z^1^>  =  "^Q^^^^  ■'   \(^^^ 


=  A(3I"  -  f^)  -  vTif-H  (2.5.5) 

z  o 


where  x,  y,  z  is  the  PA  fraiL-e. 

As  the  molecule  is  axi.illy  syimnetric  (and  therefore  n  =  0)  ,  we 
may  arbitrarily  select  the  orientf)tiou  of  x  and  y  and  therefore  choose 
x  in  the  plane  of  H   ar.d  ;:.   Thvs.  K  ^  zH  ccsS  +  xH  sine  where  9  is 

^  O  GO  O 

-♦■ 

the  angle  between  H  and  z,  and 

o 
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.2   i-l 


^i      ^(14)  =  A(3I'  -  I  )  -  Yh(I  H  COG3  +  I  H  sin9) 
''q,Z  z         '    z  o        X  o 


=  A(3I~  -  I  )  -  Yhl  H  cos8  - 
z  z  o 


^'tiH  sine 


Introducing 


2YhH       „ 
T    o   „    E 

e  qQ 


=  cos9,  V  = 


sin9 

ri 


(2.5.7) 


we  obtain 


^Q.Z(^^> 


2A 


(31^  -  I^)/2  -  Lul   -  Lv(T^  +  I  )//2  1  .   (2.5.8) 
z  z       +    — 


Matrix  elements  are  calculated  and  displayed  in  matrix  form 


=  2A 


•3    -    I>u 

-Lv 

0 

-Lv 

-1 

-Lv 

0 

-Lv 

h  +  Lu 

(2.5.9) 


The  secular  deterrainant 


I?/-  e|  i  =  0 


(2.5.10) 


results  in  the  cubic  equation 


T^  -  1(3/4  +  L")  V   1/4  -  vl/2)L^(3u^  -  1)  =  0    ,     (2.5.11) 
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where  we   have  used  T  =  E/2A  from  (2.5.7).   The  solutions  of  a  cubic 
equation 


T   -  a,T  -  a  =  0 
1     o 


(2.5.12) 


are 


T_^   =  2(aj^/3)'^cos(4./3) 


T   =  2(a  /3)  'cos((j)/3  +  120°) 


^1  "  2(3^^/3)  'cos((()/3  +  240°) 


(2.5.13) 


-3/2 


,^ 


where  cos<|)  =  (a^/2)(a-/3)   '  .   Parker  [12]  has  tabulated  T  's  for 
0     1  n 

various  angles,  external  fields,  and  spins , 

The  assignment  of  subscripts  to  the  reduced  energies  T  is  made 

n 

by  examining  the  behavior  of  T   in  the  limit  of  small  H  where  -1,  0,  1 

n  o 

are  good  quantum  numbers.   In  the  region  of  interest,  linear  combina- 
tions must  be  used  as  expressed  in  Eq.  (2.5.2)  which  can  be  written  as 
a  column  vector. 


iV 


In 


On 


'-In 


(2.5.14) 


In  our  notation,  /il'li   >  =  E    1 6   >  becomes 

'    n  n'    n 
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h   -  Lu  - 


-Lv 


Lv    0 
1    -Lv 

Lv  ?!  +  Lu 


' 

^a 

"On 

< 

=  Tn 


• 

^'in 

"On 

"-m 

. 

(2.5.15) 


2     2     2 

which  can  be  solved  for  a   's  using  a,   +  ct.  +  a  ,   =  !• 

mn  In    On    -In 

Introducing 


^ln=  (^^-Ta-Lu).   C^^  =  Lv,   C_^^ 


=  (%  -  Tn  +  Lu), 


n     in  On    In  -In    -In  On 


(2.5.16) 


the  solutions  are 


^ 


a-   =  C_  C  .  /C  ,  a^  =  C  ,  C^  /C  ,  a  -   =  C.  C,  /C   .  (2.5.17) 
In    On  -In  n   On    -In  In  n   -In    On  In  n 


The  complexity  of  these  rather  simple  appearing  equations  can 

2 
only  be  appreciated  by  expressing  them  in  temis  of  H  ,  e  qQ,  and  9; 

however,  the  substitutions  in  (2.5.7)  and  (2.5.15)  greatly  simplify 

the  calculations. 

14 
Up  to  this  point  only   N  has  been  con'3idered.   From  here  on  we 

must  add  a  14  or  15  label  to  avoid  confusion.   The  Zeeman  Hamiltonian 

"U  (15)  expressed  in  the  PA  frame  has  a  fore  similar  to  W  (14)  in  Eq. 

(2.5.6): 


ri„(15)  =  -Y(15)^H  co3ei  (15)  - 
Z  o     z 


Y(15)-hR  sine 


[l^(15)  +  I_(15))  . 

(2.5.13) 
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In  the  frame  with  H   parallel  to  the  z  axis,  the  exact  states 

o 

are  \+h>,     \-h>,    but  here  H   forms  an  angle  9  with  the  z  axis  and 
linear  combinations  of  the  exact  state  functions  must  be  used: 


U^>  =  e^+|%>  +  ^2+\-h> 


'2± 


(2.5.19) 


in  column  vector  form. 


Letting  B  -  (h)y(.15)hYl   cosS,  C 


(J2)Y(15)tlH  sine,  Eq.  (2.5.18) 


becomes 


^^(15)  =  -2BI^(15)  -  C[I^(15)  +  I_(15)] 


-B   -C 
-C    C 


(2.5.20) 


in  matrix  form,  and 


-S   -C 


-C 


1 

'                         ' 

^11- 

J 

'2± 

=  +Y(15)riH 


1± 


'2± 


(2.5.21) 


2     2 
This  equation  together  with  the  normalizatioxi  condition,  B.._j.  +  S>j+  ~   !■» 


can  be  used  to  obtain 
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1+ 


COSJ  , 


'2+ 


.  9 


.  9 


^2-  ^  ^°^ 


(2.5.22) 


or 


V 


cos 


.  Q 
siny 


^  >  = 


.  9 


cosy 


(2.5.23) 


We  can  now  form  the  product  functions  \^    .  >  Eq.  (2.5.4), 
We  next  consider  the  perturbation  term. 


l/o(15-14) 


y(15:'y(14)1i 
3 


2 


1(15) -1(14)  -3 


(t(15)>?iri(14)-^j 


v^ 


(2.5.24) 


which  comes  from  Eq.  (2.2.5).   In  the  PA  frame,  r/r  =  z,  so 


7/_  (15-14)  =  d{(1/2)[1  (15)1  (14)  +  I  (i5)I  (1'4)]  -  21  (15)1  (14)} 
U  -r  —  —  +  ZZ 

(2.5.25) 

where  D  =  Y(15)7(l'^)-h  7r  .   We  recognize  the  B  and  A  terms,  respec- 
tively, of  Eqs.  (2.2.S)  with  6^0,  and  note  all  the  other  terras  vanish 
due  to  a  sinO  f.actor.   We  recall  the  angle  6  in  Sec.  2.3  was  between 
r  and  the.  z   axis  which  is  taken  to  be  zero  here. 

We  are  interested  only  in  the  first-order  energy  shifts  due  to 

?/p.(l-5-14),  given  by  the  diagonal  matrix  elements  <E,   +.I??'  (15-14)  ! £  +>. 

^  .  n-  D         n— 
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Expressing  Eq,  (2.5.20)  in  matrix  form. 


7/j^(15-14)  =  D 


_1_ 

0  1  o" 

0  0 

+ 1 

0  0  o" 

0   l' 

1  0    o" 

1    o" 

0  0   1 
^0  0  0 

1  0 

10  0 
0   10 

0  0 

0  0     0 

0  0  -1 

0  -1. 

Diagonal  elements  can  be  calculated 


(2.5.26) 


(a   a   a  1  ) 
In,  On,  -In 


0  1  0 
0  0  1 
0  0  0 


r^,    1 

In 

('«       «>         ''  - 

0  0 

^1± 

^On 

(8^+82+) 

1  0 

\± 

ci    , 

\. 

L  -Inj 

V.      f 


=  ^^m^'On-^  V^-ln>(^l±22±> 


(2.5.27) 


ind  so  on  for  the  other  terms.   The  result  is 


^ 


AE  .  =  D 
ni 


^^«Gn°'ln  -^  -0n"-ln>(^l±^2±>  "  ^^n  "  ^-ln>^^l±  "  A^ 


=  ±D 


(a^  a,   +  a^  a  ,  ) (ot,   -  a  i  )cos9 

On  In    On  -in    rx  In    -In 


'] 


=  -AE  _ 
n-s- 


(2.5.28) 


15. 


The   K  +%  >  -%  transition  will  occar,  to  zero  order,  at 

hw  =  YtiH  which  is  of  no  interest  and  will  be  shifted  by  AE 
c    '   o  n-*- 


-  AE   = 
n- 


2AE  , .   We  therefore  calculate  oaly 
n+ 


2aE 


n+ 


D  I  .''2'(a, 


On'^ln  ^   "Cn^-ln>^"^  '  ^^"L  "  ^-In^^""^ 


] 


(2.5.29) 
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It  is  useful  to  examine  the  6  dependence  of  2AE  , .   We  note  from 

n+ 

Eqs.  (2.5.11)  and  (2.5.12)  that  a   and  a   have  the  same  value  at  9 

and  180  -  6,  and  thus  4;  and  subsequently  T_  ,  T  ,  T  also  have  the 

same  value.   Therefore,  from  Eqs.  (2.5.16),  C,  (6)  =  C  ,  (180  -  6), 

in       —in 

C„  (9)  =  C„  (180  -  9)  and  C  (9)  =  C  (180  -  9).   Using  these  equations. 
On       On  n       n 

we  find  from  Eqs.  (2.5.17)  that  a   (9)  =  a_      (130  -  9)  and 

a_  (9)  =  a„  (180  -  9).   We  now  inspect  Eq.  (2.5.29)  and  note  (a„  a,  + 
Un       Un  Un  In 

ot-  a  ,  )  and  sin9  have  the  same  value  at  9  as  at  180  -  9,  while  both 
On  -In 

2      2 
(a.   -  Qt  ,  )  and  cos9  reverse  signs,  so  the  product  does  not  change 
In    -in 

value  or  sign.   Thus  2AE   needs  to  be  calculated  for  the  first  quad- 
rant only,  as  it  merely  repeats  the  same  values  in  the  second  quadrant. 

As  the  sample  was  polycrystalline,  we  again  must  work  out  the 
pov7der  pattern  spectrum.   In  Sec.  2.2  it  was  shown  that  any  value  of 
cos9  was  equally  probable  so  the  unbroadened  spectrum  was  the  'density' 
of  u(=  cos9)  versus  h,  the  dipole  contribution  to  the  field,  or 
du/dh.   It  was  easy  to  express  u  in  terms  of  h,  and  thereby  obtain  an 
analytic  expression. 

In  principle  the  same  procedure  would  work  here,  but  it  is  not 
possible  to  express  u  in  terms  of v  (equivalent  to  h  in  Sec.  2.2).   It 
is  possible  in  both  cases,  however,  to  numerically  calculate  v  (or  h) 
for  various  9's  and  weigh  the  results  by  a  factor  depending  on  the 
probability  of  each  9  occurring. 

We  center  a  sphere  on  the  origin  and  note  that  a  given  9.  and  A0 
define  an  annular  area  on  the  surface  of  the  sphere  given  by 


AA.  =  (lergth)  (width)  =  (27Trsin8 .)  (rAS) 
-«-  X 

=  (consta'-t)  sine  (2.5.30) 


-41- 


in  the  limit  of  snail  AS.   The  intensity  of  the  signal  at  frequency  v 

due  to  molecules  v/ith  orientation  9.  is  proportional  to  the  number  of 

molecules  with  that  orientation,  or  sin8 . .   We  therefore  plot  sin9 

1  i 

versus  v    to  obtain  the  unbroadened  spectrum,  or  convolute  with  a 

o  . 

1 

Gaussian  function  to  obtain  the  iatermolecular  broadened  function. 

A  computer  is  used  to  perform  the  calculations,  taking  values  of 
8.  between  0°  and  90°  with  a  sufficiently  small  A0  to  adequately  repre- 
sent the  smooth  curve. 


2.6  Nuclear  Spin-Lattice  Relaxation 
2.6.1  General 

The  resonant  condition  is  detected  by  noting  t)ie  energy  absorbed 
in  driving  nuclear  spins  to  a  higher  energy  state  or  by  noting  th^ 
energy  given  off  when  they  return  to  a  lower  energy  state.   The  radio 
frequency  field  K  ,  discussed  in  Sec.  2.1,  is  equally  likely  to  induce 
transitions  either  way  which  would  result  in  zero  net  energy  exchange 
if  there  were  the  same  number  of  spins  in  the  two  states  involved. 

In  practice,  the  population  of  the  lower  energy  state  is  slightly 
favored  according  to  the  Soltzmann  distribution  if  the  system  is  in 
thermal  equilibrium.   Therefore,  a  resonant  rf  signal  will  give  up  a 
net  amount  of  energy  to  the  system  as  long  as  there  are  more  spins  in 
the  lower  energy  state. 

It  can  be  seen  that  the  rf  signal  will  tend  to  destroy  the  popula- 
tion difference  by  inducing  transition.^  from  the  more  highly  populated 
state  at  a  faster  rate.   The  experimental  fact  that  the  system  returns 
to  the  equilibrium  distribution,  perhaps  at  a  rate  fast  enough  to  over- 
come the  tendency  of  the  rf  signal  to  destroy  the  equilibrium  condition. 
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is  evidence  that  there  must  be  some  other  mechanism  which  couples  with 
the  spins. 

We  find  there  are  several  mechanisms  by  which  the  spin  system  can 
give  up  energy  to  the  lattice  and  thereby  return  to  equilibrium.   The 
rate  at  which  this  occurs  is  termed  the  relaxation  rate  and  its  recipro- 
cal is  the  spin-lattice  relaxation  time  T  . 

Taking  the  simple  condition  of  a  spin  system  in  a  magnetic  field 

H  where  H   is  parallel  to  the  7.   axis,  we  see  the  population  difference 

results  in  a  net  magnetization  M  parallel  to  H  which  has  some  equili— 

z  o 

brium  value  M  .   The  magne':ization  M  approaches  M  at  a  rate  propor- 
tional to  the  difference  M  -  M  ,  resulting  in  the  exponential  equation 

M  (t)  =  M  (1  -  e"*"^'^!)   ,  (2.6.1) 

z       o 

where  t  is  measured  from  the  time  at  which  the  spin  order  was  completely 

destroyed  (M  (0^  =  0). 
z 

The  spin  system  exchanges  energy  tirilth  the  lattice  in  several  ways. 
For  instance,  in  Sec.  2.2  it  was  noted  that  the  magnetic  moment  of  one 
nucleus  produces  a  magnetic  field  at  the  sites  of  neighboring  nuclei. 
If  the  nucleus  is  undergoing  seme  type  of  raction  within  the  lattice,  it 
will  produce  a  changing  magnetic  field  at  other  sites.   We  might  imagine 
this  varying  field  as  a  sum  of  its  Fourier  components.   If  some  of  the 
components  are  of  the  appropriate  frequency,,  they  can  induce  transitions 
in  neighboring  spins,  allowing  the  nonequilibrium  spin  system  to  return 
to  equilibrium.   This  relaxation  towards  equilibrium  through  a  coupling 
of  the  spins  with  the  lattice  through  the  field  produced  by  nearby 
nuclear  dipoles  is  appropriately  called  dipolar  relaxation. 

A  nucleus  having  a  nuclear  quadrupole  moment  in  an  electric  field 
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gradient  which  is  time  dependent  due  to  motion  within  the  lattice  exper- 
iences nuclear  quadrupole  relaxation.   Considering  the  magnetic  raoment 
of  a  molecule  as  a  whole,  a  time-dependent  field  is  produced  at  the 
nuclear  sites  within  the  molecule  itself  due  to  the  rotational  motion 
of  the  molecule  which  produces  spin-rotational  relaxation. 

There  are  many  other  relaxation  processes  which  are  of  no  interest 
to  the  present  experiment.   It  might  be  that  several  processes  simultane- 
ously  contribute  to  the  relaxation  and  the  observed  relaxation  rate  is 
the  sum  of  the  rates  due  to  individual  processes, 

1/T^  =  1/T^^  +  1/T^Q  -M/T^3^  +  ...   .         (2.6.2) 
v/here  T   ,  T   ,  T    are  the  dipolar,  quadrupolar,  and  spin-rotation 

J.  lJ  LW      -Lot 

relaxation  times,  respectively.   It  may  be  relatively  difficult  in  such 
a  case  to  determine  the  relaxation  times  individually. 

The  molecular  motion  responsible  for  producing  the  fluctuating 
electric  or  magnetic  fields  at  the  site  of  a  nucleus  must  have  frequency 
components  in  the  spectral  density  at  resonance  and  must  be  significant 
enough  to   induce  transitions.   In  the  case  of  dipolar  and  quadrupolar 
relaxation,  a  reorientation  on  the  order  of  one  radian  at  the  appropri- 
ate rate  is  required,  whereas  spin-rotation  relaxation  Involves  only  a 
change  in  angular  momentum.   In  order  to  relate  this  molecular  motion 
to  the  temperature  of  the  sample,  we  use  a  quantity  called  tlie  correla- 
tion time  which  is  something  like  the  time  required  by  the  molecule  to 
undergo  the  appropriate  change  in  orientation  or  angular  momentum.   The 
quadrupolar  correlation  time  :,,  is  given  [l3]  by  the  integral  cf  the 
autocorrelation  function  of  P  (cosG)  where  9  specifies  the  orientation 
of  the  molecule,  and  the  spin-rota:ional  correlation  time  t_^   it.  given 
[13]  by  the  integral  of  the  timi-:  autocorrelation  function  of  molecular 
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Figure  4.   Temperature  dependence  of  quadrupolar  relaxation  time. 

(a)  Spectral  density  of  correlation  time  versus  frequency. 

(b)  Qualitative  dependence  of  Tj^q  versus  temperature.   Fre- 
quency independence  of  T-]^-,  at  high  temperature  is  a  result  of 
the  flatness  of  the  high  temperature  curve  in  the  upper  figure 
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angular  momentum. 

As  a  sample  cools,  collisions  are  more  frequent  due  to  increased 
density.   A  molecule  therefore  requires  more  time  to  change  orientation 
through  a  large  angle  as  the  random  collisions  cancel  each  other  out  to 
a  greater  extent.   The  correlation  time  t   therefore  increases  with 
decreasing  temperature.   On  the  other  hand,  a  single  collision  can  even 
completely  reverse  the  angular  momentum  of  a  molectile,  and  therefore, 
increasing  the  frequency  of  collisions  reduces  the  spin-rotational  cor- 
relation time  r 

sr 

The  reciprocal  of  the  correlation  time  is  something  like  the  fre- 
quency of  the  changing  electric  or  magnetic  field.   If  this  frequency  is 
very  high  so  cu  r  <<  1  where  u)   is  the  Larmor  frequency  of  the  nucleus, 
there  will  be  frequency  components  over  a  very  ^vrLde  range,  resulting  in 
only  a  small  proportion  near  w  .   The  mechanism  will  be  fairly  ineffec- 
tive under  these  conditions.   On  the  other  hand,  if  the  correlation  time 
is  very  long,  it  corresponds  to  a  very  low  frequency,  and  (jj  x  »  1. 
Most  of  the  frequency  components  x^ill  be  very  low,  and  again  the  mechan- 
ism will  be  fairlv  ineffective.   Between  the  two  extremes,  where  w  t  '{j  1, 

o 

the  frequency  is  on  the  order  of  the  Larmor  frequency  and  tlie  relaxation 
mechanism  attains  maximum  effectiveness  and  T  is  a  minimum. 

Figure  4  illustrates  the  cas3  of  quadrupolsr  relaxation.   Figure 
4(a)  shows  the  spectral  density  of  the  changing  electric  field  gradient 
at  both  high  and  low  temperatures.   A  iraximum  at  w  would  occur  at  some 
intermediate  temperature.   We  note  that  if  u   is  changed  by  changing 
the  external  field,  for  instance,  there  will  bt  little  effect  in  the 
high  temperature  region,  but  in  the  low  temperature  region,  reducing  u.« 
wil]  increase  the  effectiveness  and  shorten  T, , 
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Converted  to  T  dependence  on  temperature,  the  plot  would  be  some- 
thing like  Fig.  '^(b). 

2.6.2   The  Hubbard  Relation 

We  now  restrict  the  discussion  to  the  high  temperature  range  where 

0)  T   <<  1  and  where  T,„  is  independent  of  frequency  w  .   In  this  region 
o  Q  IQ  ° 

the  spin-rotational  component  of  the  relaxation  time  is  given  by  [14] 

tT^  =  (4/3)  (I  kT/H^)(2TTC)^T^^  (2.6.3) 

Isr  o  sr 

where  I   is  the  molecular  moment  of  inertia,  T  is  the  temperature,  and 
o 

C  is  the  spin-rotational  constant. 

The  quadrupolar  relaxation  time  is  given  by  [4] 

T^J  =  (3/8)(e^qQ/il)^Q   .  (2.6.4) 

We  now  look  at  the  relation  between  t^^  and  t  .   If  t_^^  is  long 
enough  to  allow  the  molecule  to  undergo  large  angle  reorientation 
between  collisions,  we  have  t   <^  t^^.   Furthermore,  changing  t^^  does 
not  directly  have,  much  effect  on  t  .   On  the  other  hand,  we  may  inciease 
the  frequency  of  collisions  and  chereby  reduce  x ^^   to  the  point  where 
many  random  collisions  are  required  before  the  molecule  reorients 
through  a  large  angle.   In  this  region  x   >>  x^^  ar/d  a  change  in  t^^ 
will  be  accompanied  by  an  opposite  change  in  x  .   This  compleraeatary 
effect  has  been  shown  to  obey  the  cord it ion 

XX   =  t  /6kT  (2.6.5) 

Q  sr    o 

by  Hubbard  [14]  using  a  rotational  diffusion  model.   The  Hubbard  rela- 
tion can  be  written 

x*x*  =  1/6   .  (2.6.6) 

Q  sr 
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U 

where  reduced  correlation  times  are  defined  by  t*  =  T(kT/I  ) *.   The 

o 

reduced  correlation  times  are  the  correlacion  times  expressed  in  units 
of  time  required  for  a  classical  freely  rotating  spherical  molecule  to 
rotate  through  one  radian.   The  Hubbard  relation  was  based  on  spherical 
molecules;  hov;ever  Kluk  and  Powles  [l5]  have  shown  it  holds  for  linear 
molecules  as  well. 


CHAPTER  III 

EXPERIMENTAL  EQUIPMENT  AND   PROCEDURE 


3.1   Superconductive  Magnet  System 

3.1.1  Magnet 

The  heart  of  the  system  is  an  RCA  Superconductive  Magnet  Type 
SM-28A1  capable  of  producing  a  maximum  field  strength  of  100  kG  at  a 
current  level  of  91  A.   The  solenoid  itself  is  about  7  in.  long  and  7 
in.  in  diameter  with  a  bore  of  just  over  two  inches  and  is  wound  of 
Nb  Sn  superconductive  ribbon.   Overall  dimensions,  including  the  form 
(housing),  but  excluding  mounting  studs,  are  10  in.  long,  7.2  in.  out- 
side diameter,  and  a  2.03  in,  bore. 

The  magnet  is  equipped  with  a  magneto-resistive  probe  for  field 
strength  measurement,  which  also  turns  out  to  be  useful  as  a  thermometer 
during  cool-down.   With  50  mA  dc  through  the  magneto-resistive  probe  the 
voltage  varied  from  about  700  mV  at  room  temperature  to  84  mV  at  liquid 
nitrogen  temperature  and  2.84  mV  at  liquid  helium  temperature  at  zero 
field.   Thus  it  was  easy  to  determine  when  the  magnet  had  reached  liquid 
nitrop.en  temperature   during   precooling  and  when  liquid  helium  began 
collecting  during  helium  transfer. 

At  liquid  helium  temperature,  again  with  50  mA  dc,  the  voltage 
across  the  magneto-resistive  probe  varied  from  2.84  mV  at  zero  field  to 
about  20  mV  at  100  kC.   Once  calibrated,  it  was  accurate  to  better  than 
±20  G. 
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Tlie  magnet  is  equipped  with  a  superconductive  switch  across  the 
terminals  for   scmipersistent  operation.   A  300-mA  current  through  a 
resistance  wire  in  the  switch  keeps  the  switch  in  its  normal  state. 
Shutting  off  the  heater  current:  allows  the  switch  to  cool  and  become 
superconductive.   The  superconductive  path  is  then  complete  except  for 
the  resistance  of  the  copper  current  cenninais  and  contact  resistance 
where  the  ends  of  the  coil  and  the  ends  of  the  switch  are  soldered  to 
the  current  terminals. 

On  one  occasion  the  switch  heater  wire  was  accidentally  burned  out, 
making  it  impossible  to  charge  the  solenoid,  so  the  epoxy-potted  switch 
had  to  be  disassembled  and  repaired.   Unfortunately,  this  happened 
before  the  first  resonance  was  obtained  so  behavior  of  the  magnet  with 
the  original  switch  as  delivered  was  not  determined. 

The  magnet  was  advertised  to  have  a  homogeneity  of  5  G  over  a  3/4- 
in.  diameter  spherical  volume  (DSV) .   In  order  to  test  the  field  char- 
acteristics and  to  have  a  convenient  nitrogen  reference  signal,  a 
liquid  ammonia  vial  was  prepared  which  had  inside  dimensions  of  about 
7/15  in.  long  and  3/16  in.  diameter.   The  magnet  was  mounted  with  the 

bo^e  (and  hence  H  )  vertical  and  the  ammonia  sample  was  mounted  horizon- 
o 

tally. 

Variation  of  the  field  and  line  width  with  position  of  the  sample 
along  the  axis  of  the  solenoid  are  shown  in  Fig.  5.   The  abscissa  is  in 
inches  above  an  arbitrary  reference  point  on  the  axis  near  the  center 
of  the  solenoid.   The  ordinates  of  the  line-width  data  points  are  in 
gauss,  and  were  determined  by  peak-to-peak  measurements  of  the  deriva- 
tive signal  obtained  by  a  cw  spectrometer.   The  ordinates  of  the  mag- 
netic field  strength  H   are  in  gauss  above  an  arbitrary  L-efcreace  of 
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Figure  6.   Differences  between  measured  values  of  magnetic  field  and 

value  calculated  by  exponential  decay  equation.   (a)  Decay  to 
zero  field  assumed.   Smooth  curves  are  drawn  throu^'^h  data 
points  at  0.50,  0.55,  and  0.60  inches  from  Fig.  5.   (b)  Decay 
to  intermediate  value. 
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about  89  kG,  and  represent  the  zero  crossover  of  a  sometimes  quite  wide 
derivative  line. 

The  data  points  had  to  be  corrected  for  a  varying  field  drift  rate, 
the  spectrometer  frequency  sweep  rate,  direction  of  sweep,  and  the 
spectrometer  time  constant  used. 

As  the  field  measurements  are  an  average  over  the  sample  volume, 
they  cannot  be  taken  directly  as  a  point-to-point  variation  along  the 
axis  of  the  bore.   Considering  the  sample  size,  however,  it  should  be 
possible  to  move  it  about  0.45  in.  and  still  keep  it  entirely  within 
the  3/4-in.  DSV  which  was  advertized  to  have  a  homogeneity  of  5  G.   No 
such  region  was  found.   lu  appears  the  homogeneity  was  closer  to  7-9  G 
over  the  3/4-in.  DSV. 

The  peak-to-peak  amplitude  of  the  c'e.rivative  signal  was  roughly 
inversely  proportional  to  the  l:.ne  widuh  over  most  of  the  region,  ^s 
expected.  However,  tne  maxira-im  f.,n';lJtude  occurred  at  0.9  in.  on  the. 
reference  scale  and  unexpectedly  decreased  rapidly  as  the  sample  was 
raised  until  the  signal  was  hardly  discernible  above  1.05  in. 

The  0.9 -in.  position  v;as  chosen  for  the  center  of  the  field  as  the 
best  combination  of  linewidth  and  signal  amplitude. 

The  field  typically  decayed  at  a  rate  of  0.3  to  0.9  G/min.  in  the 
semipersistont  mode.   Values  of  the  field  recorded  with  the  sample  at 
the  0.55-in.  position  were  first  crudely  fitted  to  the  exponential  decay 
equation 

H(t)  =  H(0)exp(-t/T)  (3.1.1) 

by  a  least  squares  method.   Differences  between  the  data  points  and  the 
values  calculated  according  to  Eq.  (3.1.1),  H,(t)  -  H(t),  are  plotted 
in  Fig.  6(a)  as  solid  circles. 
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As  the  value  of  the  field  depends  on  both  drift  rate  and  variation 
of  the  field  with  position  of  the  sample,  an  error  in  the  assumed  drift 
rate  causes  an  erroneous  determination  of  field  variation  with  position. 
This  effect  is  apparent  in  Fig.  6(a)  where  data  points  at  positions 
other  than  0.55  in.  are  shown  as  open  circles  which  are  offset  from  a 
smooth  curve  through  the  solid  circles. 

Clearly  Eq.  (3.1.1)  is  not  adequate  over  the  11.5-h  duration  of 
the  experiment.   Hence,  the  equation 

H(t)  =  H(-)  +  [H(0)  -  H(=°)]exp(-t/T)  (3.1.2) 

was  used.   A  plot  of  H.(t)  -  H(t)  versus  time  is  shown  in  Fig.  6(b) 
where  it  is  seen  the  fit  is  far  better  and  entirely  adequate. 
Equation  (3.1.1)  was 

H(t)  =  89,411  exp(-t/l/55.5)  (3,1.3) 

and  Eq.  (3.1.2)  was 

H(t)  =  81,774  +  7639.1  exp(-t/144.5)         (3.1.4) 
where  t  is  in  hours  and  II  in  gauss. 

Some  variation  in  the  values  of  K(">)  and  x  occurred  on  succeeding 
runs,  but  generally  H(°°)  could  be  considered  constant  during  any  given 
run.   On  occasion  the  drift  late  "^a3  nearly  double  the  typical  0.8  or 
0.9  G/min  value.   During  one  run,  the  drif';  rate  increased  with  time  as 
though  H(°°)  were  decreasing. 

Maintaining  full  current  in  the  le^ds  appeared  to  have  no  effect 
on  the  drift  rate,  but  the  magnet  was  not  operated  for  extended  periods 
with  full  current  in  the  leads  due  to  increased  liquid  helium  consump- 
tion. 

In  some  runs,  the  field  V7as  deliberately  increased  beyond  the 
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desired  value,  then  decreased.  The  drift  rate  was  affected  only  the 
first  few  minutes  before  returning  to  a  behavior  similar  to  that  de- 
scribed by  Eq.  (3.1.2). 

It  appears  that  Eq.  (3.1.2)  is  normally  quite  accurate  over  opera- 
ting periods  of  up  to  23  h  used  thus  far.   At  some  time  the  value  of 
H("»)  must  decrease  due  to  the  small  resistance  in  the  solenoid  circuit, 
but  no  long  term  study  of  this  effect  was  undertaken. 

Several  field  determinations  during  each  run  using  the  ammonia 
sample  were  adequate  for  drift  determination  and  would  probably  still 
be  necessary  even  if  the  long  term  behavior  were  known. 

When  the  magnet  was  first  installed,  it  frequently  quenched  for 
reasons  which  were  not  initially  understood.   As  experience  was  gained, 
some  causes  were  identified  and  eliminated.   One  cause  of  quenching 
which  could  not  be  eliminated  is  thought  to  be  connected  with  the  re- 
paired superconductive  switch.   The  current  through  the  switch  is  the 
difference  in  the  solenoid  current  and  the  current  in  the  leads  to  the 
power  supply.   In  the  semipersistent  mode  the  switch  current  therefore 
increases  as  the  power  supply  current  is  reduced.   The  magnet  quenches 
fairly  consistently  if  the  switch  current  becomes  greater  than  about 

50  A.   This  can  be  explained  in  terms  of  resistance  in  the  switch  or 

-4 
switch  solder  connections.   A  resistance  of  10   fi  and  50  or  50  A  pro- 
duces about  as  much  heat  as  the  300  m^.  in  ihe  heater  wire  used  to  keep 

2 
the  switch  normal.   Thus  the  I^R  heat  associated  with  a  small  resistance 

at  the  solder  joints  may  be  enougii  to  drive  the  supercondacLive  switch 
normal  even  though  it  is  in  contact  with  large  copper  heat  sinks. 

The  result  was  that  the  power  supply  had  to  be  left  engaged  carry- 
ing about  30  A  to  keep  the  current  in  the  suierconductive  switch  below 
50  A.   This  of  course  increased  the  helium  ccnsumptiou  and  possibly 
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Figure  7.   Top  access  dewar  system. 
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contributed  to  electronic  noise.   It  also  made  the  system  susceptible 
to  pov/er  failure  which  would  not  be  conducive  to  long  temi  operation. 

3.1.2   Dewars 


Two  interchangeable  sets  of  dewars  were  obtained  for  use  witli  the 
superconductive  magnet.   The  first  set,  manufactured  by  Sulfrian 
Cryogenics,  Inc.,  consisted  of  separate  nitrogen  and  helium  dewars  with 
top  access.   It  was  complemented  by  a  Janis  Research  Company  Super  Vari- 
Temp  insert  assembly  which  supported  the  magnet,  provided  access  to  the 
working  space,  and  provided  sample  temperature  control.   The  system  is 
shown  schematically  in  Fig.  7. 

Tlie  other  system  consisted  of  a  Ci;yofab  unit  construction  nitrogen 
and  helium  dewar  with  a  coiranon  vacuum  space  and  equipped  witli  a  re- 
entrant  access  tail.   It  is  schematically  shovrn  in  Fig.  3. 

The  first  system  has  the  advantage  of  low  tempeiature  opei"ation 
and  built-in  temperature  control.  Liquid  helium  flow  through  the  capil- 
lary in  conjunction  with  either  electrical  heaters  in  the  bottom  of  the 
insert  or  a  heater  added  to  the  sample  holder  permits  operation  from  1.5 
to  300  K.  The  obvious  disadvantage  is  the  distance  fron.  the  top  access 
flange  to  the  working  center  of  the  magnet,  about  60  in. 

Lack  of  adequate  vertical  clearance  in  the  laboratory  made  it 
necessary  to  either  build  the  sample  probe  in  sections  or  to  shut  the 
system  dovm  and  lower  it  through  the  floor  in  order  to  insert  the 
sample.   Structural  steel  in  the  concrete  floor  prevented  operation  in 
the  lowered  position. 

The  re-entrant  dewar  had  the  significant  advantage  of  easy  sample 
access.   Tlie  center  of  the  magnetic  field  was  only  about  11  in.  up  from 
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the  outside  bottom  of  the  dewar  and  thus  tht:  sample  probe  could  be  in- 
serted or  removed  while  the  magnet  remained  at  high  field.   Temperature 
control  was  less  convenient,  partly  because  it  was  not  built  in  and 
partly  because  of  the  short  stand-off  distance  frorp  room  temperature. 

A  gas-flow  temperature  control  system  to  be  described  later  used 
in  conjunction  with  the  re-entrant  dewar  permitted  operation  over  a  90 
to  44  K  sample  temperature  range.   The  desired  sample  temperature  could 
be  attained  and  stabilized  fairly  rapidly  but  with  increasing  diffi- 
culty near  the  extremes  of  the  range. 

One  problem  required  constant  attention  during  operation.   The 
helium  space  in  the  dewars  was  connected  to  the  helium  recovery  system 
through  the  current  lead  feed  throughs  and  also  through  a  larger  port. 
Severe  thermal  oscillations  occurred  occasionally  which  could  be  stopped 
by  a  trial-and-error  adjustment  of  valves  in  the  recovery  lines.   In 
general,  the  main  valve  in  the  recovery  line  from  the  large  port  in  the 
dewar  could  not  be  opened  very  wide  without  inducing  these  oscillations. 

The  dewars  were  suspended  on  a  large  steel  plate  which  could  be 
raised  and  lowered  on  a  vertical  track  by  means  of  a  cable  and  winch 
arrangement.   The  dewars  alone  cr  the  whole  assembly  could  be  lowered 
either  to  gain  access  to  the  magnet  or  to  obtain  adequate  overhead 
clearance  for  sample  probe  installation  in  the  top  access  dewars. 

3.1.3   S^i-mple  Probe 

A  sample  probe,  shown  in  Figs.  9  and  10,  was  constructed  to  posi- 
tion and  orient  the  sample  as  well  as  to  provide  temperature  control 
and  rf  connections. 

The  outer  conductor  of  the  rf  coaxial  cable  doubled  as  structural 
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Figure  9.   Sample  probe,  side  view  of  sample  holder. 
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support.   The  crystal  was  mounted  inside  a  cylindrical  cage-like  holder 
which  had  worm  gear  teeth  cut  around  its  center.   The  angular  rotation 
of  the  crystal  could  be  controlled  to  within  +2  minutes  of  arc  relative 
to  the  initial  orientation,  and  the  initial  orientation  could  usually 
be  accurately  determined  from  the  data.   Slack  in  leads  from  the  coil 
to  the  top  of  the  solid  rf  coax  permitted  a  full  180°  crystal  rotation. 

The  entire  unit  fit  snugly  in  the  re-entrant  tail  of  the  dewar, 
leaving  a  small  space  above  the  nylon  head.   Hot  (or  cold)  gas  is  passed 
up  to  this  space  through  a  3/16~in.  stainless  steel  tube  which  was 
■"acuum  jacketed  up  to  the  top  spacer.   The  gas  flows  back  through  and 
around  the  nylon  head  to  the  space  below,  where  it  exhausts  through  a 
1/4 -in.  vacuum  jacketed  tube. 

A  uniform  temperature  is  maintained  in  the  space  above  the  top 
Teflon  spacer  by  a  combination  of  gas  flow  rate  and  insulation.   Near 
room  temperature  there  is  only  a  small  temperature  gradient  between  the 
sample  area  and  ambient  temperature  outside  the  dewar,  providing  little 
heat  loss.   At  higher  temperatures,  a  much  higher  hot  gas  flow  rate  is 
used  which  offsets  the  greater  heat  loss  through  the  insulation  withcut 
significantly  cooling  the  gas.   A  similar  situation  appJies  at  low 
temperatures, 

A  thermocouple  close  to  the  sample  was  used  to  measure  temperature 
and  was  adequate  in  view  of  the  weak  temperature  dependence  encountered. 

The  gas-flow  systems  for  operation  both  above  and  below  room  tem- 
perature are  shoivn  in  Fig.  11.   For  high  temperature  operation,  air  was 
passed  through  a  pressure  regulator,  desiccator,  and  electrical  heat 
exchanger.   Manual  adjustment  of  inlet  air  pressure  and  heat  exchanger 
voltage  provided  better  tlian  1  K  stability  with  reasonably  infrequent: 
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adjustments. 

The  low  temperature  system  operated  on  cold  nitrogen  gas  boiled 
off  by  passing  current  through  a  25-^  resistance  immersed  in  liquid 
nitrogen.   Stability  was  enhanced  by  passing  the  gas  through  copper 
tubing  wound  around  a  copper  block  which  provided  thermal  inertia.  Tem- 
perature variation  was  predominantly  provided  by  varying  the  current 
through  the  25-Q  resistance  in  the  liquid  nitrogen  supply,  with  addi- 
tional control  achieved  by  means  of  a  resistance-wire  heater  on  the 
copper  block. 

The  systems  were  operated  over  a  range  from  90  to  440  K.   Higher 
temperatures  would  have  been  readily  obtainable  with  no  modification  to 
the  system,  the  limit  being  t'ne  ability  of  the  sample  probe  to  withstand 
the  heat.   Initially,  an  rf  coaxial  cable  was  inserted  in  the  probe  all 
the  v;ay  up  to  just  below  the  sample  holder,  but  plastic  spacers  inside 
the  cable  melted  at  167°C.   Tha  nylon  worm  gtar  assembly  also  deformed 
and  seized  at  high  temperatures.   These  problems  were  corrected  by  use 
of  rigid  rf  conductors  as  shown  in  Fig.  9  and  by  slightly  increasing 
clearances  between  moving  parts. 

Somewhat  lower  temperatures  should  be  possible  with  some  modifica- 
tion to  the  system.   The  inlet  gas  tube  could  be  passed  through  the 
liquid  nitrogen  to  insure  it  is  as  cold  as  possible  prior  to  entering 
the  vacuum  transfer  tubes.   Also,  pumping  on  the  exhaust  tube  to  in- 
crease gas  flow  rather  than  passing  current  through  resistors  in  liquid 
nitrogen  would  have  the  additional  effect  cf  depressing  the  boiling 
point. 
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3 • 1.4   Operation  of  the  Superconductive  System 

It  was  generally  best  to  start  cool-dowi  about  three  days  before 
the  actual  run.   Simply  maintaining  liquid  nitrogen  in  the  liquid  nitro- 
gen dewar  v/ould  cause  the  magnet  temperature  to  lower  to  an  estimated 
81-83  K  in  that  time  period  and  an  extra  day  would  result  in  lowering 
the  temperature  2  to  3  degrees. 

It  V7as  also  possible  to   precool  rapidly  by  introduction  of  liquid 
nitrogen  into  the  helium  dewar  Ln  direct  contact  with  the  magnet,  but 
this  method  required  back  transfer  which  was  more  trouble  than  keeping 
the  outer  jacket  full  of  liquid  nitrogen. 

Overhead  clearance  limited  the  length  of  the  helium  transfer  tube 
so  when  fully  inserted  it  would  reach  roughly  one  third  of  the  way  down 
from  the  top.   Thus,  a  vacuum  jacketed  extension  was  required,  but  it 
was  also  desirable   to  be  able  uo  pull  the  extension  up  out  of  the 
liquid  once  liquid  helium  was  in  the  dewar.   This  was  especially  impor- 
tant to  avoid  passing  warm  gas  through  the  liquid  when  starting  a 
replenishment  transfer.   Therefore,  a  vacuum  jacketed  extension  was 
made  which  extended  through  the  r.op  flange  of  the  dewar  where  it  was 
terminated  with  a  ball  valve. 

This  extension  could  be  lowered  for  insertion  of  the  transfer  tube 
and  for  initial  transfer  when  it  was  desirable   to  have  the  cold  gas 
pass  up  around  tjie  magnet  from  below.   For  replenishment,  it  could  be 
raised  part  way,  and  bet\;eeu  transfers  it  could  be  raised  well  above 
the  liquid  helium  to  reduce  the  heat  path. 

During  cool-down,  a  pressure  in  the  supply  dewar  of  about  5  in. 
H„0  above  recovery  line  pressure  vould  be  the  maximum  the  recovery 
system  would  take,  but  once  l:''quid  helium  started  collecting,  40  in. 
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H„0  pressure  or  more  could  be  used.   Only  a  few  minutes  would  be  re- 
quired to  collect  about  25  liters  after  a  45-rain  cool-down  from  liquid 
nitrogen  temperature. 

Once-a-day  replenishment  was  generally  adequate.   Helium  consump- 
tion was  about  2  liters/h  during  the  time  the  field  was  being  either 
increased  or  decreased  and  was  as  low  as  0.3  liter/h  with  the  system 
shut  down.   Vibrations  fed  back  from  the  recovery  compressor,  thermal 
oscillations,  changes  in  recovery  line  pressure,  current  in  the  electri- 
cal leads,  and  possibly  other  factors  affected  the  helium  consumption 
which  usually  averaged  about  20  liters/day  including  that  lost  while 
cooling  the  warm  transfer  tubes  during  replenishment  transfers. 

Moveable  level  resistors  in  both  the  experimental  and  supply  dewars 
made  it  very  easy  to  monitor  liquid  helium  levels.   The  rather  simple 
level  indicators  exhibited  only  a  small  needle  deflection  as  the  resis- 
tor passed  from  one  side  of  the  siirface  to  the  other  so  needle  movement 
while  dipping  the  resistor  in  and  out  of  liquid  helium  was  a  much  more 
positive  indication  than  absolute  needle  position. 

A  Magnion  CFC-100  power  supply  was  used  to  charge  the  solenoid.   A 
40-turn  wire-wound  potentiometer,  irotor  driven  at  a  selected  rate,  con- 
trolled the  rate  cf  current  increase.   A  maximum  rate  of  2.0  A/min  was 
allowed,  which  required  kt  .j   -ciu  to  bring  tl-e  magnet  up  to  100  kG.   It 
was  discovered  that  the  core  of  the  40- tuTri  potentiometer  had  a  high 
coefficient  of  thermal  expansion)  ^nd  ^'onld  sometimes  shrink  so  much 
that  the  wire  would  loosen  and  cause  intermittent  contact  with  the  wiper 
which  would  result  in  the  magnet  quenching.   The  40-turn  potentiometer 
had  to  be  replaced  with  one  having  only  Iv)  turns.   With  this  potentio- 
meter the  fastest  charging  rate  v.-itbcut  e::cecdln-2  2.0  A/min  was 
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1.6  A/min  which  increased  charging  time  by  25%. 

Prior  to  placing  the  magnet  in  the  semipersistent  mode,  the 
magneto-resistive  voltage  and  power  supply  dial  reading  were  recorded. 
The  power  supply  current  could  then  be  reduced  to  about  30  A  with  the 
magnet  in  the  semipersistent  mode.   At  the  end  of  the  run,  the  raagneto- 
resistive  voltage  was  again  noted  which  would  be  somewhat  lower  than 
the  initial  reading.   The  power  supply  current  would  be  brought  up  to 
a  proportionately  lower  reading  prior  to  opening  the  superconductive 
switch.   The  power  supply  current  and  the  magnet  current  would  be  close- 
ly matched  in  this  condition,  any  difference  between  them  being  through 
the  superconductive  switch.   When  the  switch  is  opened  by  driving  it  to 
the  normal  state,  the  current  in  it  causes  a  voltage  to  de\/elop  across 
it,  normally  less  than  25  mV.   In  case  of  too  great  a  mismatch  of 
currents,  the  voltage  would  be  high  enough  to  cause  the  power  supply  to 
quench  which  would  result  in  the  magnet  discharging  through  a  2-Q   shunt 
across  the  terminals.   Under  these  conditions  the  magnet  would  be  de- 
energized  in  a  few  minutes  with  little  loss  of  liquid  helium. 

As  mentioned  already,  various  causes  would  result  in  quenching  the 
magnet,  and  it  would  very  rapidly  de- energize,  generally  boiling  off 
all  the  liquid  helium  in  the  dewar  but  not  warming  up  appreciably.   The 
recovery  system  was  not  capable  of  handling  such  a  large  flow  rate,  so 
most  of  the  helium  would  be  lost  to  the  atmosphere  through  a  safety 
valve.   A  helium  gas  bag  is  needed  to  prevent  this  loss.   It  appears 
that  the  primary  remaining  cause  of  quenching  is  excessive  current  in 
the  superconductive  switch  which  c?n  be  avoided  by  leaving  the  power 
supply  on  and  carrying  part  of  the  current. 
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3.2      Elect:roni£3c;nat:   System 

3.2.1  Electromagnet 

A  Varian  V-4012-3B  12-inch  electromagnet  with  a  3.5-in.  gap  was 
used  for  the   N  line  shape  and  relaxation  studies.   The  magnet  was  capa- 
ble of  providing  magnetic  fields  of  up  to  9000  G.   A  Varian  V--2100B 
magnet  power  supply  and  a  Varian  V-FR.-2100  Fieldial  field  regulator  were 
used  to  operate  the  magnet  in  the  field-regulation  mode. 

Drift  rates  of  up  to  about  i  G/h  were  encountered  whereas  a  change 
of  as  little  as  0.03  G  v;ould  sometimes  cause  significant  inaccuracies 
in  the  data.   It  was  sometimes  possible  to  correct  for  the  drift  rate 
once  it  was  known;  however  the  drift  rate  varied  and  even  reversed  direc- 
tion throughout  the  day,  especially  when  the  laboratory  temperature  was 
fluctuating  duo  to  intermittent  operation  of  air  conditioning. 

3.2.2  Cryostat  and  Sample  System 

A  schematic  diagram  of  the  cryostat  is  shown  in  Fig.  12.   The 
sample  cell  was  made  of  Kel-F  and  had  a  volume  of  about  1.8  cm  .   The 
space  around  it  contained  helitim  exchange  gas  for  thermal  contact  with 
the  surrounding  copper  can  which  contained  the  heater  and  theiTnometer. 
Exchange  gas  in  the  brass  can  was  used  to  control  thermal  contact  with 
the  cryogenic  fluid  outside  the  brass  can.   The  sample  gas  down  pipe 
doubled  as  the  center  rf  conductor,  and  the  sample  assembly  was  sup- 
ported by  a  double  walled  tube  v/hich  also  served  as  the  outer  rf  con- 
ductor.  The  space  between  the  walls  was  evacuated  to  therroally  isolate 
the  sample  gas  dcwn  pipe  from  the  surrounding  cryogenic  fluid.   The 
sample  gas  down  pipe  was  electrically  insulated  from  the  external 
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Figure  12.   Cryostal  used  with  electromagnet. 
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Figure  13.   Sample-gas  system. 
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plumbing  by  an  arrangement  not  shown  on  the  drawing. 

Temperature  control  was  provided  for  by  a  GaAs  thermometer  and 
resistance  heater  wire  installed  on  the  copper  can  surrounding  the 
sample  cell  containing  exchange  gas,  and  by  varying  the  exchange-gas 
pressure  outside  the  copper  can  to  vary  thermal  contact  with  the  cryo- 
genic fluid. 

Figure  13  is  a  schematic  of  the  sample-gas  system.   The  sample  gas 
was  normally  stored  in  the  high  pressure  cylinder.   Prior  to  releasing 
the  sample  to  the  rest  of  the  system,  the  pressure  was  reduced  to  approxi 
mately  one  atmosphere  by  immersing  the  cylinder  in  liquid  nitrogen.  The 
gas  was  then  moved  around  by  controlling  the  temperature  of  the  sample 
cell,  the  high  pressure  cylinder,  and  the  two  cryopumps.   Ideally, 
thermal  isolation  of  the  sample  cell  and  down  pipe  would  have  perroitted 
operation  with  liquid  helium  as  a  cryogenic  fluid  over  a  wide  tempera- 
ture range.   In  practice,  if  helium  was  used  during  the  sample  conden- 
sation stage,  the  down  pipe  would  freeze  up  and  block  before  the  sample 
was  condensed  in  the  sample  cell.   Part  of  the  problem  was  due  to  an 
uninsulated  tube  which  passed  up  through  the  cryogenic  fluid  connecting 
the  copper  sample  cell  to  room  temperature.   This  provided  some  thermal 
contact  even  when  the  exchange  gas  was  removed  from  the  outside  of  the 
copper  can.   Rather  than  modifying  the  cryostat,  pumped  liquid  and  soli.d 
nitrogen  were  used  while  condensing  the  sample  and  for  operation  above 
50  K. 

The  sample-cell  temperature  would  be  maintained  at  about  65-67  K 
while  the  remainder  of  the  sample  system  was  gradually  wanned,  thus 
condensing  the  sample  in  the  sample  cell  as  it  boiled  off  from  other 
parts  of  the  sample  system.   The  pumped  liquid  or  solid  nitrogen  could 
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then  be  used  as  the  ^ryogenic  fluid  (or  solid  below  63  K)  or  it  could 
be  replaced  by  liquid  helium  for  jouec  temperature  operation. 

3 . 3  Spectrometers 

3.3. 1   Pulse  spectronete r 

Various  pulse  spectrometer  coiabinations  were  used  depending  on 
equipment  availability  and  the  type  of  measurements  to  be  made. 

It  was  normally  not  possible  to  observe  the  free  induction  decay 

14 
(FID)  of   N  in  a  glycine  directly  without  signal  averaging.   Thus  it 

was  necessary  to  adjust  the  frequency  by  small  increments  and  signal 

averaf,e  after  each  adjustment  to  first  detect  the  signal  and  then  to 

determine  the  frequency,  s  very  tedious  process.   Under  these  conditions, 

the  Fast  Fourier  Transforn  (FFT)  technique  was  employed.   Tlie  Fourier 

transform  of  the  FID  is  the  frequency  spectrum  or  absorption  line  with 

respect  to  the  oscillator  frequency  as  the  reference.   Thus,  if  the 

oscillator  frequency  is  displaced  below  the  central  frequency,  the  peak 

of  the  absorption  line  v/ill  be  offset  above  the  origin,  and  the  amount 

of  offset  is  added  to  the  oscillator  frequency  to  determine  the  central 

frequency.   Kence  it  is  only  necessary  to  set  the  oscillator  frequency 

close  to  the  frequency  of  the  signal. 

Complex  Fourier  trans fonriation  capability  would  provide  information 

a-3  to  whether  the  offset  it;   above  or  belov  the  oscillator  frequency. 

Equipi'jent  used  in  the  prei-i^nt  work  did  not  have  complex  capability,  and 

only  the  frequency  differences  could  be  foLind.   The  ambiguity  as  to 

whether  the  offset  is  above  or  below  the  oscillator  frequency  had  to  be 

resolved  by  other  means.   The.  usual  method  was  to  change  the  oscillator 
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frequency  and  note  whetlier  the  offset  increased  or  decreased. 

It  is  not  only  unnecessary  to  set  the  oscillator  frequency  to  the 
central  frequency  of  the  line,  it  is  undesirable  without  complex  Fourier 
transform  capability.   If  the  oscillator  frequency  is  within  the  fre- 
quency spectrum  of  the  line,  pacts  of  the  spectrum  equal  distances  above 
and  below  the  oscillator  frequency  will  be  added,  thus  distorting  the 
frequency  spectrum. 

The  spectrcoieter  used  in  ccajur\ction  with  the  FFT  equipment  had 
frequency  stabilization  provided  by  a  General  Radio  Model  1061  Frequency 
Syntliesizer .   It  could  be  set  digitally  in  10  Hz  increments  and  incor- 
porated a  search  mode  by  which  any  digit  position  could  be  replaced  by 
a  continuously  variable  value.   For  example,  replacing  the  10  Hz  digital 
control  permitted  a  continuous  variation  over  a  range  from  -20  Hz  to 
109  Hz  belov7  or  above  tiie  frequency  set  on  the  other  controls. 

An  F  and  H  Instruments  Pulse  Program  Generator  PPG  45  permitted 
considerable  flexibility  in  choice  of  pulse  sequence,  width,  separation, 
repetition  rate,  etc. 

A  locally  made  Quadrature  Phase  Sensitive  Detector  was  especially 
convenient  when  it  was  desired  to  monitor  the  tuning  while  recording 
data.   I'Jhen  the  frequency  and  phase  are  properly  adjusted,  the  FID  ampli- 
tude is  a  maximum  while  the  quadrature  oiitput  is  flat,  and  a  drift  off 
resonance  is  easier  to  detect  relative  to  a  flat  response.   This  of 
course  was  not  used  when  the  oscillator  v;as  deliberately  set  off  reso- 
nance. 

Each  FID  was  recorded  by  a  Bionation  Model  801  Transient  Recorder 
and  transferred  tc  a  Fabri-Tek  Instruments,  Inc.,  Model  1072  Instrument 
Computer  for  accuraulaCioa. 
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After  accumulation  of  an  adequate  number  of  signals,  a  Digital 

Equipment  Corporation  PDP  S/E  computer  was  used  to  perform  the  FFT. 

lA 
Uliile  observing  the  NQR-perturbed   N  Zeeman  resonance  in  glycine, 

14 
the  frequency  of  the  unperturbed   N  Zeeman  line  was  periodically  deter- 
mined by  means  of  a  liquid  ammonia  sample.   From  these  data  and  the 
calculated  drift  rate  of  the  superconductive  magnet,  the  unperturbed 
frequency  was  calculated  for  each  time  at  which  data  were  recorded. 

After  performing  the  FFT,  the  displacement  of  the  line  center  from 
the  origin  was  measured  in  frequency  units  and  the  correction  was 
applied  to  the  oscillator  frequency,  resulting  in  the  actual  frequency 

of  the  NQR-perturbed  line.   Finally  the  difference  between  this  fre- 

14 
quency  and  the  frequency  ot  the  unperturbed   N  Zeeman  line  was  deter- 
mined which  was  the  frequency  shift  of  interest. 

Wien  using  the  same  spectro-netar  to  measure  T   of   N  in  liquid  and 
6-solid   N  ,  th(2  amplitude  of  the  FID  was  sampled  o^.-er  a  short  section 
near  the  origin.   Successivf!  amplitudes  were  recorded  in  successive 
channels  of  the  Biomation  Transient  Recorder.   When  enough  data  v/ere 
recorded,  they  were  transferred  to  the  Fabri-Tek  Instrument  Recorder  and 
T^  was  then  calculated  by  th^  PDP  8/E  computer  using  a  FOCAL  program.   A 
series  of  90-t-90  pulse  sequences  were  i.ised,  separated  by  at  least  6T 
in  order  to  allow  the  spin  systev.  to  reach  equilibrium  before  the  next 
sequence.   Thus,  the  first  pulse  would  yield  the  amplitude  at  equilib- 
rium and  destroy  the  z  component  of  the  magnetization,  and  the  second 
pulse  would  yield  the  recovery  amplitude  after  time  t.   The  difference 
between  the  two  amplitudes  decreases  exponentially  with  increasing  t, 
and  the  pair-wise  differences  versus  t  were  used  to   calculate  T  . 

This  method  of  ?neasuring  T,  had  the  advantage  of  requiring  only  a 

1 
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single  sampling  point  on  each  FID.   It  was  not  necessary  to  determine 
the  base  line  which  would  be  required  to  get  the  correct  amplitude,  as 
only  differences  in  amplitudes  were  required.   It  had  the  disadvantage 
of  requiring  a  long  wait  between  pairs  of  pulses  as  a  new  determination 
of  the  equilibrium  magnetization  was  required  each  time.   The  long  wait 
also  subjected  the  data  to  more  errors  due  to  magnetic  field  drift  and 
instabilities  of  the  electronics. 

Another  pulse  sequence,  t-9G-T-180  following  complete  saturation, 
was  also  used.   In  this  method,  a  Princeton  Applied  Research  TDH-9 
Waveform  Eductor  was  used  to  record  the  echo  following  a  180°  pulse. 
The  echo  was  then  displayed  on  an  oscilloscope  where  the  echo  amplitude 
could  be  determined  by  comparison  with  a  scale  on  the  screen.   The  base 
line  on  both  sides  of  the  echo  was  visible  so  the  absolute  amplitude 
could  be  determined. 

This  method  was  much  faster,  as  it  was  not  necessary  to  wait  be- 
tween pulse  sequences.   However,  it  depended  on  visual  estimation  of  the 
pulse  height  rather  than  an  electronically  determined  amplitude  averaged 
over  a  short  section. 

3.3.2   Continuous  Wave  Spectrometer 

A  continuous  wave  (cw)  spectrometer  was  used  to  detenaine  the   N 
NMR  line  shape  of  solid  "  N^  and  of  a  '  U^:  '"R-   "'N:   N  =  l:^:h   mixture. 

A  locally  constructed  Robinson  oscillator  could  be  operated  over 
the  entire  range  of  frequencies  (270  kHz  to  4  Milz)   with  the  appropriate 
rf  coil  and  capacitance  in  r.he  s?jnple  circuit.   It  could  be  frequency 
swept,  the  sweep  being  controlled  by  vriCage  variable  capacitors.   A 
S-kHz  sweep  at  270  kEz  turned  o'l.t  to  be  nonlinear,  and  therefore 
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magrietic  field  sweep  was  normally  used. 

A  Princeton  Applied  Research  Model  HR-8  Tjock-In  Amplifier  was  used 
to  detect  the  signal  and  a  Hewlett  Packard  Model  201C  Audio  Amplifier 
and  a  Ling  Electronics  Model  TPlOO  Power  Amplifier  provided  field  modu- 
lation. 

The  derivative  line  shape  was  recorded  by  the  Fabri-Tek  Instrument 
Computer  which  was  also  used  to  integrate  the  line  shape  to  obtain  the 
absorption  line  shape.   Both  derivative  and  integrated  line  shapes  were 
recorded  on  paper  by  an  X-Y  recorder. 


CHAPTER  IV 
RESULTS  A}TD   DISCUSSION  -  GLYCINE 


4.1   Glycine.  Structure 

3+ 

Glj'cine,  NH„CH  COOK  (or  NH  CII  00  ),  is  the  simplest  of  the  amino 

acids.   It  crystallizes  in  one  of  three  forms,  a,    6,  or  y>  depending  on 
the  preparation  conditions.   Crystals  grox^m  from  an  aqueous  solution 
form  Ci  glycine,  the  most  common  and  stable.   Addition  of  ethyl  alcohol 
to  the  solution  results  in  6  glycine,  and  y  glycine  is  obtained  by  addi- 
tion of  acetic  acid  or  ammonium  hydroxide  to  make  the  solution  acidic 
or  basic. 

Structures  of  the  a,  S,  and  y  forms  were  reported  by  Albrecht  and 
Corey  [l5],  litaka  [l7],  and  litaka  [18],  respectively,  and  Marsh  119] 
reported  a  refinement  of  the  structure  of  a  glycine. 

Crystals  studied  in  the  present  work  were  of  the  a  glycine  form, 
grown,  cut,  and  mounted  at  the  University  of  Nottingham.   Reference  [16] 
has  some  excellent  drawings  of  a  glycine,  though  reference  [19]  probably 
contains  more  accurate  data. 

The  molecule  is  shown  in  Fig.  14.   Two  enantiotuorp.^ious  forms  of  Lh& 
molecule  occur.   The  second  form  would  be  obtained  by  reversing  the  sign 
of  all  the  numbers  in  parent?ieses  in  Fig.  14,  i.e.,  putting  atoms  shown 
above  the  0  ,  0   ,  C   plane  below  and  vice  versa.   Only  one  form  occurs 
in  3  and  y  glycine,  whereas  both  are  found  in  a  gij'cine. 

The  two  forms  described  above  would  be  mirror  imaftes  of  each,  other 
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Figure  14. 


Glycine  molecule.   (a)  Viewed  normal  to  the  0^,  Ott   C 
plane.   liond  distances  are  in  A.   Numbers  in  pirenthesll  are 
distances  in  A  above  (+)  or  below  (-)  the  plane.   (b)  View 
parallel  to  the  Oj,  Ojj,  Cji  plane  showing  the  molecule  is 
rattier  flat  except  for  the  protons. 
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In  a  reflection  plane  parallel  to  the  0^,  0   ,  C   plane.   A  point  in- 
version  of  one  form  would  also  result  in  the  other  form,  in  v;hich  case 
any  line,  such  as  the  lines  representing  bonds  in  Fig.  14,  would  be 
antiparallel  to  the  corresponding  line  in  the  other  form. 

The  ct-glycine  crystal  is  monoclinic  with  four  molecules  per  unit 
cell.   Unit  cell  dimensions  are  [l9] 

a  =  5.1020  ±  0.0008  A 

b  =  11.9709  +  0.0017  A 

c  =  5.4575  t   0.0015  A 

B  =  111°  42.3'  ±  1.0' 

A  typical  crystal  is  shown  in  Fig.  15,  together  with  a  set  of 

crystal  axes  Z  =  X,  Y,  Z,  v/hich  will  be  used  in  the  discussion.   Figure 

15  is  similar  to  Fig.  1  in  reference  [ij  for  purposes  of  consistency. 

>■ 
Figure  16  shows  the  arrangement  of  molecules  in  ci  glycine.   The 

molecules  are  identified  by  A^,  A^ ,  B  ,  B   corresponding  to  D,  C,  A,  3, 

respectively,  in  reference  [16].   The  N  atom  in  B   is  bound  to  0   in 

B^  and  to  0^  in  an  adjacent  B   molecule,  binding  3   and  B^  layers  in  a 

double  layer.   The  double  B  layer  is  loosely  bound  to  the  double  A  layer 

by  van  der  Waals  forces,  which  is  the  reason  for  the  (010)  cleavage 

plane. 

4-2   Electric  Field  Gradient  Tensors 

There  are  4  inequivalent  N  sites;  however,  B,  and  B,  are  related 
by  point  symmetry.   Components  of  the  EFG  tensors  at  the  N  site  in  B, 
and  B^  are  therefore  antiparallel,  ^.nd  being  a  second  rank  tensor, 
parallel  and  antiparallel  cannot   be  distinguished.   Thus,  EFG  tensors 
are  equivalent  at  the  two  sites.   Considering  the  A  molecules,  there 
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Figure  15.  Typical  a-glycine  c; cyr. t:al .  The  {011}  and  {120}  faces  appear 
most  frequently  although  {210}  faces  also  appear.  The  {010} 
faces  are  readily  obtrained  );y  cleavage. 
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Figure  16.   Unit  cell  of  a   glycine.   Axes  are  defined  in  Fig.  15.   Mole- 
cules At  and  Bt  are  of  the  enantiomorphic  form  shown  in  Fig, 
14 J  and  are  related  by  point  inversion  symmetry  to  A2  and  B2. 
Moleculea  B]^  and  B2  form  a  double  layer  which  is  loosely 
bound  to  tha  A-^,  A^  layer. 
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Figure  17.   A  vector  r-,  undergoing  a  180°  rotation  about  the  Y  axis  would 
bring  it  to  r^ ,  whereas  a  reflection  in  the  X,Z  plane  would 
7> 


bring  it  to  r-,,  where  f^  =  -To.   A  tensor  component  has  no 
positive  or  negative  sense,  and  therefore  the  180°  Y-axis 
rocation  and  the  X,Z  plane  reflection  are  equivalent.   Kence 
the  EFG  tensors  from  sites  A  and  B  are  related  by  a  180° 
Y-axis  rotation  or,  equivalently ,  an  X,Z  reflection  plane. 
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Rotation  axis  - 


Figure  13.   Symmetry  of  rotation  patterns  for  rotation  axes  in  or  normal 
to  the  X,Z  plane.   Let  A'  and  B'  be  projections  of  A  and  E, 
respectively,  in  the  rotation  plane  for  any  axis  in  the  X,Z 
plane.   T?ie  projections  are  symmetric  about  the  Y  axis  and 
about  the  X,Z  plane  (some  detail  omitted  below  the  X,Z  plane 
to  rriduce  clutter).   For  a  Y-axis  rotation,  rotation  patterns 
from  site.s  A  and  B  merge  as  projections  in  the  rotation  plane 
are  parallel. 
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are  two  inequivalent  EFG  tensors,  one  for  site  A  and  one  for  site  B. 

Reference  to  Fig.  16   shows  that  a  180°  rotation  of  the  A  layers 
about  an  axis  parallel  to  the  Y  axis  would  make  the  A  layer  equivalent 
to  the  B  layer.   The  two  EFG  tensors  are  therefore  related  by  such  a 
180°  rotation,  or  equivalently  by  a  reflection  in  the  X,  Z  plane,  as 
shown  in  Fig.  17. 

It  is  seen  from  Fig.  18  that  the  rotation  patterns  obtained  from 

sites  A  and  B  are  symmetric  for  any  axis  in  the  X,  Z  plane  when  the  inl- 

tial  orientation  is  chosen  '>7ith  H   in  the  X,  Z  plane  or  normal  to  it. 

o 

Site  A  and  a  clocla^ise  rotation  results  in  the  same  pattern  as  site  B 
and  a  counterclockwise  rotation.   Furthermore,  this  requires  the  two 

patterns  to  coincide  when  H  is  in  the  Diane  or  normal  to  it.   A  rota- 

'^  o 

tion  about  the  Y  axis  would  result  in  a  single  curve,  as  the  projections 

y 

of  both  EFG  tensors  on  the  rotation  plane  v;ould  coincide. 

Tliese  considerations  a.cfect  the  choice  of  rotation  axes.   A  Y-axis 
rotation  would  result  in  a  single  curve  for  (2Av)  ,  taking  no  advantage 
of  the  fact  that  there  are  two  inequivalent  sites,  other  than  alignment 
information.   A  single  curve  would  confirm  Chat  the  rotation  axis  was 
parallel  to  the  Y  axis,  and  two  curves  would  mean  the  axis  was  not  quite 
parallel  to  Y. 

A  rotation  axis  in  the  X,  Z  plane  would  result  in  two  curves,  sym- 
metric about  the  3  =^  0  or  B  =  90°  orientations.   Thus,  the  0-^0  orien- 
tation could  be  accurately  determined,  and  if  the  curves  were  not  sym- 
metric in  amplitude  or  offset,  it  would  be  due  to  the  rotation  axis 
being  scmewbat  out  of  the  X,  Z  plane.   Again,  the  symmetry  of  the  two 
sites  would  provide  only  orientation  information,  and  three  such  rota- 
tion patterns  would  bi:i  needed  to  solve  the  problem  completely. 
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Using  the  convention  that  upper  signs  go  with  site  A  and  lower  with 
site  B  where  the  signs  differ,  X,  Y,  and  Z  rotations  would  result  in  the 
following  rotation  patterns: 


=  A^  +  B^cos29^  t   C^sinZe^ 


(2Av)^  -  ^-^+   (By  +  C^) "^005  2  (9.^  -  6^) 
=  A,^  +  3^00329^  +  C.^3in29Y 


(2Av),  =  A^  +  (bJ  -I-  C?)^i 


^)W2(9^:p5^) 


=  A^  +  B^cos2S^  ±  C^sin29^ 


(4.2.1) 


^ 


where  tan25  =  C/B. 

Any  of  the  six  curves  will  fur»tish  the  diagonal  elements  of  the 
EFG  tensors,  wnich  are  the  same  for  both  sites.   One  off-diagonal  ele- 
ment comes  from  each  rotation,  and  off-diagonal  elements  from  the  two 
sites  have  the  same  magnitude,  differing  in  sign  only  in  two  cases.   The 
complete  tensors,  in  terms  of  coefficients  in  equation  (4.2.1),  are 


1 

K 


(4.2.2) 


-S  ^^x  -^'z 


where  other  choices  could  have  been  made  for  the  diagonal  elements. 
Here  K  =  3eQ/2L,  as  in  (2.4.11)  and  will  be  suppressed  for  simplicity 
in  the  following. 

The  syrnmetry  of  the  EFG  tensors  at  sites  A  and  B  is  displayed  in 


(4.2,2),  where  the   only  differences  in  the  elements  are  (J> 


XY 


and 
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A      B 

(|)*   =  ~<i>„y-      For  a  general  choice  of  crystal  axes,  the  elements  of  the 

L  tit  1.  cj 

A   tensor  would  not  hcve  the  same  magnitude  as  those  of  the  B  tensor, 
unless  by  coincidence. 

Suppose  a  Z  =  U,  V,  W  frame  is  chosen,  and  a  V-axis  rotation  is 
performed.   The  rotation  patterns  for  sites  A  and  B  would  be 
<2A^'),,  =  \  +  B.^cos20^  +  C^sin29^ 

(2Av)^  =  A^  +   ByCos29^  +  C^sin2ey   ,  (4.2.3) 

for  sites  A  and  B  respectively,  the  prime  indicating  site  B.   The  cor- 
responding tensors  would  be 


4>- 


uv 


\v 

\-^\ 

-^w 

^v^^; 

(4.2.4) 


-c 
_v 

-c 

V 


the  upper  and  lower  values  belonging  to  the  site  A  and  site  B  tensors, 
respectively,  as  indicated  by  the  primes.   Here  i|;'s  have  been  used  to 
avoid  confusion  in  distinguishing  between  the  Z  and  Z      frames.   There 
are  four  unknowns,  compared  to  only  two  which  would  result  in  tensor 
(4.2.2)  if  a  single  Y  rotation  were  performed.   If  the  tensor  (4.2.4)  is 
transformed  to  the  Z   frame,  it  miifjt  have  the  symmetry  of  (4.2.2)  which 
permits  solving  for  the  foiir  unknowns. 
The  transformation  is 


*k.i-  > 


-I 


/  a,  .ij.'.  ,a.,  =  /  a  a,  .4>.  , 
.'-.  ki^n  t1  .-.  kx  Ij  xj 
1,3  i,J 


(4.2.5) 


where  k  ,  1  stand  for  X,  Y,  Z,  and  1  ,  j  ^  U,  V,  W,  and  a   is  the  direc- 
tion cosine  of  the  k  (=  X,  Y,  or  Z)  axis  with  respect  to  the  i(--U,  V,  or 
W)  axis.   Note  a.,  -a,.. 
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Using  primes  to  indicate  elements  of  the  B  tensor  in  (4.2.2)  we 
have  <i,^^  t  <i>'^^   =  0,  where  the  upper  sign  holds  for  kl  =  X^   or  YZ,  and 
the  lower  for  kl  =  XX,  YY,  ZZ,  and  XZ.   This  results  in  six  equations: 

^xx  -  *ix  =  J.  ^xi^hj  -  KP   =  °        (^-^-s) 


("yy  "  *YY  "   ^  ^v-t^v-i^'^.  -  *'-.)  =  0  (4.2.7) 


>.   Yi  Yj^^ij    '*'ij' 


^ZZ  -  *ZZ  =  J     ^Zi^Zj^^ij  -  '^li>  =  0  (^-2.8) 


^XZ  -  *^  =  J.    ^Xi^Zj^'^ii  -  ^IP    =  0  (^-2.9) 

*XY  -^  *iY  =  .^.  ^Xi^Yj'^'ij  -"   *lj>  =  °  (4.2.10) 


i.J 


v' 


*YZ  ■*■  *Y7  ■"  I  ^vi-^z^C'^'i-i  +  ■4':.)  =0  (4.2.11) 

Equations  (4.2.10)  and  (4.2.11)  can  be  solved  for  (if),  +  i|)'  )  and 

uv   uv 

^"''VW  "*"  "^W^  which  are  the  only  anknowns  appearing.   Any  two  of  the  first 
four  equations  can  be  solved  for  (4*^^,  -  ^j;^^)  and  (v.^^  -  i-^)  ,    the  only 
unknowns  in  these  equations.   Finally,  (il/^.  +  w^^^  and  (4>        -   •>/,')  can 
be  solved  for  ^p^   and  -4-,^^,  and  likewise  for  i>^^   and  ^^J^^   using  the  other 
two  equations. 

Thus,  a  single  rotation  actually  overspecif i&s  the  solution,  as 
only  four  of  the  six  equations  are  required.   The  only  disadvantage  is 
that  less  information  about  ulignr,.ent  accuracy  is  orovided.   Misalign- 
ment will  cause  the  solution  to  depend  on  which  four  of  the  six  equa- 
tions are  chosen. 

As  an  example,  calling  Anoersscn's  [l]  third  rotation  a  V-axls 
rotation,  the  followj.iig  tersor  results: 
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-322.4 


^UV 


-199.3 


1307.3 

^"uv 

-997.0 

'^UV 

1540.0 

"^VW 

^UV 

-400.0 

-199.3 

-1217.6 

-997.0 


V. 


vw 


-907.3 


(4.2.12) 


(a^.)  = 


Here  the  U,  V,  W  axes  have  been  chosen  such  that  the  direction  cosines 

of  X,  Y,  Z  with  respect  to  U,  V,  W  are 

.7839     -.5313      .3212 

.6208      .6709     -.4056     .        (4.2.13) 

0        .5173      .8558 

The  equations  (4.2.6)  through  (4.2.9)  can  be  used  with  values  from 

(4.2.12)  and  (4.2.13)  to  solve  for  (^jj^^^   -   -^^^j^,)  and  (-^^^  -  '^^^)  .   One  of 

the  first  three  equations  is  a  linear  combination  of  the  other  t\vo,  and 

therefore  any  two  of  the  first  three  give  consistent  results.   The 

values  of  tjj,,,,,  ■4).'„,  i^^,,  and  ij;'  ,  are  listed  in  Table  1  together  with 
UV    uV    VW         VW 

the  equations  used  in  obtaining  the  values. 

Table  1.   Values  of  the  unknoOTi  EFG  tensor  elements  calculated 
using  the  equations  listed  together  with  Eqs.  (4.2.10)  and  (4.2.11). 


Solution 
number 


Equations 
used 


■UV 


"^UV 


w 


'VI* 


I 

(4.1.6,  7, 

8) 

--247 

-281 

-615 

-286 

II 

(4.1.6,  9) 

-258 

-269 

-5S7 

-314 

III 

(4.1.7,  9) 

107 

-634 

-74 

-S27 

IV 

(4.1.8,  9) 

-278 

-249 

-615 

-236 

The  values  in  the  table  were  calculated  by  transforming  frcia  the 
U,  V,  W  frame  to  the  X,  Y,  Z  frame,  using  only  the  syr.unetry  properties 
of  the  tensor  in  the  latter  frame.   The  actual  transformations  can  now 
be  performed  using  Table  1  values  in  tensors (4. 2. 12)  and  (4.2,13). 
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The  following  four  tensors  result,  corresponding  to  the  four  sets 
of  values  in  Table  1: 

-814 


^26  +837  _g^o 
+837  598  ±548 
fill     ±548     -1024 


(4.2.14) 


II 


426 

+827 

-827 

+827 

573 
623 

±554 

-827 

±554 

-999 

-1049 


(4.2.15) 


III 


-53 

906 


-827 


+534 

598 

±859 


-827 

±859 

-545 
-1504  j 


(4.2.16) 


IV 


^52     _„,_ 
400     +^^^ 


q:843 


572 
624 


-827 


:538 


(4.2.17) 


-827     ±538     -1024  ^ 
It  is  noted  that  if  averages  are  taken  where  tlie  A  and  B  site 
elements  differ,  all  the  diagonal  elements  and  the  X,  Z  off-diagonal 
elements  are  consistent  among  all  four  tensors.   Furthermore,  it  seems 
clear  that  (4.2.16)  should  be  disregarded.   The  spread  in  values  iu 
(4.2.16)  arises  from  division  by  the  difference  of  two  nearly  equal 
numbers.   Tensor  (4. 2.]  4')  is  the  most  consistent  and  the  X,  Y  and  Y,  Z 
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elements  fall  between  corresponding  values  in  (4.2.15)  and  (4.2.17) 
Therefore,  the  coir.plete  tensor  would  be 

426     +837     -827 
+837      598     1548  (4.2.18) 

-827     ±543     -1024 
using  only  the  third  rotation  reported  by  Andersson  [l]. 
According  to  (4.2.18)  a  Z-axis  rocation  should  give 


(2Av)   =  512  +  841cos2(9   q:  47.9°) 


(4.2.19) 


whereas  Andersson  [l]  reported 


(2Av)^  =  395  -r-  870cos2(9^  +  39°) 


(4.2.20) 


The  differences  are  probably  due  to  uncertainty  in  orientation  of 
the  axes  which  Andersson  e_t  al.  [l]  suggested  may  have  accounted  for 
the  discrepancies  they  observed. 

Using  two  or  three  rotation  patterns  to  solve  the  proljlem  would 
result  in  even  greater  overspecif ication,  and  make  the  problem  of  de- 
ciding which  equations  to  use  even  more  difficult. 

4.3  Data 


In  the  present  experiment,  two  rotation  patterns  were  obtained  at 
about  5°  C  and  one  pattern  each  at  74°  C,  148°  C,  and  167°  C. 

The  c  axis  was  chosen  for  the  first  series  for  several  reasons. 
First,  cleavage  along  the  (010)  and  (010)  planes  results  in  a.   roughly 
hexagonal  crystal,  facilitating  shaping  and  aligning  the  crystal.   Next, 
Andersson  et  a_l.  [l]  had  performed  a  c-axis  rotation,  so  there  was  a 
reference  to  use  in  initially  obtaining  a  nitrogen  resonance  in  glycine, 
and  a  direct  comparison  could  be  made  with  that  work.   Finally,  it  was 
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not  realized  at  the  time  that  a  single  rotation  would  suffice  if  the 
axis  were  not  in  the  a  ,  c  plane  or  normal  to  it. 

The  first  rotation  pattern.  Fig.  19,  produced  results  deceptively 
close  to  those  reported  by  /mdersson  ejt  a^.  [l].  It  was  thought  that 
the  alignment  was  quite  good  and  that  only  a  part  of  the  rotation  pat- 
tern would  be  required,  making  use  of  symmetry  to  obtain  all  four  com- 
plete curves. 

The  signal  was  extremely  weak,  and  results  shown  in  Fig.  19  were 
the  first. obtained  after  many  fruitless  efforts.   V7ith  all  equipment 
finally  working,  emphasis  was  on  taking  data,  leaving  the  detailed 
analysis  until  later.   Therefore,  it  was  not  recognized  at  the  time  that 
the  curves  were  not  quite  symmetric. 

A  second  rotation  pattern.  Fig.  20,  was  obtained  at  74°  C,  which 
appeared  to  have  a  reduced  offset  and  amplitude. 

The  third  pattern,  Fig.  21,  was  extremely  difficult,  both  electron- 
ically and  mechanically,  as  the  goniometer  would  not  turn  at  this  tem- 
perature.  Some  resonances  could  not  be  found  where  expected  and  some- 
times doubtful  signals  were  recorded  if  they  appeared  near  the  expected 
frequency.   An  effort  was  made  to  define  one  of  the  four  curves  quite 
accurately  and  rely  on  a  fevr  points  and  symmetry  for  the  others. 

The  data  showed  none  of  the  behavior  reported  by  Blinc  et^  al^.  [2] 
with  regard   to  the  shift  in  the  principal  axis  direction.   A  reduction 
in  the  quadrupole  coupling  constant  appeared  to  be  the  major  change. 
Therefore,  one  more  run  was  attempted  at  167°  C  with  understandably 
greater  difficulty.   This  run  continued  until  hot  gas  from  the  gas-flow 
temperature  system  melted  the  insulation  in  the  rf  coaxial  cable,  caus- 
ing it  to  short.   Enough'  data  were  obtained  to  plot  seme  of  the  curves. 
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so  no  more  c-axis  rotation  measurements  were  attempted.   Data  are  shown 
in  Fig.  22. 

After  a  design  change  to  withstand  the  higher  temperature,  a  second 
crystal  was  inserted,  mounted  with  the  rotation  axis  normal  to  the  (120) 
crystal  plane.   The  signal  was  weaker  than  during  the  previous  runs, 
possibly  due  to  the  change  in  design  of  the  equipment  or  to  the  quality 
of  the  crj'Stal.   Use  of  Fast  Fourier  Transform  (FFT)  techniques  aided, 
but  the  only  rotation  pattern  obtained  was  at  5°  C,  shown  in  Fig.  23. 
Each  suspected  data  point  had  to  be  repeated  several  times  for  verifi- 
cation, and  small  angular  rotations  were  used  initially  until  the  trend 
of  the  curves  was  established. 

The  signal  could  not  be  distinguished  at  even  moderately  higher 
temperatures,  and  the  experiment  was  terminated.   The  c-axis  rotation 
patterns  demonstrated  normal  temperature  dependence  and  did  not  confirm 
the  change  in  the  orientation  of  the  PA  frame  reported  by  Blinc  et^   al. 
[2]. 

Figures  19  through  23  show  311  the  data  together  with  smooth  curves 
calculated  by  a  least  squares  fit,  except  for  some  of  the  curves  in 
Figs,  21  and  22.   Curves  (3)  and  (4)  in  those  two  figures  were  deter- 
mined quite  accurately,  but  the  few  data  points  thought  to  lie  on  curves 
(1)  and  (2)  turned  out  to  be  unreliable.   No  consistent  set  of  curves 
could  be  found  utilizing  all  data  points,  and  use  of  all  data  caused  a 
very  poor  fit  to  curves  (3)  and  (4).   Therefore  only  the  latter  two 
cur\'es  were  considered.   Cur\Tes  (1)  and  (2)  were  estimated  from  curves 
(3)  and  (4)  and  the  known  syrametiy. 
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A.4   Data  Analysis 

The  data  were  first  put  in  the  form  used  by  Vclkoff  e^  £l.  [9], 
which  required  a  96°  phase  change  in  the  rotational  orientation  as  well 
as  combining  pairs  of  curves  co  give  the  total  splitting  rather  than 
the  shift  from  the  central  unperturbed  Zeeman  frequency.   The  c-axis 
rotations  became 

T  =  5°  C:   (2Aj)^  =  428  +  917cof2(6   -  38.8°) 


(2Av)^  =  407  +  855cos2(0   f  38.8°) 


T  =  74°  C:   (2Av)t  ■-■   ^24  +  882cos2(0„  -  33.8°) 
(2Av)^  =  408  +  812cos2(e  +  38.8°) 


T  =  148°  C:   (2Av)J  =  325  +  734ccs2(92  -  38.3°) 
(2Av)2  =  312  +  676cos2(e2  +  38.8°) 


T  =  167°  C:   (2Av)J  -  261  +  6140022(6^  -  38.3°) 

(2Av)^  =  244  +  558cos2(9  +  38.8°)   .      (4.4.1) 

Equations  (4.4.1)  should  have  the  same  constants  for  the  A  and  B 
sites  if  the  axis  were  exactly  in  the  a  ,  c  plane,  so  an  average  is  taken 
for  each  temperature.   The  equations  become 

T  =  5°  C:   (2Av)„  =  417.5  +  886.0cos2(e  +  38.8°) 


T  =  74°  C:   (2Av)„  =  416.0  -f-  847.0co32(9  +  38.8°) 
z  z 


T  •-=  ]4S°  C:   (2Av)^  =  318.5  +  705.0ccs2(9^  +  38.8°) 
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Figure  24.   Relationship  of  U,  V,  W  axes  to  X,  Y,  Z  axes.   The  U  axis  is 
normal  to  the  (120)  plane  and  a  =  33.378°. 
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T  =  167°  C:   (2Av)„  =  252.5  +  586.0cos2(9  +  38.8*=*)     (4.4.2) 

where  the  upper  (lower)  sign  pertains  to  site  A  (3). 

A  Z  frame  is  used  for  the  rotation  in  the  (120)  plane,  defined  by 
U  being  norxial  to  the  plane  and  W  being  antiparallel  to  the  Z  axis  as 
sho^^m  in  Fig.  24.   The  curves  obtained  are 

T  =  5°  C:   (2Av)y  =  185.2  +  1444.5cos2(e  +  21.89°) 

(2Av)y  =  -646.7  +  258.10032(6^  +  23.65°)  .   (4.4.3) 

The  two  curves  should  cross  at  9„  =  90°  as  H  would  be  in  the  a  .  c 

U  o 

crystal  plane  and  the  two  EFG  tensors  are  mirror  inages  of  each  other 
in  that  plane.   The  curves  cross  at  0  =  91.25°,  indicating  a  1.25° 
error  in  rotational  orientation.   The  equations  (4.4.3)  then  become 

T  =  5°  C:   (2Av),^  =  185.2  +  1444. 5cos2(0„  +  23.14°) 

(2Av)^  =  -646.7  +  258.1cos2(ajj  +  24.91°)  .   (4.4.4) 

Equations  (4.4.2)  and  (4.4.4)  are  expanded  to  give 

T  =  5°  C:   (2Av)f;  =  185,2  +  998.3cos29„  -  1044.0sin2e„  (4.4.5) 
ij  u  u 


(2^v)J  =  -645.7 -!i66.5cos2Sj,  -  197.2sin29    (4.4.6) 
(2Av)   =  417.5  -*-  190.3cos2e„  ±  865.3sin29_   (4.4.7) 


T  =  74°  C:  (2Au)^  =  416.0  +  lSi.9ccs29^  t   827.23in29    (4.4.8) 


=  148°  C:  (2u^>)   =  316.5  -t-  i51.4co32e^  ±  S83.6sin20,   (4.4.9) 


T=167°  C:  (2Av)^  =  252.5  +  l.o.Scos2G^  t   572.3sin20^  .  (4.4.10) 
Equations  (4.4.5)  and  (4.4.6)  -an  be  'j.^.ed  alone,  or  in  conjunction 
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with  (4.4.7)  to  solve  the  complete  EFG  tensor  at  T  =■•  5°  C.   These 
equations  lead  to  the  following  tensors  in  the  Z  and  E   frames,  where 
again  the  upper  value  or  sign  goes  with  site  A  and  the  lower  value  or 
sign  with  site  B  if  they  differ: 


^*ki^i: 


607.8 
+865.3 

-370.4 
1293.4 


+865.3 

in  .1 


-835-0 


(4.4.11) 


-C. 


W 


1133.5 
-480.2 


1044.0 
197.2 


-r' 

1044. 
197. 

0 
2 

-813.1 


-813.2 


U 


(4.^.12) 


./' 


In  (4.4.11)  use  has  been  madn  of  the  relations  C  =  C'  and  C  =  'C', 
the  primes  again  indicating  the  B  sits. 

The  direction  cosines  of  Z  =  X,  Y,  Z  with  respect  to  Z  =  U^  V,  W 


are 


(^ui>- 


cosa 

since 

0 


sxna 
-cosa 
0 


0 
0 

-1  j 


-  ^^'^^ 


(4.4.13) 


where  a  =  38.378°  as  in  Pig.  24.   The  tensors  (4.4.11)  and  (4.4.12)  are 
related  by 


'^kl 


i,J 


-1 

'ki'ij'^jl 


<f'i  1  ~  Z  ^u -'^- •^.•■.  =  ^  a,  _. ^-i  ,_'l'j  .    "        (4.4.14) 


.^.  ki  Ik  ij 


Again  recalling  the  EFG  tensors  are  synnietric,  equations  (4.4,14) 
lead  to  twelve  equations-,  six  for  each  site.   The  set  of  six  rt  present- 
ing site  A  is 
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?  2 

*XX  "   '^UU'^°^   ""  "^  2-;^j^sinacosa  +  'l^^sin  a  (4.4.15) 

?  2 

{J>        =   t|j      sin  a  -   2:!;      sinacosa  +  'I'yyCOS  a  (4.4.16) 


hz  =  "^WW 


(4.4.17) 


<J>        =   ^  ,r  =  -'I'TrTTCOsa  -   >i;,_-sina  (4.4.18) 

^XZ  ZX  UW  \r.'I 

'f'xY  "  '^YX  "  ^'^UU  ~  *VV^^'"'"'°^""^'''UV'^^''  ""  ~   ""-"^  ^"^        (4.4.19) 
*YZ  =  *ZY  =  -"^Ul.^^""  ^   '"w*^"^"    •  (4.4.20) 

Adding  primes  to  the  <i>' s   and  ijj's  would  result  in  the  equations  for 

the  B  site. 

The  solution  may  be  obtained  from  the  single  U  rotation  by  making- 

use  of  the  symmetry  of  (4.4.11).   Six  equations  result: 

(4^^.21) 
(4.4.22) 

<^zz-*zz  =  *ww-^:m  =  °  ^'-'-^'^ 

hz  -   *iz  =  -^V.  -  '^'uw^'^"^^  -  ^*w  -  '^w>^^""  =  °        ^'-'-'^^ 
^XY-^^^Y  ~-   (^mj-'^'iu- V-*w^^^^"''°^^''^"*uv-' V^"^"^"-"^^"^  ^  "^ 

(4.4.25) 

*YZ  ^  *;Z  =  -^*UW  ^  W^^^'^  ■"  ^"^VW  +  ^W^^^^^  =  °     '         (^•^•^•^> 

The  small  correction  for  rotational  orientation  made  la   equation 
(4.4.4)  assured  equation  (4.4.23)  would  be  very  nearly  true,  which  in 
tura  caused  (4.4.21)  and. (4.4.22)  to  be  consistent.   Because  of  the 
choice  of  axes,  ( "-l^nr  ~  '^™»)  appears  only  in  (4.4.21)  and  (4.4.22)  and 
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either  equation  can  be  used  to  evaluate  it  without  ambiguity.   Also  due 

to  the  choice  of  axes,  ('J'jyT,  ~  '!'['„■)  appears  only  in  (4.  A.  24),  resulting 

in  a  single  value.   The  quantities  (.i,,..  -^  \p'.,)    and  (iii,„,  +  i>'  ,)  appear 

Uv     Uv  UW     UK 

in  (4.4.25)  and  (4.4.26)  respectively,  making  their  evaluation  simple. 
Even  with  a  different  choice  of  axes,  the  latter  two  quantities  would  be 
unambiguous. 

Inserting  values  from  (4.4.12)  and  recalling  a  =  38.378°,  the  un- 
knowns in  (4.4.12)  are 


W 


-^  =  -^29.1 


C*  =  -,h'   =  -37.5 


S  =  -^uw  =  -''^-^ 


^V  =  "^uw 


and  (4.4.12)  becomes 


-1118.9 


(4.4.27) 


(4'..)  ,, 


-370.4 
1293.4 

429.1 
37.5 

448.3 
1118.9 


429.1 
37.5 

1183.5 
-480.2 

1044.0 
197.2 


448.3 
1118.9 

1044.0 
197.2 

-813.1 
-8i3.2 


U  rFrom  U 


(4.4.28) 


J 


'rotatioiv 

The  transformation  to  the  Z  frame  (4,4.14)  can  now  be  performed, 
resulting  in 

+854.6     -999.6  ' 

156.9     ±540.1  ,  (4.4.29) 

±540       -813.2     "~  ~  " 


(\l^Z 


646.3 
+854.6 
-999.6 


„  cFrora  U  -i 
Z !      . 
^rotation'' 


V7hich  should  be  compared  with 
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<\l>Z  = 


607.8     .11365.3 


+365.3 


-S 


227. 


+c. 


-835.0 


(4.4.11) 


(■From  Z  -v 
'rotation'' 


The  differences  are  probably  due  to  small  errors  in  crystal  orien- 
tation. 

In  order  to  solve  the  problem  using  both  rotations  rather  than 

known  symmetry,  values  of  4.^^,  '}>v-v»  "^77  ^^'^   '^w  f'^"™  (A. 4. 11)  are  put 

AJv     11     Zt/t  Xl 

in  equations  (4.4.15)  through  (4.4.20),  together  with  known  values  of 
'4.'s  from  (4.4.12).   Equations  (4.4.15),  (4.4.16),  and  (4.4.19)  each  con- 
tain only  one  unknown,  t^   t  \p'      for  site  A  or  B,  respectively,  which 
are  found  to  be 


\p       =  389.6,  475.6,  357.2 


^p^   =  -2.0.  -9.1,  -24.4 


(4.4.30) 


depending  on  which  of  the  three  equations  is  used.   There  is  now  con- 
siderable variation,  due  again  to  orientation  errors,  but  the  averages. 


il)   =  410.8,    V      =  -11.8 


(4.4.31) 


agree  within  reason  with  the  values  429.1  and  37.5,  respectively,  found 
from  the  single  rotation. 

The  four  equations  (4.4.18),  (4.4.20)  and  the  corresponding  equa- 
tions for  site  B  have  six  unknowns,  and  cannot  be  solved.   Therefore, 
two  more  relations,  (f   =  (jj*   and  *   =  -(^'   must  be  used.   This  has 
already  been  done  in  solving  the  problem  from  the  single  U  rotation,  and 
the  values  can  be  taken  from  (4.4.28)  and  (4.4.29). 

The  resulting  two  tensors  are 
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i'i)     ) 


-370.4 
1293.4 

410.8 
-11.8 

448.3 
1118.9 

410.8 
-11.8 

1183.5 
-480.2 

1044.0 
197.2 

448.3 
1118.9 

1044.0 
197.2 

-813.1 
-813.2 

(4.4.32) 


U  |-From  bothN 
^rotations-' 


and 


(*kl^Z  = 


607.8 
qFS65.3 
-999.6 


+865.3 

227.2 

±540.1 


-999.6 


^540.1 


(4.4.33) 


From  bot;h\ 


835.0  ,  E  r"'°  °°''-") 
^         rotations-' 


In  order  to  illustrate  the  effect  of  small  alignment  errors,  the 
tensors  (4.4.29)  may  be  transformed  to  a  set  of  axes  Z"  =  X",  Y" ,  Z" 
such  that  X"  and  Z''  are  out  of  the  a  ,  c  plane  0.5°  and  therefore  Y"  is 
0.707°  from  the  b  axis.   The  resulting  tensors  are 


631 
661 

-867 
841 

-987 
-1012 

-867 
S41 

191 
143 

540 
-540 

-987 
-1019 

540 

-823 

(4.4.34) 


The  order  of  niagnitudes  of  the  shifts  caused  by  such  a  small  orien- 
tation error  in  comparison  with  the  differences  between  (4.4.29)  and 
(4.4.33)  indic=.tGS  a  cumulatlvi  alignment  error  on  the  order  of  1°  in 
obtaining  the  desired  rotation  axes. 

Following  the  method  of  Vjlkoif  ^  _al.  [9],  the  tensors  (4.4.29) 

2 
and  (4.4.33"^  may  be  diagonalized  to  obtain  e  qQ/h,  n ,  and  the  direction 

cosines  of  the  principal  axes.   These  results  are  displayed  followed  by 

those  of  Andersscn  et^  £il.  [l]  and  Biinc  e_t  al.  [2]  for  comparison. 

Andersson's  group  worked  ;o  t  a  few  degrees  higher  temperature  than  the 

present  v/orlc  but  the  temperature  dependence  is  very  small  over  this 
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<.y)     = 


(4.4.35) 


interval,  and  may  be  ignored. 

Blinc's  group  worked  at  a  significantly  higher  temperature,  but 
their  results  are  included  for  comparison. 

The  tensors  (4.4.29)  from  the  single  U  rotation  result   in 

n  =  .492,    e^qQ/h  =  1188  kHz  (4.4.35) 

'   .513     ±.844     -.158  ' 
.420     +.086      .904 
+.749     -.529     -T.399  _ 

at  T  =  5°  C. 

The  tensors  (4.4.33)  from  the  Z  rotation  with  two  elements  obtained 

from  the  U  rotation  results  in 

n  =  .517,   e%Q/h  =  1189  kHz 

.523     t.834     -.179 
.427     -.074 
+.738     -.547 


(a^j)  . 


.901 
+.395  J 


(4.4.37) 


(4.4.38) 


at  T  =  5°  C. 

Andersson  et^  £l.  [ij  reported 


n  =■ 


.5  ±   .1,        e^Q/h  -   1200  ±  100  VJKz 


(4.4.39) 


(=y)  - 


.53     t.B5  -.17 

.43     T-0?>      -^^  (4.4.40) 

^  ±.74     -.54     +.40  , 

at  room  temperature. 

Blinc's  publishec  EFG  r.crsor  -an  be  used  tc  find  the  direction 
cosines  of  the  principal  axe  J  J.n  Chc:  X,  Y,  Z   frame,  and  his  results 

are  [2] 

n  =  0.61  +  0.03,   e\o/h=  745  +  20  kHz       (4.4.41) 
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.216 

3.391 

-.895 

(a,j)  = 

.902 

+  .430 

.030 

.  +^373 

-.814 

+  .446 

(4.4.42) 


at  T  =  140°  C. 


Edmonds  and  Speight  [3]  reported 


(4.4.43) 


n  =   0.51,    e'qQ/h  =  1250  kH^ 
at  T  -  77  K. 

At  higher  temperatures,  only  incomplete  EFG  tensors  could  be  ob- 
tained, calculated  from  equations  (4.4.8),  (4.4.9),  and  (4.4.10).   They 


are 


T  =  74"  C; 


^\lh   = 


598 
+327 


+827 
234 


+C 


X 


-832 


(4.4.44) 


T  =  148°  C: 


'^\lh. 


470     +689 
+f^.89      167 


-C, 


TC 


X 


+^X 
-637 


(4.4.45) 


T  =  167-"  C; 


378 

+572 

-s' 

kl-^E  ■" 

:o/2 

127 

+^X 

-s- 

^S 

-505 

(4.4.46) 


If  a  constant  asymnittry  parameter  and  constant  direction  cosines 

2 

are  assumed,  but  e  qO/h  decreases  with  temperature  to  1170,  925,  and 

750  kHz  at  T  =  74°  C,  145'-''  C,  and  167°  C,  respectively,  the  tensor 
(4.4.33)  would  become 
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74°  C: 

598 

1-851 

-984 

%l\   = 

-F85]. 

224 

±532 

-984 

±532 

-822 

148'^  C: 

473 

+673 

-778  ' 

^\l>Z  = 

+673 

177 

±420 

-778 

±420 

-650 

(4.4.47) 


(4.4.48) 


T  =   167'^   C: 


^^'.l\ 


384  +546  -631 

T-546  143  ±341 


-631 


:34i 


-527 


Finally,  Blinc's  EFG  tensor  in  the  X,  Y,  Z,  frame  is  [2] 


T  =  140°  C: 


'  -587 

±670 

203 

«ki>.:  ' 

±670 

540 

±492 

203 

±492 

47  j 

(4.4.49) 


(4,4.50) 


There  is  quite  good  agreeraerit  between  the  present  results,  (4.4.35) 
through  (4.4.38),  and  those  reported  by  Andersson  e^t  aA.  [l],  (4.4.39) 
and  (4.4.40),  and  by  Edmonds  and  Speight  [3],  (4.4.43),  even  though  the 
temperatures  are  different.   This  indicates  a  very  weak  temperature 
dependence  over  the  interval  froin  77  K  to  room  temperature,  and  confirms 
the  orientation  of  the  EFG  principal  axes  to  within  about  1°,  which  is 
the  order  of  the  variation  between  (4.4.3G),  (4.4.38)  and  (4.4.40).   The 
direction  of  the  largest  principal  axis  is  about  7°  from  the  direction 
of  the  C-N  bond. 

A  comparison  of  the  obsei/ved  data  at  T  =  74°,  l48°,  and  167°  C, 
(4.4.44),  (4.4.45),  and  (4.4.46),  and  the  results  predicted  by  a  reduced 


e  qQ/h,  (4.4.4?),  (4.4.43),  and  (4. 4.-^3),  indicates  the  major  effect  of 
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increasing  the  temperature  is  a  reduction  of  the  quadrupole  coupling 
constant. 

The  orientation  of  the  principal  axes  observed  by  Blinc  (4.4.42), 
is  completely  different  than  the  orientation  reported  by  ^Jidersson 
(4.4.20)  and  that  found  in  the  present  work,  (4.4.36),  (4.4.38).   The 
latter  orientation  supports  a  hindered  rotation  of  the  M  group,  and 
there  should  be  no  reason  for  it  to  cease  or  become  less  symmetric  at 

higher  temperatures. 

2 
While  a  simple  reduction  in  e  qQ/h  at  higher  temperatures  does  not 

completely  explain  the  EFG  tensors  found  at  T  =  148°  C  and  T  =  167°  C, 

there  is  no  relation  between  these  tensors  and  that  reported  by  Blinc 

(4.4.20). 

4.5   Conclusions 

We  conclude  that  the  principal  axis  of  the  EFG  tensor  lies  about 
7°  from  the  direction  of  the  C-N  bond  at  room  temperature  and  below,  and 
that  the  major  effect  of  increasing  temperature  is  a  decrease  in  the 
quadrupole  coupling  constant  of  about  37%. 

The  reason  for  the  disagreement  with  Blinc  et_   aj^.  [2]   is  not 
understood.   One  possibility  is  that  a  different  form  of  the  crystal  wa.s 
inadvertently  used  by  Blinc  et^  aJ^.  [2].   Another  is  that  the  high  field 
approximation  is  not  adequate  in  the  4,130  G  or  4,280  G  fields  used  by 
Blinc.   The  unperturbed  Zeeman  frequencies,  1.286  and  1.316  l-'iEZy    are  of 
the  same  order  as  the  maximum  quadrupolar  splitting  observed  of  over 
1.0  ^!H7. .   Also  a  simple  error  in  orienting  the  crystal  could  take  place. 


CHAPTER  V 
RESULTS  AND  DISCUSSION  -  NMR  LINE  SHAPE  OF  ^^N  IN 
a  NITROGEN  WITH  O2  IMPURITY 


5.1   Structure  of  a  Nitrogen 

The  structure  of  a  nitrogen,  stable  below  35.6  K,  is  cubic  and 
contains  a  basis  of  four  molecules  per  unit  cell.   A  unit  cell  is  shown 
in  Fig.  25.   The  molecules  are  oriented  parallel  to  the  body  diagonals 
and  may  or  v.ay  not  be  displaced  slightly  from  the  centrosymmetric  posi- 
tions along  the  body  diagonals,  a  question  which  has  not  been  completely 

resolved  to  date.  y 

o 

The  intramolecular  nuclear  separation  is  about  1.1  A  whereas  the 

intennolecular  nearest  neighbor  distance  is  3.57  A  assuming  the  molecules 

o 
are  at  the  centrosymmetric  positions,  or  3.42  A  if  they  are  displaced 

0.17  A  L20J.   Other  displacements  would  result  in  different  distances; 

however  the  r   factor  in  Eq.  (2.2.3)  results  in  the  intramolecular 

local  magnetic  dipolar  field  being  about  30  times  stronger  than  the 

intermoiecular  j.ocal  field. 

The  molecules  undergo  libraticnal  motion  which  causes  a  reduction 

in  the  amplitude  of  a  in  Eq .  (2.2.14)  by  a  factor  of  about  0.86  due  to 

motional  averaging  [2l]. 

The  librat.ional  notion  results  in  a  quadrupolar  dominated  relaxation 

14    -   -, 
time  of  about  100  seconds  at  4.2  pC  in  '  N   1.22.1,  and  a  much  longer  relax- 

15  15 

ation  time  is  expected  in   N   as  the  spin  %  '^  N  nuclei  have  no  quadru- 
polar interaction.   There.'rore,  in  order  to  conventi^ntly  observe  the 
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Figure  25.   Unit  cell  of  a  N-.   The  four  basis  molecules  ara  at  the 
comers  of  the  small  nested  cube,  each  oriented  parallel 
one  of  the  four  body  diagonals.   Distances  are  in  A. 


to 
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spectnim  in   N,,  paramagnetic  0  was  added  to  shorten  the  relaxation 

time  during  some  runs. 


5.2  Data 

15  1^ 

Using  a  99.5%   N  sample, with  the  remainder  being   N,  about  0.1% 

0  was  added  initially. 

At  th<i  lowest  field  of  625  G.  the  spectrum  vary  closely  resembled  a 
classical  Pake   powder  pattern  doublet  with  some  asymmetry.   The  asym- 
metry increased  markedly  with  increasing  H   up  to  the  limit  of  the 
magnet,  9.0  kG.   The  data  are  shown  in  Figs.  25  through  29. 

The  data  at  625  G  and  950  G  are  shown  in  Fig.  26.   The  derivative 
line  shape  as  recorded  by  the  CW  spectrometer  is  sho\«i  as  a  dashed  line 
in  Fig.  26(a)  with  the  absorption  signal  obtained  by  Integration  super- 
imposed as  a  solid  line.   Figure  25(b)  shows  the  data  as  a  solid  line 
superimposed  on  broadened  and  unbroadened  theoretical  Pake  doublets  frora 
Fig.  1,  which  are  sho'.%m  by  dotted  and  dashed  lines  respectively. 

Examination  of  Fig.  26(b)  can  leave  little  doubt  that  tl^e  spectrum 
Is  anything  but  the  expected  Pake  doublet,  although  the  high  frequency 
wing  appears  to  be  somewhat  too  close  in. 

Figures  27  and  23  shovv'  the  effect  of  increasing  the  magnetic,  field 

R  .   The  derivative  spectrum  is  again  superimposed  on  the  absorption 
o 

spectrum  at  9000  G  which  shows  the  low  frequency  wing,  although  low  in 

am.plitude,  has  a  well-defined  beginning. 

The  high  frequency  wing,  which  was  closer  in  than  predicted  by  the 

Pake  doublet  theory  at  625  G,  moves  even  closer  as  H   i.^  increased  until 

-^  o 

it  cannot  be  distinguished  above  2  kG. 

Other  structure  is  nearly  masked  by  the  basic  doublet  spectrum. 
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N2  with  0.1%  O2  added  at  625  and 


Figure  26.   NMR  spectrum  of  '"N  in  a 

950  G.   (a)  Derivative  spectrum  shovni  by  dashed  line  at 

625  G.   (b)  Theoretical  unbroadened  (  )  and  broadened 

(  •••  )  Pake  doublet  superimposed  on  data  at  625  G.   (c)  Data 
at  950  G  showing  lov  frequency  (high  field)  wing  extending. 


Curves  were  obtained  by  sweeping  Hq,  where  h  -  Hq 


H  = 


H^  -  w/y  from  Eq.  (2.2.14). 
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Figure  27, 


tMR  spectrum  of   N  in  ct   N2  with  0.1%  O2  added  at  2000, 
3000,  and  4000  G. 


-115- 
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Figure  28.   N'MR  spectrum  of   N  in  a   "No  with  0.1%  0^  added  at  5000, 

7000,  and  9000  G.   (a)  5000  G.   (b)  7000  G.   (c)  Derivative 
spectrum  shown  by  dashed  line  at  9000  G  showing  low  fre- 
quency wing  has  a  well-defined  beginning. 


7000  G 


c.   0.5%  0, 
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b.   0.1%  0, 


a.   <  0.02%  0, 


2     h(G) 


Figure  29.   MMR  spectrum  of   N  in  a   N2  as  a  function  of  ©2 
concentration. 
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Thus  a  small  hump  appears  between  the  two  main  peaks  at  low  fields,  but 

is  no  longer  distinguishable  above  5  kG.   It  has  either  lost  relative 

intensity  or  is  combined  with  the  low  frequency  peak  at  high  H  . 

o 

The  peak-to-peak  separation  remains  nearly  constant,  although  it  is 
a  little  less  in  medium  fields  than  in  the  highest  and  lowest.   This 
effect  is  probably  due  to  the  fact  that  the  spectrum  is  a  combination  of 
a  Pake  doublet  and  a  Less  intense  spectrum.   When  the  two  peaks  are 
added,  the  resultant  peak  lies  between  the  two,  giving  the  appearance  of 
shifting  the  more  intense  peak. 

The  Pake  doublet  should  be  centered  about  the  unperturbed  Zeeman 

frequency  \>    ;  hovjever  v   lies  between  the  center  and  the  lov;  frequency 
o  o 

peak,  meaning  the  entire  Pake  doublet  was  shifted  to  higher  frequencies 
by  about  0.5  G  or  0.25  kHz.   This  shift  is  not  displayed  in  the  data 
where  the  zero  point  \7as  taken  to  be  midway  between  the  two  peaks. 

Following  the  experiments  on  the  0.1%  0„  sample  the  gas  was  passed 
through  an  Oxisorb  [23]  cartridge  several  times  to  remove  the  0„.   As  a 
check  on  the  success  of  removal  the  spin-lattice  relaxation  time  T^  was 
measured  at  77.4  K  and  compared  with  T^  of  high  purity   N^  [24].   These 
measurements  placed  an  upper  limit  of  0.02%  0„^  remaining,  the  limit 
being  set  by  the  accuracy  of  the  data.   The  Oxisorb  specifications  guar- 
antee less  than  0.1  ppm  0^  remaining  for  certain  pressures  and  flow 
rates;  however  not  all  of  the  sample  could  be  withdrawn  from  the  outlet 
side  of  the  Oxisorb  unit  each  time.   The  sample  was  cycled  repeatedly 
through  the  Oxisorb  until  it  was  calculated  that  less  than  1  ppm  0^  re- 
mained. 

Figure  29(a)  is  tae  spectiura  at  7000  G  and  T  =  4.2  K  obtained  with 

the  0„  removed.   At  4.2  K,  T.  appeared  to  be  on  the  order  of  seven  hours, 
2  J.    ' 
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makins  it  very  difficult  to  obtain  a  good  spectrum  in  the  presence  of 
saturation.   By  leaving  the  magnet  on  continuously,  one  or  tv;o  traces 
could  be  obtained  each  day,  and  all  were  reproducible  to  the  extent  that 
the  two  peaks  were  narrower  than  the  corresponding  data  with  0.1%  0^  and 
the  relative  intensities  were  unchanged  by  the  removal  of  the  O^.   No 
other  features  could  be  discerned. 

Figure  29(c)  shows  the  effect  of  adding  0.5%  O2  where  the  broaden- 
ing has  iTiasked  any  structure. 

5_.  3  Data  Analysis 

In  searching  for  an  explanation  of  these  observations,  we  look  to 
the  effects  of  the  constituents  of  the  sample  other  than   N  (O2  and 
"'*N)  and  to  the   N  itself.   The  literature  contains  samples  of  spectra 
showing  effects  due  to  paramagnetic  ions  or  0^  which  have  a  deceptive 
similarity  to  the  present  data  [25],  [26].   References  [27]  and  [28] 
examine  the  case  of  ortho-para  conversion  in  spin  %  solid  H^,  and  again, 

the  resulting  spectra  are  deceptively  similar  to  the  present  data.   Both 

14 

of  the  above  nfficti  are  analyzed  below  as  is  the  effect  of  the   N 

atoms  which  is  also  examined  theoretically  and  experimentally  in  Chap.  6. 


5.3.1  Paramagret.ic  0   Impurity 


The  most  obvious  cause  of  t.ha  asymmetry  and  other  anomalous  features 
of  the  line  shape  observed  is  the  paramagnetic  0- .   Piott  and  McCormick 
[25]  have  examined  theoretically  and  experimentally  the  effect  of  O2 
impurity  on  the  N>rR  line  shape  of  protons  in  solid  methane.   They  observed 
a  central  line  and  anomalous  peaks  in  the  wings,  one  of  vjhich  is  dominant 
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Broadening  of  the  central  line  was  attributed  to  protons  which  are 
second  and  more  distant  neighbors  of  0   molecules,  whereas  the  peaks  in 
the  wings  are  due  to  first  nearest  neighbors. 

The  large  central  peak  and  the  anomalous  peak  in  the  wing  formed  an 
asymmetric  pair  which  was  first  compared  with  the  asymmetric  pair  ob- 
served in  the  present  data  in  the  4  to  7  kG  range.   This  identification 
of  the  peaks  had  to  be  ruled  out  when  measurements  at  lower  fields 
(625  G)  proved  the  peaks  to  be  the  expected  Pake  doublet  peaks. 

Piott  and  McCormick  give 

h  =  g3r~^<S„>(l  -  3cos^0)  (5.3.1) 

as  the  local  field  at  a  proton  site  due  to  the  electron  spins  in  a 
neighboring  0  molecule.   lliey  use 

H   =  gpH  S„  +  D'[sJ,  -  (1/3)S(S  +  1)1  (5^,3.2) 

e      o  Z       Z 

as  the  Hamiltonian  to  calculate  <S  >.   There  are  12  nearest  neighbors 

and  steric  considerations  limit  *^'S  >  to  chree  values,  resulting  in  36 

values  of  h. 

Using  a  computer,  they  calculated  the  derivative  1  ine  shape  for  100 
randomly  oriented  crystallites  due  to  the  nearest  neighbors  of  0^  mole- 
cules in  the  field  produced  by  the  C„  molecules.   A  simple  integraclon 
of  their  derivative  line  shape  is  shcwi  in  Fig.  30  as  a  solid  line.   The 
dashed  line  qualitatively  represents;  the  location  of  the  central  peak  as 
it  cannot  be  accurately  calculated  from  the  data  in  ref.  [25]. 

The  identification  of  the  two  peaks  in  ref.  [25]  with  the  two  in 
the  present  data  has  already  been  ru.ied  out.   There  is  a  possibility 
that  the  present  data  area  superposition  of  a  spectrum  due  to  the  nearest 
neighbors  of  0^  molecules  (hereafter  referred  to  as  the  paramaj^natic 
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spectrum  or  line  shape)  and  the  Pake  doublet.   Two  such  spectra  are 
qualitatively  indicated  in  Fig.  31  at  H   =  950  G  and  H   =  7000  G. 

^  °  O  Q 

In  order  to  predict  the  behavior  of  the  paramagnetic  line  shape  it 

is  easier  to  folloi^.'  the  work  of  Kroon  [25]  wlio  examined  the  proton 

spectrum  in  NH,Fe(SO,),  experimentally  and  theoretically  and  obtained 
4     4  2 

+++ 
the  proton  spectra  shown  in  Fig.  32.   The  paramagnetic  Fe    ion  pro- 
vides the  local  field  rather  than  0   impurity  and  there  is  no  central 
line  as  all  of  the  protons  are  nearest  neighbors  of  an  Fe    ion.   Kroon 
arrived  at  an  analytic  expression  for  the  line  shape  Vi7hich  is  more  con- 
venient to  use  than  the  random  crystallite  calculations. 

Considering  that  the  magnetic  moment  of  an  electron  is  on  the  order 
of  1000  times  that  of  a  nucleus,  it  might  be  expected  that  the  local 
field  produced  by  an  electron  at  a  nearby  nucleus  would  be  as  much  as  a 
icilogauss,  which  would  spread  the  powder  pattern  spectrum  over  such  a 
region  as  to  make  it  unobservable.   However,  if  the  electron  occupies 
all  of  its  states  over  a  time  short  with  respect  to  T^  and  T  of  the 
protons,  an  average  iTicment  must  be  used.   The  classical  average  is  given 
by 

^D  =  ^'^D^\/^''^  (5.3.3) 

and  has  a  fixed  direction  p.^rallel  to  the  external  field.   Here  p   is 
the  large  magnetic  moment  of  the  paramagnetic  dipole.   The  fact  that  the 
paramagnetic  spectrum  is  observable  is  experimental  evidence  that  the 
electron  is  changing  states  rapidly  enough  to  make  Eq.  (5.3.3)  valid. 

The  parallel  component  of  the  magnetization  due  to  the  ion  is 

H,,  -  Uj^r"^(3cos-'e  -  1)  (5. 3. A) 

and  the  resulting  unbroadr-ned  line  shape  is  given  by 
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Figure  31.   Possible  superpositiors  of  Pake  doublet  and  paraEsagnetic 
spectrum  which  could  explain  ciata  at  550  and  7000  G. 
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Figure  32.   Calculated  (  )  and  observed  ( 


)  MT-IR  spectra  of 


protons  in  NH^Fe(S0/^)2.   Data  are  taken  from  Fig.  29  of 
ref.  [25]. 
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g(H)  =  --zr-7= >  -t1"m  ^-T  ^^-^'^^ 

where  A,  =-  3y^r   ,   This  ifi  the  same  form  as  Eq.  (2.3.17)  using  Kroon's 

1  L> 

notation,  and  corresponds  to  the  | 0>  ->  |+>  transition  in  Fig.  1.   Under 
conditions  applicable  to  Fig.  1,  some  of  the  magnetic  moments  are  paral- 
lel to  the  external  magnetic  field  and  some  are  antiparallei  depending 
on  their  states,  resulting  in  the  two  components  of  the  spectrum.   In 
the  case  of  a  paramagnetic  ion,  an  average  of  the  states  is  used  which 
has  a  fixed  direction  parallel  to  the  field,  resulting  in  only  one  com- 
ponent of  the  Pake  doublet. 

The  spectruni  in  Fig,  ?2  looks  similar  to  some  of  the  spectra  obtaine 
in  the  present  work,  especially  those  at  4000  G.   However,  it  is  deceptive 
as  the  asyiunetry  between  the  two  components  of  the  Pake  doublet  observed 
at  4000  G  causes  it  to  resenible  the  one.  asymmetric  component  in  Fig.  32. 
Thus,  a  direct  correspondence  between  the  present  data  and  Fig.  32  is 
ruled  out.   In  addition,  the  spectrum  observed  by  Kroon  [26]  is  reversed 
with  respect  to  the  present  data.   The  spectrum  in  Fig.  32  in  fact  re- 
sembles that  in  Fig.  30  which  was  calculated  by  a  different  method  as 
might  be  expected. 

We  now  compare  the  model  of  a  paramagnetic  line  superimposed  on  a 
classical  Pake  powder  pattern  doublet  to  the  observed  data. 

First,  peak-to-peak  separation  of  the  Paice  doublet  should  not  be 

affected;  however  superposition  of  two  peaks  close  together  might  cause 

an  apparent  shift  of  both  peaks.   The  observed  separation  decreased 

about  5%  as  H   increased  from  525  G  to  4000  G,  then  increases  about  10% 
o 

as  H   further  increases  to  9000  G.   This  is  not  inconsistent  with  the 
o 

model,  but  cannot  be-  even  qualitatively  predicted.   The  separation  is 
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also  affected  by  changes  in  broadening  of  the  individual  lines  in  the 
Pake  doublet. 

Next,  it  has  been  noted  the  low  frequency  wing  extends  at  a  rate  of 

-4 
about  4.6  X  10   G  per  G  increase  in  H  .   By  Eq.  (5.3.5),  the  width  of 

the  paramagnetic  line  is  given  by  A  ,  where 

A,  =  3p^r"^  -  u^H  /r\T    .  (5.3.6) 

1      D        Do 

-21 
In  the  present  case  we  calculate  p   ~  1.676  x  10    erg/G  taking 

V  --=  8.990  X  10^  Hz  and  H   =  17,768  G  from  ref.  [29]  since  y  =  vhl/H  . 

o  o 

The  proper  value  of  r  is  somewhat  in  doubt  as  mentioned  earlier,  but  is 

o 
close  to  3.5  A.   We  thus  calculate 

A,  =  1.13  X  lO"'^  H  (5.3.7) 

1  o 

approximately,  as  the  paramagnetic  line  width. 

Ttie  data  permit   loc.itj.on  of  only  the  lov;  frequency  extent  of  the 

paramagnetic  line,  which  extends  at  about  four  times  the  rate  in  Eq. 

(5.3.7)  calculated  t-'or  the  entire  line.   It  the  paramagnetic  line  as  a 

whole  shifts  with  changing  field,  the  line  width  may  change  at  a  rate 

—  h 

greater  or  less  than  4 , 6  x  10   G  per  G  increase  in  H  .   The  low  Ere- 

o 

quency  wing  observed  in  the  data  therefore  may  vel],  be  due  to  the  para- 
magnetic spectrum. 

The  high  frequency  wing  observed  at  low  H  does  not  fit  the  above 
theory,  as  Eq.  (5.3.7)  predicts  zero  paramagnetic  line  width  as  H 
approaches  zero;  howev(;r  th-^  data  indicate   about  5.5  G  line  width  if 
the  high  and  low  frequency  wings  are  counted  as  the  limics  of  the  para- 
magnetic spectrum. 

If  the  param.agnetic  spectrum  is  of  relatively  low  intensity  at  low 

H  ,  it  could  account  for  the  s"nai  I  peak  between  Z'aa   Pake  doublet  peaks. 
o 
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Anci  if  it  increases  in  relative  intensity  as  well  as  moving  to  lower 
frequencies,  as  H   increases,  it  could  explain  the  intensity  asymmetry 
of  the  Pake  doublet  peaks.   The  intensity  of  the  Pake  doublet  spectrum 
will  decrease  somewhat  due  to  broadening,  but  so  will  the  intensity  of 
the  paramagnetic  spectrum  which  widens  with  increasing  H^.   It  is  there- 
fore rather  unlikely  that  the  paramagnetic  spectrum  superimposed  on  the 
Pake  doublet  can  cause  much  intensity  asymmetry. 

Finally,  the  entire  Pake  spectrum  is  shifted  about  0.5  G.   Kroon 
[26]  does  predicc  a  component  of  the  local  field  due  to  the  paramagnetic 
ion  which  is  independent  of  orientation  of  H  and  which  would  therefore 
shift  the  entire  paramagnetic  spectrum.   It  is  given  by 

C--  -I   u^r^^P^Ccosej^)  (5.3.8) 

where  g  is  the  angle  between  the  z  axis  of  a  rectangular  coordinate 
system  in  the  crystal  and  a  vector  from  the  nucleus  being  observed  to 
the  dipole  or  paramagnetic  ion.   Averaging  over  the  twelve  nearest  neigh- 
bor positions  results  in  zero  net  shift,  some  parts  of  the  spectrum 
being  shifted  one  way  and  some  parts  the  other  way. 

In  order  for  the  Pake  doublet  to  be  shifted,  a  relatively  unifona 
field  would  have  to  exist  throughout  the  sample.   Possibly  a  very  unifom 
distribution  of  0  could  cause  such  a  field,  however  the  effect  should 

increase  vrith  H  which  was  not  the  case, 
o 

Fiott  and  McCormick  [25]  measured  the  shift  of  the  anomalous  peak 
from  the  center  of  the  line  as  s  function  of  temperature  and  compared  it 
with  the  calculated  shift.   At  ^.2   K,  the  shift  is  about  30  to  35  G,  mucl 
larger  than  that  predicted  by  Eq.  (5.3.7)  even  after  considering  that 
the  value  of  r  is  different.   However,  the  peak  in  Fig.  30  taken  from 
reference  [25]  i;s  at  aboul;  -3,  un.its  no^  ^iveu.   Units  are  also  om.itted 
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on  the  experimental  spectram,  and  not  much  information  is  furnished  with 
regard   to  temperature,  field  scrength  H  ,  etc.,  so  it  is  difficult  to 
draw  any  qu.antitative  conclusions  from  ref.  [25]. 

Line  width  variation  as  a  function  of  0   ccncentratio.i  is  shown  ia 
Fig.  29  and  is  consistent  with  expectations.   A  decreasing  0   concen- 
tration results  in  fewer   N  nuclei  close  enough  to  an  0   molecule  to 
be  influenced  by  it,  and  broadening  approaches  that  due  to  internolecular 
dipole-dipole  interactions  between   N  nuclei. 

The  main  feature  of  the  spectrum  at  low  0„    concentration  is  the 
asymmetry  of  the  peak  intensity.   Because  the  relaxation  time  was  on  the 
order  of  several  hours  it  was  difficult  to  obtain  good  resolution  of  Che 
line  shape ;  however ,  the  narrowirig  and  as^j^mmetry  were  reproducible. 

If  the  spectrum  is  a  superposition  of  a  paramagnetic  line  and  a 
Pake  doublet  as  discussed  earlier,  it  is  still  a  possible  mechanism  here, 
as  the  enhancement  due  to  a  relatively  faster  relaxation  rate  could 
offset  the  reduction  due  to  the  lower  concentration. 

In  summary,  the  paramagnetic  0   could  account  for  the  lov/  frequency 
wing  and  part  of  the  asymmetry  of  the  Pake  peaks.   Also,  the  Increasing 
broadening  with  increasing  H   and  with  increasing  0^  concentration  are 
to  be  expected.   It  do3S  rot  fully  explain  the  peak  asyrinietry ,  the  shift 
of  the  Pake  doublet,  and  the  behavior  of  the  high  frequency  wing. 


5.3.2   Ort>io-P ara  Conversion 

r 

As  hv/kT  is  on  the  order:  o'  10    in  the  present  experiment,  the 
|+>,  |0>  and  |0>,  !->  pojmlation  differences  are  equal  to  an  excellent 
approximation  according  :o  the  Boltkmfnn  distribution.   Therefore,  the 
Pake  doublet  shovm  in  Fig.  1  is  symmetric. 
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Hydrogen  molecules,  consisting  of  two  spin-'^  protons,  also  have 
a  triplet  and  singlet  state  siirdlar  to   M„  discussed  in  Sec.  2.2. 
Triplet  and  singlet  state  molecules  are  called  ortho  and  para  molecules, 
respectively.   UTien  hydrogen  is  cooled  at  a  nominal  rate  to  below  1.2  K, 
ortho-para  conversion  does  not  occur  rapidly  enough  to  maintain  thermal 
equilibrium.   Hardy  and  Berlinsky  [27]  and  Berlinsky  and  Hardy  [28]  have 
shown  that  the  conversion  rates  from  the  three  ortho  states  depend  on  the 
state  of  the  ortho  molecule,  the  magnetic  field  H  ,  and  the  angle  between 
the  molecular  axis  and  H  .   Thus,  the  |+>,  | 0>  and  J0>,  [->  population 
differences  are  not  equal  which  results  in  an  asymmetric  Pake  doublet. 

A  family  of  curves  is  shown  in  ref.  [28]  which  closely  resembles 
the  data  in  Figs.  26-28.   However,  the  asymmetry  increases  with  decreas- 
ing field  H  in  ref.  [28],  counter  to  the  field  dependence  in  the  present 
work. 

It  is  noted  that  the  theory  in  refs.  [27.28]  depends  on  a  net 
ortho-para  conversion.   The  spin  state  populations  are  set  at  high  tem- 
perature, and  an  energy  difference  between  the  ortho  and  para  states 
causes  this  net  conversion  at  low  temperature. 

In  the  case  of   N„,  the  crtho  and  para  states  are  separated  by  on 
the  order  of  a  microdegree  [30],  and  hance  o->-tho-para  conversion  v/ould 
not  be  significant  unless  the  sanple  were  cooled  to  about  1  \iK. 

The  theory  of  ortho-para  conversion  presenued  in  refy.  [27,23]  must 
therefore  be  discounted  on  both  experimental  and  theoretical  grounds  as 
applied  to  the  present  case. 


14 
5.3.3 N  Impurity 

■"A  14       15   14 

The  sample  containe«'  about  Oo5%   N  in  the  foiiii  of   N".,  and   N-  N 


mc 
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lecules.   In  the  mixed  molecule,  the   N  Zceman  energy  if  perturbed  by 


the   N- "  N  dlpole-dipole  interaction  which  is  dependent  on  the  state  of 
the   N  molecule.   The  perturbation  is  on  the  same  order  as  the   N„ 
dipole-dipole  interaction,  and  the  resulting  spectrum  overlaps  the  Pake 
pov;der  pattern  spectrum  due  to  ths   M„  molecules.   This  effect  has  been 
worked  out  theoretically  in  Sec.  2.5  and  experimental  results  are  pre- 
sented in  Chap.  VI. 

15   14 
Due  to  the  low  concentration  of   N-   M  molecules,  a  superposition 

of  the  mixed  nolecule  spectrum  and  the   N„  Pake  spectrum  cannot  differ 

significantly  from  the  pure  Pake  spectrum. 

The   M  nuclei  in  either   N   or   N-   N  mixed  molecules  undergo 

nuclear  quadrupole  transitions,  some  of  which  are  of  the  same  frequency 

as  the   N  Zeeman  transitions  for  certain  orientations  and  external 

magnetic  field  strengths.   The  transition  frequencies  are  shown  in  Fig. 

3.3  as  a  function  of  H   and  the  angle  between  H   and  the  m.oiecular  axis. 

o  o 

It  can  be  seen  that  only  a  small  percent  of  the  transitions  are 

near  the   N  Zeeman  transition  frequency  and  therefore,  considering  the 

14  14  .  .     .     .  c 

small  concentration  of   N  nuclei,  the   N  spectrum  m  Liie  vicinity  ot 

the   N  Zeeman  spectrum  must  be  undetectable. 

There  is  another  interaction,  however,  knowa  as  level  crossing 

which  could  have  a  much  greater  effect  on  the  spectrum  than  indicated  by 

14  r   1 

the  low  concentration  of   N  in  the  sample.   Smith  and  Coccs  l31J  have 

studied  this  effect  in  thiourea,  SC(^)H  )  ,  in  an  external  magnetic  field 
of  such  a  magnitude  as  to  make  the  proton  Zeeman  frequency  about  equal 
to  the  one  of  the  nitrogen  quadrupole  transition  frequencies.   In  this 
region  spin  exchange  can  occur  which  shifts  the  nitrogen  resonance  fre- 
quency.  The  effect  is  schematically  indicated  in  Fig.  34  which  has  been 
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10  HgCkG) 


Figure  33.  Quadrupole  transition  frequencies  of  ^-^N  and  Zeeinan  frequencj 
of  l^N  as  a  function  of  magnetic  field  K^  and  angle  9  betweer 
Hq  and  the  raolecular  axis.  The  dashed  lines  indicate  transi- 
tions which  are  forbidden  in  the  limit  of  low  or  high  magneti 
fields. 
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Figure  34.   Qualitative  effect  on  energy  levels  when  the  proton  Zeeman 

transition  energy  matches  a  -'-^N  quadrupolar  transition  energy 
in  SC(NH2)2-   The  +1,  -1,  and  0  levels  of  l^N  are  labeled  N^., 
and  Nn,  respectively,  and  the  proton  levels  are  labeled 


N_. 

+h   and  -Ij.   The  effect  is  illustrated  for  the  case  Z  ^  A,  but 

should  occur  at  Z  %  B  and  Z.  ^  C,    also.   The  diagram  is  a 

modification  of  Fig.  8  in  ref.  [31j. 
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exaggerated  for  clarity.   As  H   is  increased,  the  proton  Zeenian  energy 
Z  increases  until  it  approaches  the  energy  difference  A  between  the 
nitrogen  -Kl  and  -1  states,  N^  and  N_.   ^^fhen  Z  is  slightly  below  A,  the 
effect  is  to  'repel'  the  N   and  N_  levels  which  increases  A.   They 
observed  the  N  ,  N'   transition  B  and  N^ ,  N^  transition  C,  and  found  3 
decreased  while  C  increased.   As  H  was  further  increased  beyond  the 
value  where  A  and  Z  were  equal,  the  effect  reversed.   Similar  effects 
should  occur  in  the  region  where  the  proton  Zeeman  energy  approximately 
equals  the  N_,  N   or  N_|_,  N   transition.   An  expected  consequence  veri- 
fied experimentally  [3l]  is  that  level  crossing  permits  cross  relaxation 
which  greatly  enhances  the  relaxation  rates  of  both  the  nitrogen  and 

protons. 

Relative  to  the  present  experiment,  the  experiment  in  ref.  [31]  is 
much  simpler,  thus  facilitating  theoretical  analysis  and  experimental 
verification.   First,  a  single  c-.-ystal  was  available  which  permitted  a 
knowledge  and  choice  of  orientation  of  H  .   Next,  only  two  types  of  mole- 
cules exist  with  nitrogen  frequencies  separated  enough  to  permit  obser- 
vation of  only  one  type.   Thus,  all  the  nitrogen  nuclei  had  the  same 
three  energy  levels  and  level  crossing  occurred  at  well-defined  magnetic- 
fields,  permitting  observation  of  the  effect  at  values  of  H^  changing 
continuously  trom  below  to  above  the  crossing.   Third,  effects  were 
intramolecular  with  the  same  relative  oriencation  of  protons  and  nitro- 
gen.  Finally,  the  proton  Zeeman  resonance  was  a  single  line. 

In  contract,  the  sample  used  in  the  present  experiment  contained 
"'"\  ,  "''S   and  ■'^N-^^N  molecules.   There  are  three   N^  dipole-dipola 
perturbed  Zeeraan  levels,  three  "^N  quadrapole  levels  strongly  field 
dependent,  and  six   N-"' N  levels  with  transitions  in  the  vicinity  Oi- 
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the   N   transitions.   There  are  four  molecules  in  the  unit  cell,  each 
with  a  different  orientation,  and  each  has  12  nearest  neighbors  which 
must  be  considered.   Finally,  the  sample  was  polycrystalline  with  effec- 
tively a  continum  of  orientations,  and  a] 1  of  the  energy  levels  are 
orientational  dependent. 

The  twelve  levels  result  in  seven  transitions  at  very  low  external 
fields  and  perhaps  eight  transitions  at  higher  fields.   At  every  magnetic 
field  and  orientation,  a  minimum  of  five  transitions  are  within  about  a 

two-kHz  spread,  and  at  certain  values  of  H   and  orientations,  as  many  as 

o 

7  transitions  are  in  the  same  neighborhood.   Thus  the  possibilities  for 

14 

level  crossing  are  numerous.   Even  though  the   N  concentration  is  only 

on  the  order  of  0.5%,  cross  relaxation  can  cause  an  enhancement  of  the 
level-crossing  effect  by  increasing  the  relaxation  rate.  It  therefore 
could  have  a  significant  effect  on  the  line  shape. 

There  is  an  intuitive  objection  to  level  crossing  as  the  major 
mechanism  causing  the  anomalous  line  shape,  however.   The  observed 
effects  such  as  the  peak  intensity  as>T:nmetry,  the  high  and  low  frequency 
wings,  and  the  .-^mall  hump  between  the  Pake  doublet  peaks  all  change 
monotonically  with  Increasing  H  .   This  would  be  more  plausible  if  the 
various  transitions  were  all  either  approaching  tov/ards  or  receding  from 
each  other,  or  some  of  each.   Ho>.'Cver,  in  the  actual  sample,  transitions 
approach,  cross,  and  teciide  and  we  would  expect  shifts  first  one  way 
then  the  other  and  a  spect.rum  vit  h  a  irore  complex  magnetic  field  depen- 
dence. 

5.3.^   Conclusions 

The  spectruu;  consists  of  a  Pake  powder  pattern  doublet  with 
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broadening  and  a  low  frequency  wing  due  to  paramagnetic  0„.   The  asym- 
metry in  peak  intensity  has  not  been  explained.   Other  systems  such  as 

CO,  D       C^ and  0   impurity  in  a   N- ,    N^ ,    N-   N  mix  might  be 
helpful  in  identifying  the  various  components  of  the  observed  spectra. 


CHAPTER  VI 
RESULTS  AND  DISCUSSION  -  N^IR  SPECTRUM  OF 
l^N  IN  SOLID  l-'^N-l^^N 


6.1   Sample 

15 

The  nitrogen  sample  contained  a  statistical  mix  o!:  3J%   N  which 

resulted  in  proportions  of   N^,    N-   N,    N  molecules  in  a  1:4:4 
ratio.   In  the  vicinity  of  the   N  Zeeman  spectrum  there  were  twice  as 
many   'N  nuclei  participating  in  the   '  N-   N  mixed  molecule  spectrum  as 
in  the   N   (Pake  doublet)  spectrum.   We  therefore  calculate  the  two 

spectra  separately  then  superimpose  50%  of  the   N  Pake  doublet  ^spectrum 

15   14 
on  the   N-  N  mixed  spectrum. 


6.2   Data 

The  data  are  presented  in  Figs.  35  and  36.   A  relaxation  time  on 
the  order  of  an  hour  was  encountered  makin.^.  it  difficult  to  obtain  a 
good  spectrum  and  at  external  fields  below  5000  G  the  signal-to--noise 
ratio  was  very  poor.   The  relaxation  time  was  not  shortened  by  addition 
of  a  paramagnetic  impurity  because  of  the  effect  of  0^  on  the  line  shape. 

The  intensity  of  the  signal  is,  of  course,  expected  to  increase 

with  increasing  K   as  was  observed.   Absolute  intensity  measurements 
o 

were  not  made;  however  the  intensity  increase  was  more  pronounced  when 
the  field  was  increased  from  4000  G  to  5000  G    than  at  higher  or  lower 
fields.   At  lower  fields,  results  were  very  poor,  but  at  higher  fields 
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6000 


5000  G 


Figure  35.   N!-m  spectra  of  "N  in  a  N  enriched  to  33%   N  at  5000  and 
6000  G. 
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Figure  36.   :J>fH  spectra  of  ''"'n  in  a  N   enriched  to  33%   N  at  7000  and 
9000  G. 
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5000  G 


1  kHz 

K — H 


Av  =  110  Hz 


Av  =  60  Hz 


Figure  37.   Calculated  ?nill  spectra  of   'N  in  a  N2  enriched  to  33%  ■'-^N  at 
5000  G.   Ths  lower  figure  is  the  sura  of  the  triplet  (  •••  ) 
resulting  from  the  l-N-l'+N  molecules  and  50%  of  the  Pake 

double-.  (  )  resulting  from  the  ^-^N2  molecules.   Gaussian 

broadening  of  60  and  110  Hz  is  shown. 
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reproducibility  was  fairly  good.   Two  curves  at  6G00  G  and  7000  G  are 
shewn  as  an  illastraticn.   The  main  features  change  rather  uniformly 
with  increasing  applied  magnetic  field, 

6. 3 Computer  Simulation 

A  FOCAL  program  (Appendix  A)  was  used  to  calculate  the  powder  pat- 
tern spectrum  predicted  by  Eq.  (2.5.29)  broadened  by  a  Gaussian  function, 

Figure  37  shows  the  calculated  spectrum  at  5000  G.   The  lower  figure 

15  14 
shows  the  individual  components  obtained  from  the   N-  N  triplet  (dottec 

curve)  and  the   N  Pake  doublet  (dashed  curve).   The  solid  line  is  the 

sum  of  the  triplet  and  50%  of  the  Pake  doublet  intensity. 

Effects  of  varying  the  intermolecular  broadening  are  shown  by  the 

two  solid  curves  in  lig.  37.   Comparison  with  Fig.  35  indicates  the  line 

width  is  about  100  Hz.   Brookeman.  McEnnan,  and  Scott  [21]  calculated 

the  line  width  of  dipolar  broadened  NQR  lines  in  solid  a   nitrogen  and 

found  it  to  be  on  the  order  of  ! 00  Hz.   A  similar  broadening  in  the  case 

of   N  NMR  is  expected  although  the  calculations  are  not  the  same. 

6.4  Data  Analysis 

Figure  38  is  a  comparison  of  the  data  with  the  computer  simulations 
Intensity  differences  between  the  data  and  the  simulations  are  not  sig- 
nificant. 

The  theory  predicts  the  main  features  of  the  spectrum  such  as  the 
number  and  position  of  the  peaks.   However,  the  relative  intensities  are 
not  correctly  predicted.   The  asymmetry  of  the  Pake  doublet  discussed  in 
Chap.  V  and  the  level- crossing  mechanism  mentioned  in  Sec.  5.3.3  may  be 
causes  of  the  unexpected  peak  intensities. 
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6.4.1  Asymmetric  Pake  Doublet 

In  this  chapter  and  in  Sec.  2.5,  it  was  assumed  the  Pake  doublet 
would  be  symmetric.   Noting  the  data  in  Chap.  V  do  not  bear  this  out,  we 
would  expect  some  similar  effect  here.   The  narrow  but  asymmetric  peaks 
observed  irtth  very  low  0   concentrati.on  in  pure   N   (Fig.  29)  could 
well  account  for  the  two  peaks  observed  in  the  data  at  7900  G  and  9000  G 
in  Figs.  36  and  38.   The  simulation  has  two  peaks  also  but  one  is  of 
such  low  intensity  that  it  is  not  resolved.   .'Vn  analytical  formula  for 
the  asymmetry  in   N„  is  not  yet  known  so  it  cannot  be  included  in  the 
computer  calculations  at  this  Uiine.   Furthermore,  the  shift  in  the  entire 
Pake  spectrum  obseirved  with  0.1%  0„  added  to   N^  is  not  understood 
quantitatively.   Thus  the  experimental  spectrum  is  thought  to  be  the 
calculated  spectrum  modified  by  the  asyioDietric  Pake  doublet  with  perhaps 
other  effects  due  to  leve]  crossing  also  present, 

6.4.2  Level  Crossing 

Level  crossing  was  discusse.a  in  Sec.  j.3.3.   It  should  be  a  much 
more  significant  effect  in  the  present  case  with  337o   N  as  concentra- 
tions of  the  various  types  nuclei  are  of  the  saiue  order. 

It  was  noted  in  Chap.  V  that  the  level  crossings  varied  in  a  com- 
plex way  with  field,  orientation,  type  of  nuclei  and  uolecule,  etc.,  but 
the  spectrum  smoothly  went  from  a  nearly  sj-nroetric  Pake  doublet  to  a 
highly  asymmetric  doublet.   In  the  present  case  the  change  between  5000 
G  and  7000  G  is  not  as  smooth.   Furthermore,  it  was  noted  in  Sec.  6.1 
that  the  intensity  did  not  increase  smoothly  as  the  roagnetic  field  in- 
creased. 

We  note  that  level  crossing  increases  mere  rapidly  at  approximately 
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the  same  magnetic  fields.   Referring  to  Fig.  33,  there  is  an  increase  in 

14 
level  crossing  above  3/ SO  G  due  r.o  the  i  ,  0  and  1,-1   N  quadrnpole 

transitions  overlapping  for  some  orientations  and  another  increase  above 

15  14 

4700  G  where  the   N  Zoeman  line  overlaps  the  1,0   N  transition.   The 

15   14 
actual  situation  is  a  liLtle  more  complex  as    N-   N  dipole-dipole  per- 
turbation has  not  been  included  in  Fig.  33,   Thus,  level  crossing  may  be 

15 
more  effective  here  than  in  thc^  case  of  che    N  Pake  doublet  discussed 

in  Chap.  V. 


6.5   Conclusions 

The    N  MMR  spectrum  observed  in  the  33%   N  sample  is  a  corabi.na- 

tion  of  a  triplet  due  to  dipole-dipole  perturbation  of    N  in    N-   N 

15 

molecules  superimposed  on  the  Pake  doublet  resulting  from  the   N^ 

molecules.  Differences  between  tlie  calculated  and  observed  line  shapes 
are  urobably  due  to  the  asymmetry  of  the  N  Pake  doublet  spectrum  and 
to  level-crossing  effects. 

It  would  be  informative  to  examine  the  spectrum  as  a  function  of 

15  15   14 

N  concentration,  especially  if  the  proportion  of    N-   N  mixed  mole- 

15 
cules  to   N„  molecules  wera  increased  above  that  of  a  statistical  mix. 

Study  of  a  single  crystal  specimen  would  probably  provide  a  definitive 

explanation  of  the  complex  behavior  observed  in  the  powder  samples. 

However,  it  is  exceedingly  difficult  to  obtain  a  good  single  crystal  of 

a  nitrogen  because  of  trie  disruption  of  the  a-B  phase  transition. 


CHAPTER  VII 
RESULTS  AND  DISCUSSION  -  NUCLEAR  SPIN-L/TTICE  RELAXATION 
IN  LIQUID  AND  SOLID  ^^N  Al^D    ^^N 


15       14 
7.1   Previous   N   and   N„  RelaxatLoii  Time  Studies 


Relaxation  time  studies  in  nitrogen  are  of  interest  because  they 
afford  the  opportunity  to  isolate  two  different  relaxation  processes  in 

nearly  the  same  dynamic  environment.   Molecules  of  the  two  stable  iso-- 

15       14 
topes   N„  and   N„  have  only  a  small  isotopic  mass  difference  and 

therefore  nearly  the  same  molecular  dynamics,  but  the  dominant  relaxation 
mechanism  is  different  in  each  case.   Relaxation  via  dipolar  interactions 
and  anisotropic  chemical  shift  is  negligible  for  both  isotopes.   Spin- 
rotational  relaxation  is  therefore  dominant  in  "  N  which  has  no 

14 
quadrupolar  interaction,  whereas  quadrupolar  relaxation  in   N.,  dominates 

the  spin-rotational  relaxation  [24j.   Hence  relaxation  time  measurements 

15       14 
in   N_  and   N^  yield  T,    and  T, ^,  respectively. 
2        2        Isr      IQ     "^ 

14 
Relaxation  times  had  been  measured  in   N  from  4,2  K  to  145  K 

[22,24,32]  and  in  "'"^N  from  77  K  to  130  K  [24].   Krynicki,  Rahkamaa,  and 

Powles  [24]  found  that  the  data  tend  towards  the  Hubbard  relation  [14] 

T*  T*  =  1/6  (2.6.6) 

Q  sr 

as  the  temperature  is  reduced  to  77  K.   As  the  sample  temperature  is 

lowered  from  150  K,  t„  at  first  is  less  than  t   but  they  become  equal 

Q  sr        -         - 

at  about  77  K.   Thus  t  and  t   are  approaching  the  i'   >>  x        condition 

Q      sr  Q     •'sr 

which  is  necessary  to  the  Hubbard  relation  as  discussftd  in  Sec.  2.6. 
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The  present,  wcrk  extends  the  relaxation-time  measurements  in  liquid 

15 

^^N   from  77  K  dovm  to  the  triple  point  63  K,  and  in  S-solid   N^  from 

63  K  down  to  33  K,  near  the  a-B  phase  transition. 


7.2  Data 

The  observed  spin-lattice  relaxation  times  are  shown  in  Fig.  39  as 
solid  circles.   The  open  circle  at  38  K  represents  the  average  value  of 
55  seconds  observed.   It  is  possible  this  value  was  affected  by  the 
residual  0  impurity  and  such  a  long  T   is  also  more  susceptible  to 
magnet  instabilities,  therefore  it  is  not  used  in  the  analysis.   The 
data  are  in  excellent  agreement  at  77  K  with  those  of  ref.  [24]  which 
are  plotted  as  x's,  and  confirm   the  small  departure  from  linearity 
noticeable  in  the  data  from  ref.  [24].   The  solid  line  in  Fig.  39  is  a 
smooth  curve  through  the  data  and  the  dashed  line  is  a  smooth  curve 
through  ■'"'^N  data  from  refs.  [22,24,32].   It  should  be  noted  that  relaxa- 
tion times  for   N  and   N  are  in  milliseconds  and  seconds,  respective- 
ly, so  it  can  be  seen  that  the  spin-rotation  interaction  is  negligible 

14  .    -,  . 

in   N.,  as  asserted  earlier. 

Equations  (2.6.3)  and  (2.6.4)  were  used  to  calculate  the  correlatioi 
times  which  are  plotted  in  Fig.  40.   The  spin-rotational  constant  C  in 
Eq.  (2.6.3)  was  taken  as  22  kHz  [33]  and  the  most  recent  value  of  the 
quadrupole  coupling  constant  a  qQ/h  =  5.39  ±  C.05  KHz  [34]  was  used  ia 
Eq.  (2.6.4).   The  smooth  curves  and  solid  circles  v/ere  calculated  from 
the  smooth  curve?-  and  solid  circles,  respectively,  from  Fig.  39.   Final- 
ly, the  reduced  correlation  tines  were  calculated  from  the  solid  circles 
and  s;r.ooth  curves  in  Fig.  40  and  a  plot  of  t*  versus  t*^  on  a  log-log 
scale  is  shoOT  in  "i^,  41.   I'le  Hubbard  r(;l?tion  is  shov/n  as  a  dashed 
line  for  compari.scn. 
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Krynicki,  Rahkamaa,  and  Powles  [2A]  noted  an  approach  towards  the 

Hubbard  relation  as  they  lowered  the  temperature  ro    77  K  and  the  present 

data  confirm  this  tendency  as  the  temperature  is  lowered  to  the  triple 

point.   We  note  that  a  temperature  independent  constant  other  than  1/6 

in  the  product  t*  t*   would  result  in  a  line  parallel  to  the  dashed  line 
Q   sr 

but  offset  from  it.   The  experimental  curve  in  the  liquid  appears  to 

approach  the  dashed  curve  both  in  slope  and  in  value  of  th.e  constant. 

The  discontinuity  at  the  triple  point  is  due  to  the  discontinuity  in 

T   .   Any  discontinuity  in  T    which  inay  exist  within  the  limits  of 

experiaiental  accuracy  cannot  be  Itirge  enough  to  significantly  modify  the 

graph.   The  slope  of  t*  versus  t '^  should  "ary  depending  on  how  well  the 

condition  t*   <<  i*  is  met.   As  t*  decreases  significantly  across  the 
sr     Q  Q 

gap,  the  condition  r*   <<  x*   is    farther  from  being  realized  on  the  solid 

sr     Q 

side  of  the  gap,  which  -..'oul  d  predict  a  changs  in  the  slope.   A.ny  such 
discontinuity  in  the  slope  is  not  discernible  against  scatter  in  the 
experimental  results.   The  significan.;  observaticn  is  that  the  slope  is 
approximately  the  same  across  the  triple  point  gap,  and  at  lower  tempera- 
tures in  the  solid  the  product  of  the  reduced  correlation  times  ap- 
proaches even  closer  to  a  temperature  independent  constant  than  in  the 
liquid.   The  value  of  the  constant  is  apparently  less  than  1/6,  however. 


7.3   Crystalline  Structure  of  B-Solid  Nitrogen 

In  order  to  understand  the  behavior  of  data  above  and  below  the 
triple  point,  we  first  look  at  the  structure  of  g  nitrogen.   It  is  known 
from  x-ray  diffraction  [20, 35]  to  be  hexagonal  close  pc.cked  vith  two 
molecules  per  unit  cell  in  space  group  P6^/mmc(rt^,  )  .   The  c/?.  ratio  is 
very  close  [38]  to  uhe  ideal  close  packing  value  of  1,633.   Each  voc  Lecul<; 
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is  surrounded  by  12  nearest  neighbors  forming  a  cage  structure  as  showa 
in  Fig.  42.   The  anisotropic  part  of  the  intermolecular  potential  ex- 
hibits six  primary  minina  produced  by  groups  of  four  molecules  at  the 
corners  of  the  square  faces  forming  the  cage  structure.   Three  of  these 
minima  are  above  the  horizontal  midplane  of  the  figure  and  three  below. 
While  one  end  of  the  molecular  axis  points  towards  the  center  of  one  of 
these  square  faces,  the  other  end  intersects  a  triangular  face.   For  an 
ideal  c/a  ratio  the  polar  angle  0  between  the  c  axis  and  a  line  drawn 
from  the  central  molecule  to  the  center  of  one  of  the  square  faces  is 
54°44'. 

As  the  potential  wells  are  rather  shallow,  the  molecule  executes 
frequent  reorientational  jumps  aiaong  the  six  potential  minima.   Jumps 
between  adjacent  minima  while  maintaining  the  54*^44'  polar  angle  produce 
an  angular  reorientation  by  48^''ll'  .   In  addition,  reorientation  between 
minima  may  occur  by  tilting  the  molecular  axis  through  a  polar  angle  by 
70"32'.   There  is  thus  a  twelve-fold  model  for  reorientation;  however, 
for  homonuclear  moiecuies.  ttic  six  orientations  following  a  tilt  are 
equivalent  to  the  six  without  a  tilt.   Tb.e  4S'-*il'  jumps  are  presumably 
more  probable  than  the  70^32'  -jiicps.   From  the  standpoint  of  an  NMR  ex- 
periment, rapid  reorientatioi?  amo^.g,  these  12  configurations  is  indis- 
tinguishable from  a  continuous  precession  about  the  c  axis  at  a  fixed 
polar  angle  6.   In  particular,  the  electric  field  gradient  tensor  is 
averaged  by  such  motion  accordtr.g  to  the  factor  <P_(cose)>  which  vanishes 

for  the  'magic'  angle  of  54^44*, 

14 

This  picture  is  confirmed  by  the  fact  that  the   N  nuclear  quadru- 

pole  coupling  constant  in  6  nitrogen  is  reduced  by  three  orders  of 
magnitude,  from  about  5  iffiz  to  3  kHz  [22].   The  small  quadrupolar 
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Fig-u're  42.  Crystal  structure  of  p.  nitrogen.  Open  circles  defining  the 
cage  structure  denote  hep  positions  of  molecules.  Only  the 
centrai  raolecule  is  illustrated.  All  molecules  precess  about 
the  c  axis  at  an  angle  9  -  54°  44'.  At  this  orientation  the 
molecular  axis  during  the  course  of  precession  passes  approx: 
iTiately  through  the  center  of  each  of  the  square  faces  of  the 
.surrounding  cage. 
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interaction  in  S  nitrogen  is  consistent  with  a  deviation  of  only  a  few 
minutes  of  arc  from  the  54*^44'  value,  and  has  been  shown  [34]  to  be 
correlated  with  the  variation  of  the  c/a  ratio  during  thermal  expansion. 

7.4  Data  Analysis 

We  now  can  provide  an  insight  into  the  behavior  of  t   at  the  triple 
point.   The  molecules  in  the  B  solid  are  relatively  free  to  reorient 
through  large  angles,  but  are  restricted  in  translational  freedom.   Fix- 
ing the  molecules  on  lattice  sites  ensures  each  one  has  its  own  space  in 
which  to  reorient  wich  little  oostruction  from  Che  other  molecules. 
Upon  melting,  the  collisions  due  to  the  unleashed  translational  degrees 
of  freedom  destroy  this  cooperativp  behavior  and  the  molecule  is  more 
hindered  in  its  rotation  even  thoup^h  the  average  d.msity  is  less  j,n  the 
liquid  phase,  and  t   increases  accordingly. 

Several  recent  ir.cleci.lar  c■.yn?.^ics  ccmputf^r  simulations  are  consis- 
tent with  this  picture  [37,33,39]    In  particular  time  traces  of  the 
orientation  of  a  molecule  in  the  liquid  near  the  triple  point  reveal  [39] 
that  molecules  tend  to  undergo  several  small-angle  reorientations  before 
reorientating  through  an  angle  on  the  order  of  one  radian. 

That  there  should  be  at  most  a  snail  discoiitlnuity  in  t   across 
the  triple  point  is  not  surprising,  considering  the  relative  freedom 
enjoyed  by  molecules  in  3  nitrogen  which  is  a  plastic  crystal.   Molecules 
can  readily  undergo  changes  in  angular  monontura  whether  or  nor.  restricted 
in  translational  freedom.   The  more  unexpected  result  is  that  at  lower 
temperatures  in  the  solid  the  product  of  reduced  correlation  times 
approaches  a  temperature,  indepexident  constant  even  more  closely  than  in 
the  liquid.   The  value  of  the  constanc  though  is  apparently  less  than  1/6.. 
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The  Hubbard  relation  was  based  on  an  isotropic  rotational  diffusion 
model  which  at  first  thought  would  not.  seem  tc  be  valid  in  the  case  of 
the  large  angle  reorientational  jurps.   The  theoretical  problem  has  not 
yet  been  solved.   However,  Ivanov  [40]  has  discussed  the  theory  of  rota- 
tional Brownian  movement  as  a  randooi-walk  process  including  the  question 
of  how  large  the  reorientational  jumps  may  become  before  their  discrete 
nature  is  distinguishable  from  a  truly  isotropic  rotational  diffusion. 
He  introduces  a  parameter  c  which  is  related  to  the  ratio  of  the  Raman 
and  infrared  line  widths  which  varies  from  3  in  the  case  of  isotropic 
rotational  diffusion  to  unity  when  the  discrete  nature  of  the  jumps  is 

dominant. 

For  the  angles  48°11'  and  70°32' ,  the  value  of  Ivanov's  parameter  q 
is  2.6  and  2.2,  respectively,  and  hence  rather  close  to  the  rotational 
diffusion  lindt,  especially  for  the  surely  more  predominant  smaller 
jumps.   The  general  conclusions  of  Ivanov's  analysis  seem  to  apply  here 
and  it  is  therefore  reasonable  to  expect  a  temperature-independent  slope 

fo^  T*  vs.  T*   in  g  nitrogen,  although  the  constant  may  not  be  1/6. 
Q      sr 

7.5   Conclusions 

We  have  observed  that  the  Hubbard  relation  t*  t *^  =  1/6  is  approach 
ed  even  more  closely  in  liquid  nitrogen  as  the  temperature  is  reduced 
from  77  K  to  the  triple  point.   This  contirms  the  effect  noted  by 
Krynicki  et  al.  [24]  at  higher  cemperatures. 

We  next  noted  that  the  decrease  in  t.,  upon  solidification  at  the 
triple  point  and  consequently  the  fact  that  the  product  of  the  reduced 
correlation  times  is  less  than  1/6  were  explained  by  the  j^reater  reori- 
entational freedom  in  3  solid  than  in  liquid. 
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Finally,  that  the  Hubbani  relat5.ca  should  be  obeyed  at  all  was 
explained  by  noting  the  rotational  diffusion  model  is  quite  good  even  at 
the  relatively  large  angle  reorientations  encountered  in  B  nitrogen. 
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APPENDIX  A 

FOC.\L  PROGRAM  TO  CALCULATE  THE  POUT)ER  PATTERM  SPECTRUM 

OF  15j^  i^j  a  N2  ENRICHED  TO  33%  ^^N 


C-SK    FOCAL    31969 

C^l.05    A    ••RU>J    VUM^35:h"HI;A    "JIDrHCHZJ-Vl.-S    V  =  l'l/9.?»7 
01.10     S    PI  =  3.1-'i159;A5K    "H"H;S    L  =  H/756fJ 

p,i.30    T    "  ••,"FIMISHED'S  •  ;QUIT 
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-^p.!'?)   s  Y=rsr'3fx);s  z=FCOS(X) 

C5P.P'^  S  C=(-A"t/P)/FSQT(A1  T3/27);S    n  =  2*KSQl{  A  I /3) 

0P.P5  S  F=FATr>3CFSQTCl-Ct2)/C)/3 

n2.27  I  (Al^P.3;S    E=E-2*PI/3 

f^2.3;^  F  1  =  1  »1*3;D    4.P1                                                                                          -^ 

.:iP./i'^  S  D=l  .•'552')67*C  l-3*Z;t2)*61  .29  +  256;D    '!».6 

i,12./4-3    S    D=---L'+5i2;G0T0    A. 6 

nA.^5  S  T=-FSG>J(An)*R*FCOS(V:) 

(Vj.l.'i  5  Hl  =  Cl/2  +  L*Z  +  T)/CL*Y/l  .Al/4214) 

^4.20  S  H2=(l/^+L*Z+T)/( l/2-L*Z+T> 

a/1.31  S  K1=FS0TC1/(1+H1  t2  +  H2»2)) 

r^H,^'^    S    K2=-Hl*Ki;S    K3  =  H2*K1 

04.5:^    S    n=-fil.29+(2-CKlt2-K3t25*Z~K2*(Kl+K3)*Y*1.4  1^2lA)-25^ 

0/4.55  T  ^P.O'^*"       "."I".  I 

C^/!.57  T  7,7. -!6."  -VT-'T,"       "  ,  "W  (  D  ♦•256  )  /  1  22«  53,  f 

n/1.6^  F  ^J  =  5!  2  +  D"/<*V»512+D-k'J*U;D    5.0 

^4.70  S  E=E<-2*PI/3 

PSS.ifl    5    R=FEXPC-C(M-512-.D5/V)tP/?^*Y*2ri03;D    7.f» 

ri7.in    5    R  =  H  +  FN)E'*KM> 
r-t7.2'^    S    s  =  F'\JEWcr>}*R) 
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